BZOJ3331 压力【必经点，圆方树】

题目描述

一个连通无向图，有Q个从$s_i$走到$t_i$的运输，问1~n每个点作为必经点的次数。
N≤100000,M,Q≤200000

题目分析：

必经点就是割点和起点终点，缩点双（新建一个点，点双的所有点都连在这个点上），重新建图构成一棵树，然后树上差分。

其实就是圆方树模板题。

但是我第一次打点双震惊了，退栈的终止条件不能再像边双那样用stk[top]!=u来判断，而应该用stk[top]!=v，在点双里面u,v在栈中是不一定连续的。就这个调了我1.5h。。。

Code：

#include<bits/stdc++.h>
#define maxn 200005
#define maxm 400005
using namespace std;
int n,m,Q,dfn[maxn],low[maxn],stk[maxn],top,tim,sz;
int ans[maxn],f[maxn],fa[maxn]; bool vis[maxn];
int fir[maxn],nxt[maxm],to[maxm],tot;
vector<int>G[maxn],Pf[maxn];
inline void line(int x,int y){nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y;}
void tarjan(int u){
	dfn[u]=low[u]=++tim,stk[++top]=u;
	for(int i=fir[u],v;i;i=nxt[i]){
		if(!dfn[v=to[i]]){
			tarjan(v),low[u]=min(low[u],low[v]);
			if(low[v]>=dfn[u]){
				G[++sz].push_back(u),G[u].push_back(sz);
				do G[sz].push_back(stk[top]),G[stk[top]].push_back(sz); while(stk[top--]!=v);
			}
		}
		else low[u]=min(low[u],dfn[v]);
	}
}
int find(int x){return !f[x]?x:f[x]=find(f[x]);}
void dfs(int u,int ff){
	fa[u]=ff,vis[u]=1;
	for(int v: Pf[u]) if(vis[v]) ans[find(v)]--,ans[fa[find(v)]]--;
	for(int v: G[u]) if(v!=ff) dfs(v,u),ans[u]+=ans[v];
	f[u]=ff;
}
int main()
{
	int x,y;
	scanf("%d%d%d",&n,&m,&Q),sz=n;
	while(m--) scanf("%d%d",&x,&y),line(x,y),line(y,x);
	tarjan(1);
	while(Q--) scanf("%d%d",&x,&y),ans[x]++,ans[y]++,Pf[x].push_back(y),Pf[y].push_back(x);
	dfs(1,0);
	for(int i=1;i<=n;i++) printf("%d\n",ans[i]);
}