BZOJ4316: 小C的独立集【Tarjan+DP+仙人掌】

4316: 小C的独立集

如果这是一棵树，那么很好做，设$F[i][0/1]$就可以了。

我们考虑每一个环，环的最末端会对最前端有影响。

最末端是0，无所谓，最末端为1，那么最顶端只能是0。

那我们先处理环外的点，然后考虑一个环，强制最末端为0/1。

为1就强制为$-\infty$ ，分别DP就可以了。

#include<cstdio>
#include<algorithm>
using namespace std;
const int MAXN=5e4+5,MAXE=6e4+5;
struct Edge{
	int tot,lnk[MAXN],son[MAXE<<1],nxt[MAXE<<1];
	void Add(int x,int y){son[++tot]=y;nxt[tot]=lnk[x];lnk[x]=tot;}
}E;
int n,m,Ans,cnt,F[MAXN][2],Dfn[MAXN],Low[MAXN],Fa[MAXN];
#include<cctype>
int read(){
	int ret=0;char ch=getchar();bool f=1;
	for(;!isdigit(ch);ch=getchar()) f^=!(ch^'-');
	for(; isdigit(ch);ch=getchar()) ret=ret*10+ch-48;
	return f?ret:-ret;
}
void DP(int x,int y){
	int f[2],t[2];f[0]=0,f[1]=0;
	for(int i=y;i!=x;i=Fa[i]){
		t[0]=f[0]+F[i][0],t[1]=f[1]+F[i][1];
		f[0]=max(t[0],t[1]);f[1]=t[0];
	}
	F[x][0]+=f[0];
	f[0]=0,f[1]=-1e9;
	for(int i=y;i!=x;i=Fa[i]){
		t[0]=f[0]+F[i][0],t[1]=f[1]+F[i][1];
		f[0]=max(t[0],t[1]);f[1]=t[0];
	}
	F[x][1]+=f[1];
}
void Tarjan(int x,int fa){
	Dfn[x]=Low[x]=++cnt;Fa[x]=fa;F[x][1]=1,F[x][0]=0;
	for(int j=E.lnk[x];j;j=E.nxt[j])
	if(E.son[j]!=fa){
		if(!Dfn[E.son[j]]) Tarjan(E.son[j],x),Low[x]=min(Low[x],Low[E.son[j]]);
		else Low[x]=min(Low[x],Dfn[E.son[j]]);
		if(Dfn[x]<Low[E.son[j]]) F[x][1]=F[E.son[j]][0]+1,F[x][0]=max(F[E.son[j]][1],F[E.son[j]][0]);
	}
	for(int j=E.lnk[x];j;j=E.nxt[j])
	if(Fa[E.son[j]]!=x&&Dfn[x]<Dfn[E.son[j]]) DP(x,E.son[j]);
}
int main(){
//	freopen("4316.in","r",stdin);
//	freopen("4316.out","w",stdout);
	n=read(),m=read();
	for(int i=1;i<=m;i++){int x=read(),y=read();E.Add(x,y),E.Add(y,x);}
	for(int i=1;i<=n;i++) if(!Dfn[i]) Tarjan(i,0),Ans+=max(F[i][0],F[i][1]);
	printf("%d\n",Ans);
	return 0;
}