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<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<vector>
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]<dep[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;j<dcc[i].size();j++)
{
if(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;
}