BZOJ4316: 小C的独立集
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ4316: 小C的独立集相关的知识,希望对你有一定的参考价值。
圆方树预备知识。
首先搞出DFS树,然后设状态$f[node][OK][OKanc]$为别表示在当前点选/不限和环顶节点选/不选。
剩下的就是TreeDP常规的东西。
P.S.课件上说第一个状态设为父亲的,亲测没有设成自己来得方便
//BZOJ 4316 //by Cydiater //2017.2.12 #include <iostream> #include <queue> #include <map> #include <ctime> #include <cstdlib> #include <cstdio> #include <cstring> #include <string> #include <algorithm> #include <iomanip> #include <cmath> #include <bitset> #include <set> #include <vector> using namespace std; #define ll long long #define up(i,j,n) for(int i=j;i<=n;i++) #define down(i,j,n) for(int i=j;i>=n;i--) #define cmax(a,b) a=max(a,b) #define cmin(a,b) a=min(a,b) #define Auto(i,node) for(int i=LINK[node];i;i=e[i].next) const int MAXN=2e5+5; const int oo=0x3f3f3f3f; inline int read(){ char ch=getchar();int x=0,f=1; while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)f=-1;ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} return x*f; } int LINK[MAXN],N,M,len=0,fa[MAXN],dep[MAXN],anc[MAXN],f[MAXN][2][2],fat[MAXN]; bool istop[MAXN]; struct edge{ int y,next; }e[MAXN]; namespace solution{ inline void insert(int x,int y){e[++len].next=LINK[x];LINK[x]=len;e[len].y=y;} inline void Insert(int x,int y){insert(x,y);insert(y,x);} void DFS(int node,int father,int deep){ dep[node]=deep+1;fa[node]=father; Auto(i,node)if(e[i].y!=father&&!dep[e[i].y]) DFS(e[i].y,node,deep+1); } void Prepare(){ N=read();M=read(); up(i,1,M){ int x=read(),y=read(); Insert(x,y); } DFS(1,0,0); up(node,1,N)Auto(i,node){ if(dep[e[i].y]<dep[node]&&e[i].y!=fa[node]){ istop[e[i].y]=1; anc[node]=e[i].y; int tmp=node; while(tmp!=e[i].y){ fat[tmp]=e[i].y; tmp=fa[tmp]; } } } } int TreeDP(int node,int OK,int OKanc){ if(f[node][OK][OKanc]!=-1)return f[node][OK][OKanc]; f[node][OK][OKanc]=0; if(OK){ if(anc[node]&&OKanc)return f[node][OK][OKanc]; Auto(i,node)if(dep[e[i].y]==dep[node]+1){ if(fat[e[i].y]==node&&istop[node]) f[node][OK][OKanc]+=TreeDP(e[i].y,0,1); else f[node][OK][OKanc]+=TreeDP(e[i].y,0,OKanc); } f[node][OK][OKanc]++; }else{ Auto(i,node)if(dep[e[i].y]==dep[node]+1){ if(fat[e[i].y]==node&&istop[node]){ f[node][OK][OKanc]+=max(TreeDP(e[i].y,0,0),TreeDP(e[i].y,1,0)); } else f[node][OK][OKanc]+=max(TreeDP(e[i].y,0,OKanc),TreeDP(e[i].y,1,OKanc)); } } return f[node][OK][OKanc]; } } int main(){ //freopen("input.in","r",stdin); using namespace solution; Prepare(); memset(f,-1,sizeof(f)); cout<<max(TreeDP(1,1,1),TreeDP(1,0,0))<<endl; return 0; }
以上是关于BZOJ4316: 小C的独立集的主要内容,如果未能解决你的问题,请参考以下文章