终于学了圆方树
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;
}