圆方树模板

P4630 [APIO2018] Duathlon 铁人两项

纯圆方树的模板，只要会建图，就会做

在求割点（或桥）的基础上 加一个缩点的操作就行，需要注意的是，一个割点可以同时属于多个点双，所以对割点要多次连边

#include <bits/stdc++.h>
#define inf 0x7fffffff
#define ll long long
#define int long long 
//#define double long double
#define re register int
#define void inline void
#define eps 1e-8
//#define mod 1e9+7
#define ls(p) p<<1
#define rs(p) p<<1|1
#define pi acos(-1.0)
#define pb push_back
#define P pair < int , int >
#define mk make_pair
using namespace std;
const int mod=19940417;
const int M=1e9;
const int N=4e5+5;//?????????? 4e8
struct node
{
	int ver,next;
}e[N];
int head[N],tot=1;
int cnt,dfn[N],low[N],sz[N],w[N],sum,num;
int instack[N],top;
vector < int > g[N];
ll ans;
int n,m;
void add(int x,int y)
{
	e[++tot].ver=y;
	e[tot].next=head[x];
	head[x]=tot;
}
void addedge(int x,int y)
{
	add(x,y);add(y,x); 
}
void tarjan(int x,int in_edge)
{
	dfn[x]=low[x]=++sum;
	instack[++top]=x;
	num++;
	for(re i=head[x];i;i=e[i].next)
	{
		int y=e[i].ver;
		if(!dfn[y])
		{
			tarjan(y,i);
			low[x]=min(low[x],low[y]);
			if(low[y]>=dfn[x])
			{
				cnt++;
				int z;
				do
				{
					z=instack[top--];
					w[cnt]++;
					g[z].pb(cnt);g[cnt].pb(z);
				}while(z!=y);
				w[cnt]++;
				g[x].pb(cnt);g[cnt].pb(x);
			} 
		}
		else if(i!=(in_edge^1))  low[x]=min(low[x],dfn[y]);
	}
}
void dfs(int x,int pre) 
{
//	cout<<x<<endl;
	sz[x]=(x<=n);
	for(re i=0;i<g[x].size();i++)
	{
		int y=g[x][i];
		if(y==pre)  continue;
		dfs(y,x);
		ans+=2*sz[x]*sz[y]*w[x];
		sz[x]+=sz[y];
	}
	ans+=2*sz[x]*(num-sz[x])*w[x];
//	if(ans!=0)  cout<<sz[x]<<endl; 
//	cout<<
}
void solve()
{
	cin>>n>>m;
	cnt=n;
	for(re i=1;i<=n;i++)  w[i]=-1;
	for(re i=1;i<=m;i++)
	{
		int x,y;
		scanf("%lld%lld",&x,&y);
		addedge(x,y);
	}
	for(re i=1;i<=n;i++)  if(!dfn[i])
	{
		num=0;
		tarjan(i,0);
		dfs(i,i);
	}
//	for(re i=1;i<=cnt;i++)  printf("%d ",sz[i]);puts("");
	cout<<ans<<endl;
}
signed main()
{
//	freopen("Ain.txt", "r", stdin);
//	freopen("Aout.txt", "w", stdout);
    int T=1;
//    cin>>T;
    for(int index=1;index<=T;index++)
    {
//        printf("Case #%lld: ",index);
        solve();
//        puts("");
    }
    return 0;
}
/*

5
2 1 2 2 2
1 
1 5 2 1




*/