BZOJ 1123: [POI2008]BLO 图的连通性

1123: [POI2008]BLO

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 1319 Solved: 598
[Submit][Status][Discuss]

Description

Byteotia城市有n个 towns m条双向roads. 每条 road 连接 两个不同的 towns ,没有重复的road. 所有towns连通。

Input

输入n<=100000 m<=500000及m条边

Output

输出n个数，代表如果把第i个点去掉，将有多少对点不能互通。

Sample Input

5 5

1 2

2 3

1 3

3 4

4 5

Sample Output

8

8

16

14

8

HINT

Source

题解：
注意，这道题所想的数对是有序的，即（1，2）和（2，1）算两对，以及删掉的点与其他点构成的数对也算在答案内。

不难发现，如果一个点不是割点，那么删掉它之后并不会有不连通的点，所以它的贡献就只是（n-1）*2。
考虑割点的情况。如果一个点是割点，删掉它之后得到了若干个联通块，不能互相到达的点对数目就是这些联通块数目两两的乘积。
在tarjan找割点的时候就把答案统计出来，感觉还是挺巧的。
记联通块数目为x,每个联通块的大小为A1，A2，A3，…，Ax。
答案的一半为A1*(A2+A3+…+Ax)+A2*(A3+A4+…+Ax)+…+Ax-1*Ax 所以我们dfs时就统计每个子树的大小，得到siz数组，找到一个割点u之后，如果割掉u，v及它的子树和其他部分不连通，那么可以用siz[v]*size，size表示割点u之后，其他不连通的子树的大小。
最后用n-1-size表示非u的子树部分，可以看做一个联通块A1，此时的size就是A2+A3+…+Ax.
这样处理完，一定可以得到所有的答案。

[其实ans数组不只记录了割点的答案，非割点也是答案，所以可以一起输出]
注：如果输出是lld，就一定要转化为long long。之前因为这个一直WA。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;

inline int read(){
    int x=0,f=1;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;
}

const int N = 100000 + 10;
const int M = 500000 + 10;

int n,m;
struct node{
    int pre,v;
}e[M<<1];

int head[N],num=0;
void addedge(int from,int to){
    e[++num].pre=head[from],e[num].v=to;
    head[from]=num;
}

int indx=0,low[N],dfn[N];
int siz[N],size=0;
ll ans[N];
bool iscut[N];
void tarjan(int u,int f){
    low[u]=dfn[u]=++indx;
    siz[u]=1;ll size=0;
    int son=0;
    for(int i=head[u];i;i=e[i].pre){
        int v=e[i].v;
        if(v==f) continue;
        if(!dfn[v]) {
            ++son;
            tarjan(v,u);
            siz[u]+=siz[v];
            low[u]=min(low[u],low[v]);
            if(low[v]>=dfn[u]) {
                iscut[u]=true;
                ans[u]+=1LL*siz[v]*size;
                size+=siz[v];
            }
        }
        else low[u]=min(low[u],dfn[v]);
    }
    if(f==0&&son==1) iscut[u]=0;
    ans[u]+=1LL*size*(n-size-1);
}

int main(){
    n=read(),m=read();
    for(register int i=1;i<=m;++i){
        int u=read(),v=read();
        addedge(u,v);addedge(v,u);
    }
    tarjan(1,0);
    for(int i=1;i<=n;++i){
        if(iscut[i]){
           ans[i]=1LL*(ans[i]+n-1)*2;
           printf("%lld\n",ans[i]);
        }
        else printf("%lld\n",1LL*(n-1)*2);
    }
    return 0;
}