【BZOJ4435】[Cerc2015]Juice Junctions
Description
你被雇佣升级一个旧果汁加工厂的橙汁运输系统。系统有管道和节点构成。每条管道都是双向的，且每条管道的流量都是1升每秒。管道可能连接节点，每个节点最多可以连接3条管道。节点的流量是无限的。节点用整数1到n来表示。在升级系统之前，你需要对现有系统进行分析。对于两个不同节点s和t，s-t的流量被定义为：当s为源点，t为汇点，从s能流向t的最大流量。以下面的第一组样例数据为例，1-6的流量为3，1-2的流量为2。计算每一对满足a<b的节点a-b的流量的和。
Input
第一行包括2个整数n和m(2<=n<=3000,0<=m<=4500)——节点数和管道数。
接下来m行，每行包括两个相异整数a,b(1<=a,b<=n)，表示一条管道连接节点a,b。
每个节点最多连接3条管道，每对节点最多被一条管道连接。
Output
输出一个整数——每对满足a<b的节点a-b的流量之和。
Sample Input
6 8
1 3
2 3
4 1
5 6
2 6
5 1
6 4
5 3
1 3
2 3
4 1
5 6
2 6
5 1
6 4
5 3
Sample Output
36
题解：这题还真是神啊~
考虑，每一个点对对答案的贡献只能为1,2,3，我们分别讨论这三种情况。
贡献为1或2或3：用并查集看一下是否联通就行了。
贡献为2或3：跑完Tarjan后，它们一定在同一个边双里。
贡献为3：这个的做法就比较奇葩了，我们任意删除一条边，如果无论删除那条边它们都在一个边双里，那么贡献为3。这个怎么判断呢？先枚举删除那条边，将每一种情况下i所在的边双的编号hash起来，然后比较两个hash是否相等就行了。
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int maxn=3010;
typedef unsigned long long ull;
ull hs[maxn];
int n,m,ans,cnt,tot,top,sum;
int to[9010],next[9010],head[maxn],dep[maxn],low[maxn],ins[maxn],sta[maxn],bel[maxn],f[maxn];
void add(int a,int b)
{
to[cnt]=b,next[cnt]=head[a],head[a]=cnt++;
}
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;
}
void tarjan(int x,int fa,int ban)
{
dep[x]=low[x]=++tot,ins[x]=1,sta[++top]=x;
for(int i=head[x];i!=-1;i=next[i])
{
if(i==fa||i==ban||i==(ban^1)) continue;
if(!dep[to[i]]) tarjan(to[i],i^1,ban),low[x]=min(low[x],low[to[i]]);
else low[x]=min(low[x],dep[to[i]]);
}
if(dep[x]==low[x])
{
sum++;
int t;
do
{
t=sta[top--],ins[t]=0,bel[t]=sum;
}while(t!=x);
}
}
int find(int x)
{
return (f[x]==x)?x:(f[x]=find(f[x]));
}
int main()
{
n=rd(),m=rd();
int i,j,a,b;
for(i=1;i<=n;i++) f[i]=i;
memset(head,-1,sizeof(head));
for(i=1;i<=m;i++)
{
a=rd(),b=rd(),add(a,b),add(b,a);
if(find(a)!=find(b)) f[f[a]]=f[b];
}
memset(dep,0,sizeof(dep));
memset(low,0,sizeof(low));
tot=top=sum=0;
for(i=1;i<=n;i++) if(!dep[i]) tarjan(i,-1,-1);
for(i=1;i<=n;i++) for(j=i+1;j<=n;j++)
{
if(find(i)==find(j)) ans++;
if(bel[i]==bel[j]) ans++;
}
for(i=0;i<cnt;i+=2)
{
memset(dep,0,sizeof(dep));
memset(low,0,sizeof(low));
tot=top=sum=0;
for(j=1;j<=n;j++) if(!dep[j]) tarjan(j,-1,i);
for(j=1;j<=n;j++) hs[j]=hs[j]*233+bel[j];
}
for(i=1;i<=n;i++) for(j=i+1;j<=n;j++) if(hs[i]==hs[j]) ans++;
printf("%d",ans);
return 0;
}