【BZOJ3331】[BeiJing2013]压力
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
1 2
1 3
2 3
1 4
4 2
4 3
Sample Output
2
1
1
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
题解：显然先用Tarjan求缩块。。。怎么求呢。。。基本功不扎实又去学了一发。
最后我们会得到一个树形结构，但是。。。怎么得到呢。。。其实对于每个块新建一个点连向块中的所有点即可。
然后就是一个类似于树的东西了，怎么统计树上有哪些路径必经一个点呢？差分即可。
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int maxn=200010;
int n,m,q,top,tot,sum,cnt;
int sta[maxn],low[maxn],HEAD[maxn],NEXT[maxn<<1],TO[maxn<<1],head[maxn<<1],next[maxn<<2],to[maxn<<2];
int s[maxn<<1],fa[19][maxn<<1],Log[maxn<<1],dep[maxn<<1],Q[maxn<<1];
inline void ADD(int a,int b)
{
TO[cnt]=b,NEXT[cnt]=HEAD[a],HEAD[a]=cnt++;
}
inline void add(int a,int b)
{
to[cnt]=b,next[cnt]=head[a],head[a]=cnt++;
}
void tarjan(int x)
{
dep[x]=low[x]=++tot,sta[++top]=x;
for(int y,i=HEAD[x],t;i!=-1;i=NEXT[i])
{
y=TO[i];
if(!dep[y])
{
tarjan(y),low[x]=min(low[x],low[y]);
if(low[y]>=dep[x])
{
sum++;
do
{
t=sta[top--],add(sum,t),add(t,sum);
}while(t!=y);
add(sum,x),add(x,sum);
}
}
else low[x]=min(low[x],dep[y]);
}
}
void dfs(int x)
{
Q[++Q[0]]=x;
for(int i=head[x];i!=-1;i=next[i]) if(to[i]!=fa[0][x]) fa[0][to[i]]=x,dep[to[i]]=dep[x]+1,dfs(to[i]);
}
inline int lca(int a,int b)
{
int i;
if(dep[a]<dep[b]) swap(a,b);
for(i=Log[dep[a]-dep[b]];i>=0;i--) if(dep[fa[i][a]]>=dep[b]) a=fa[i][a];
if(a==b) return a;
for(i=Log[dep[a]];i>=0;i--) if(fa[i][a]!=fa[i][b]) a=fa[i][a],b=fa[i][b];
return fa[0][a];
}
inline int rd()
{
int ret=0,f=1; char gc=getchar();
while(gc<'0'||gc>'9') {if(gc=='-')f=-f; gc=getchar();}
while(gc>='0'&&gc<='9') ret=ret*10+gc-'0',gc=getchar();
return ret*f;
}
int main()
{
//freopen("bz3331.in","r",stdin);
n=rd(),m=rd(),q=rd(),sum=n;
memset(head,-1,sizeof(head)),memset(HEAD,-1,sizeof(HEAD));
int i,j,a,b,c;
for(i=1;i<=m;i++) a=rd(),b=rd(),ADD(a,b),ADD(b,a);
cnt=0,tarjan(1),dep[1]=1,dfs(1);
for(i=2;i<=sum;i++) Log[i]=Log[i>>1]+1;
for(j=1;(1<<j)<=sum;j++) for(i=1;i<=sum;i++) fa[j][i]=fa[j-1][fa[j-1][i]];
for(i=1;i<=q;i++)
{
a=rd(),b=rd(),c=lca(a,b);
s[a]++,s[b]++,s[c]--,s[fa[0][c]]--;
}
for(i=sum;i;i--) a=Q[i],s[fa[0][a]]+=s[a];
for(i=1;i<=n;i++) printf("%d\n",s[i]);
return 0;
}