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
*/