圆方树学习笔记[APIO2018] Duathlon 铁人两项

版权声明：LeoJAM Presents https://blog.csdn.net/fcb_x/article/details/82150370

终于学了圆方树

Achen和Hycc两个坑逼骗我圆切树这个算法

还好发现切点这个东西没意义

圆方树并不只是解决仙人掌的算法，只是它可以解决更多问题，主要涵盖无向联通图问题。

精髓是把BCC（点双）缩成一个菊花图

一下是我魔改的：

主要是为了体现菊花所以我画的很毒瘤

然后这就成了一棵树了

然后DP就好了

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
typedef int INT;
#define int long long
using namespace std;
const int N=1e6+1000;
inline void read(int &x){
    x=0;
    char ch=getchar();
    int f=1;
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    x*=f;
}
struct Front_star{
    int u,v,nxt;
}e[N<<2];
int cnt=0;
int first[N];
void add(int u,int v){
    cnt++;
    e[cnt].u=u;
    e[cnt].v=v;
    e[cnt].nxt=first[u];
    first[u]=cnt;
}
int head[N];
void Insert(int u,int v){
    cnt++;
    e[cnt].u=u;
    e[cnt].v=v;
    e[cnt].nxt=head[u];
    head[u]=cnt;
}
int n,m;
int dfn[N];
int low[N];
int tot=0,Id;
stack<int>S;
int val[N];
int siz[N];
int ans=0;
void Tarjan(int u){
    dfn[u]=low[u]=++tot;
    S.push(u);
    siz[u]=1;
    val[u]=-1;
    for(int i=first[u];i;i=e[i].nxt){
        int v=e[i].v;
        if(!dfn[v]){
            Tarjan(v);
            low[u]=min(low[u],low[v]);
            if(low[v]>=dfn[u]){
                Insert(u,++Id);
                val[Id]=1;
                int x;
                do{
                    x=S.top();
                    Insert(Id,x);
                    S.pop();
                    siz[Id]+=siz[x];
                    val[Id]++;
                }while(x!=v);
                siz[u]+=siz[Id];
            }
        }
        else low[u]=min(low[u],dfn[v]);
    }
}
void dfs(int u,int sum){
//	cout<<u<<'\n';
    if(u<=n){
        ans+=(sum-1)*val[u];
    }
    ans+=(sum-siz[u])*siz[u]*val[u];
    for(int i=head[u];i;i=e[i].nxt){
        int v=e[i].v;
        ans+=(sum-siz[v])*siz[v]*val[u];
        dfs(v,sum);
    }
}
INT main(){
//	freopen("P4630.in","r",stdin);
    read(n);
    read(m);
    Id=n;
    for(int i=1;i<=m;i++){
        int u,v;
        read(u);
        read(v);
        add(u,v);
        add(v,u);
    }
    for(int i=1;i<=n;i++){
        if(!dfn[i]){
            Tarjan(i);
            dfs(i,siz[i]);
        }
    }
    cout<<ans;
}