[bzoj3331][强联通][树上差分]压力

Description

如今,路由器和交换机构建起了互联网的骨架。处在互联网的骨干位置的
核心路由器典型的要处理100Gbit/s的网络流量。他们每天都生活在巨大的压力 之下。
小强建立了一个模型。这世界上有N个网络设备,他们之间有M个双向的 链接。这个世界是连通的。在一段时间里,有Q个数据包要从一个网络设备发
送到另一个网络设备。 一个网络设备承受的压力有多大呢?很显然,这取决于Q个数据包各自走
的路径。不过,某些数据包无论走什么路径都不可避免的要通过某些网络设 备。 你要计算:对每个网络设备,必须通过(包括起点、终点)他的数据包有
多少个?

Input

第一行包含3个由空格隔开的正整数N,M,Q。 接下来M行,每行两个整数u,v,表示第u个网络设备(从1开始编号)和
第v个网络设备之间有一个链接。u不会等于v。两个网络设备之间可能有多个 链接。
接下来Q行,每行两个整数p,q,表示第p个网络设备向第q个网络设备发 送了一个数据包。p不会等于q。

Output

输出N行,每行1个整数,表示必须通过某个网络设备的数据包的数量。

Sample Input

4 4 2

1 2

1 3

2 3

1 4

4 2

4 3

Sample Output

2

1

1

2

HINT

【样例解释】

设备1、2、3之间两两有链接,4只和1有链接。4想向2和3各发送一个数据

包。显然,这两个数据包必须要经过它的起点、终点和1。

【数据规模和约定】

对于40%的数据,N,M,Q≤2000

对于60%的数据,N,M,Q≤40000

对于100%的数据,N≤100000,M,Q≤200000

题解

做完lyd的书这题简直就是傻逼题…
点双+树上倍增+树上差分没了..
模板打错没谁了

#include
#include
#include
#include
#include
#include
using namespace std;
inline int read()
{
    int f=1,x=0;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
inline void write(int x)
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}
vector<int> dcc[100005];
struct node{int x,y,next;}a[400005];int len,last[200005];
void ins(int x,int y){len++;a[len].x=x;a[len].y=y;a[len].next=last[x];last[x]=len;}
int dfn[200005],low[200005],sta[200005],id,cnt,tp,root;
bool cut[200005];
void tarjan(int x)
{
    low[x]=dfn[x]=++id;sta[++tp]=x;
    int flag=0;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(dfn[y]==-1)
        {
            tarjan(y);low[x]=min(low[x],low[y]);
            if(low[y]>=dfn[x])
            {
                flag++;
                if(root!=x || flag>1)cut[x]=true;
                cnt++;int i;
                do
                {
                    i=sta[tp--];
                    dcc[cnt].push_back(i);
                }while(i!=y);
                dcc[cnt].push_back(x);
            }
        }
        else low[x]=min(low[x],dfn[y]);
    }
}
int fa[200005][25],dep[200005],bin[25];
void pre_tree_node(int x)
{
    for(int i=1;bin[i]<=dep[x];i++)fa[x][i]=fa[fa[x][i-1]][i-1];
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(y!=fa[x][0])
        {
            fa[y][0]=x;dep[y]=dep[x]+1;
            pre_tree_node(y);
        }
    }
}
inline int lca(int x,int y)
{
    if(dep[x]y])swap(x,y);
    for(int i=17;i>=0;i--)if(bin[i]<=dep[x]&&dep[fa[x][i]]>=dep[y])x=fa[x][i];
    if(x==y)return x;
    for(int i=17;i>=0;i--)if(bin[i]<=dep[x]&&fa[x][i]!=fa[y][i])x=fa[x][i],y=fa[y][i];
    return fa[x][0];
}
int sum[200005];
void meg(int x,int fa)
{
    for(int k=last[x];k;k=a[k].next)
        if(a[k].y!=fa)
        {
            int y=a[k].y;
            meg(y,x);sum[x]+=sum[y];
        }
}
int n,m,T,belong[200005],ans[200005];
int main()
{
//  freopen("a.in","r",stdin);
//  freopen("a.out","w",stdout);
    bin[0]=1;for(int i=1;i<=17;i++)bin[i]=bin[i-1]<<1;
    n=read();m=read();T=read();
    for(int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        ins(x,y);ins(y,x);
    }
    memset(dfn,-1,sizeof(dfn));
    for(int i=1;i<=n;i++)if(dfn[i]==-1)root=i,tarjan(i);
    int gg=cnt;
    for(int i=1;i<=n;i++)if(cut[i])belong[i]=++gg;
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=cnt;i++)
        for(int j=0;jif(cut[dcc[i][j]])ins(belong[dcc[i][j]],i),ins(i,belong[dcc[i][j]]);
            else belong[dcc[i][j]]=i;
        }
    pre_tree_node(1);
//  for(int i=1;i<=len;i++)printf("%d %d\n",a[i].x,a[i].y);
//  for(int i=1;i<=n;i++)printf("%d ",belong[i]);
    while(T--)
    {
        int x=read(),y=read();
        if(belong[x]==belong[y])
        {
            if(x==y)ans[x]++;
            else ans[x]++,ans[y]++;
            continue;
        }
        int LA=lca(belong[x],belong[y]);
        sum[belong[x]]++;sum[belong[y]]++;
        ans[x]++;ans[y]++;
        if(cut[x])ans[x]--;if(cut[y])ans[y]--;
        sum[LA]--;sum[fa[LA][0]]--;
    }
    meg(1,0);
    for(int i=1;i<=n;i++)
    {
        if(cut[i])write(sum[belong[i]]+ans[i]);
        else write(ans[i]);
        puts("");
    }
    return 0;
}

你可能感兴趣的:(bzoj,树上差分,强联通)