bzoj4316 小C的独立集
Posted suxxsfe
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj4316 小C的独立集相关的知识,希望对你有一定的参考价值。
https://darkbzoj.tk/problem/4316
求一个仙人掌的最大独立集
先把他建出圆方树来,每个环选一个点当做“这个环的根”,作为对应方点的父亲,其他换上的点作为这个方点的儿子
考虑用 (f(u,1/2)) 来表示 (u) 的子树中,(u) 这个点选/不选的最大独立集大小
如何转移?圆点很好转移,(f(u,0)=sum max(f(v,1),f(v,0)),f(u,1)=1+sum f(v,0))
至于方点如何转移,要弄清它在此的实际意义,考虑一下这个方点的父亲(环的根)的转移情况
对于 (f(u,0)),它会为 (f(fa,1)) 产生贡献,那么就要选上环的根,则根旁边的两个点就都不能选。就是环的根以下的部分,不选和环的根相邻的两点,能选出的最大独立集大小
对于 (f(u,1)),它只有可能为 (f(fa,0)) 产生贡献,所以就不用管和不和环的根相连了。就是环的根一下的部分,能选出的最大独立集大小
至于如何求,在 DP()
函数里用正反两遍 dp 解决
其实这种仙人掌转圆方树的题似乎都是这个套路,就是转完圆方树以后,圆点的转移一般比较简单,而方点的可能还需要个 dp 啥的
比如这个题也是这样,只不过有些难我没写出来,dp 不好啥都玩蛋啊:https://www.luogu.com.cn/problem/P4244
#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<map>
#include<iomanip>
#include<cstring>
#define reg register
#define EN puts("")
inline int read(){
register int x=0;register int y=1;
register char c=std::getchar();
while(c<‘0‘||c>‘9‘){if(c==‘-‘) y=0;c=std::getchar();}
while(c>=‘0‘&&c<=‘9‘){x=x*10+(c^48);c=std::getchar();}
return y?x:-x;
}
#define N 100006
#define M 300006
struct graph{
int fir[N],nex[M],to[M],tot;
inline void add(int u,int v){
to[++tot]=v;
nex[tot]=fir[u];fir[u]=tot;
}
}G,T;
int n,m,bcccnt;
int dfn[N],low[N],dfscnt;
int stack[N],top;
void tarjan(reg int u,int fa){
dfn[u]=low[u]=++dfscnt;stack[top++]=u;
for(reg int v,i=G.fir[u];i;i=G.nex[i]){
v=G.to[i];
if(v==fa) continue;
if(!dfn[v]){
tarjan(v,u);
low[u]=std::min(low[u],low[v]);
if(low[v]>dfn[u]) T.add(u,v);
}
else if(low[u]>dfn[v]){
low[u]=dfn[v];bcccnt++;
T.add(v,bcccnt);
for(reg int j=top-1;stack[j]^v;j--) T.add(bcccnt,stack[j]);
}
}
top--;
}
int tmp[N];
int f[N][2],dp[N][2];
inline void DP(reg int u){
tmp[0]=0;
for(reg int i=T.fir[u];i;i=T.nex[i]) tmp[++tmp[0]]=T.to[i];
if(tmp[0]==2){
int son1=tmp[1],son2=tmp[2];
f[u][0]=f[son1][0]+f[son2][0];//son1 son2 都不能选
f[u][1]=std::max(f[son1][0]+f[son2][0],std::max(f[son1][0]+f[son2][1],f[son1][1]+f[son2][0]));
return;
}
dp[1][0]=f[tmp[1]][0];dp[1][1]=0;//第一个不能选,所以是 0
for(reg int v,i=2;i<=tmp[0];i++){
v=tmp[i];
dp[i][0]=std::max(dp[i-1][0],dp[i-1][1])+f[v][0];
dp[i][1]=dp[i-1][0]+f[v][1];
}
f[u][0]=dp[tmp[0]][0];f[u][1]=dp[tmp[0]][0];
dp[tmp[0]+1][0]=dp[tmp[0]+1][1]=0;
for(reg int v,i=tmp[0];i;i--){
v=tmp[i];
dp[i][0]=std::max(dp[i+1][0],dp[i+1][1])+f[v][0];
dp[i][1]=dp[i+1][0]+f[v][1];
}
f[u][1]=std::max(f[u][1],std::max(dp[1][0],dp[1][1]));
}
void dfs(int u){
f[u][1]=1;
for(reg int v,i=T.fir[u];i;i=T.nex[i]){
v=T.to[i];
dfs(v);
if(u<=n){
f[u][0]+=std::max(f[v][1],f[v][0]);f[u][1]+=f[v][0];
}
}
if(u>n) DP(u);
}
inline void debug_tarjan(){
puts("finished tarjan");
printf("bcccnt=%d
",bcccnt);
for(reg int i=1;i<=bcccnt;i++){
printf("i=%d
",i);
for(reg int j=T.fir[i];j;j=T.nex[j]) printf("%d ",T.to[j]);
puts("");
}
}
int main(){
n=read();m=read();
for(reg int u,v,i=1;i<=m;i++){
u=read();v=read();
G.add(u,v);G.add(v,u);
}
bcccnt=n;
tarjan(1,1);
// debug_tarjan();
dfs(1);
printf("%d",std::max(f[1][0],f[1][1]));
return 0;
}
以上是关于bzoj4316 小C的独立集的主要内容,如果未能解决你的问题,请参考以下文章