bzoj4316: 小C的獨立集
阿新 • • 發佈:2018-11-27
首先這是個仙人掌,設\(f[i][0/1]\)表示當前節點\(i\),選或不選的最大獨立集
如果某條邊是樹邊,那麼直接樹形dp的轉移即可
考慮如果它的某棵子樹恰好是一個環該怎麼辦
我們列舉這個環的頂端選或者不選,然後從這個環的底部開始遍歷一遍即可
//minamoto #include<bits/stdc++.h> #define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i) #define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i) #define go(u) for(register int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v) using namespace std; #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++) char buf[1<<21],*p1=buf,*p2=buf; template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;} int read(){ int res,f=1;char ch; while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1); for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0'); return res*f; } const int N=2e5+5; struct eg{int v,nx;}e[N];int head[N],tot; inline void add(int u,int v){e[++tot]={v,head[u]},head[u]=tot;} int n,m,fa[N],f[N][2],dfn[N],low[N],tim,u,v; void dp(int u,int v){ int t0,t1,f0=0,f1=0; for(register int i=v;i!=u;i=fa[i]) t0=f0+f[i][0],t1=f1+f[i][1],f0=max(t0,t1),f1=t0; f[u][0]+=f0,f0=0,f1=-1e9; for(register int i=v;i!=u;i=fa[i]) t0=f0+f[i][0],t1=f1+f[i][1],f0=max(t0,t1),f1=t0; f[u][1]+=f1; } void dfs(int u,int fat){ fa[u]=fat,dfn[u]=low[u]=++tim,f[u][1]=1,f[u][0]=0; go(u)if(v!=fat){ if(!dfn[v])dfs(v,u),cmin(low[u],low[v]); else cmin(low[u],dfn[v]); if(low[v]>dfn[u])f[u][1]+=f[v][0],f[u][0]+=max(f[v][0],f[v][1]); } go(u)if(fa[v]!=u&&dfn[u]<dfn[v])dp(u,v); } int main(){ // freopen("testdata.in","r",stdin); n=read(),m=read(); fp(i,1,m)u=read(),v=read(),add(u,v),add(v,u); dfs(1,0);printf("%d\n",max(f[1][0],f[1][1]));return 0; }