BZOJ4316 小C的独立集

仙人掌DP，套路还是1023的套路，dp式也还是比较水，f[i][0]表示i的子仙人掌的最大独立集，i不选，f[i][1]表示i的子仙人掌的最大独立集，i选，做环的时候强制令第一个点选或不选做两次

不过交到1040WA了是什么情况

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<cstring>
#include<ctime>
#include<vector>
#include<stack>
#include<queue>
#include<set>
#include<bitset>
#include<map>
using namespace std;
#define MAXN 50010
#define MAXM 120000
#define INF 1000000000
#define MOD 1000000007
#define ll long long
#define eps 1e-8
struct vec{
int to;
int fro;
};
vec mp[MAXM];
int tai[MAXN],cnt;
int dfn[MAXN],low[MAXN],tim;
int fa[MAXN];
int dep[MAXN];
int f[MAXN][2];
int a[MAXN][4];
int n,m;
void dp(int rt,int ed){
int i,k;
a[1][0]=f[ed][0];
a[1][1]=-INF;
a[1][2]=-INF;
a[1][3]=f[ed][1];
for(k=fa[ed],i=2;k!=fa[rt];k=fa[k],i++){
a[i][0]=f[k][0]+max(a[i-1][0],a[i-1][1]);
a[i][1]=f[k][1]+a[i-1][0];
a[i][2]=f[k][0]+max(a[i-1][2],a[i-1][3]);
a[i][3]=f[k][1]+a[i-1][2];
}
f[rt][0]=max(f[rt][0],max(a[i-1][0],a[i-1][2]));
f[rt][1]=max(f[rt][1],a[i-1][1]);
}
void dfs(int x){
int i,y;
dfn[x]=low[x]=++tim;
dep[x]=dep[fa[x]]+1;
f[x][1]=1;
for(i=tai[x];i;i=mp[i].fro){
y=mp[i].to;
if(y==fa[x]){
continue ;
}
if(!dfn[y]){
fa[y]=x;
dfs(y);
}
low[x]=min(low[x],low[y]);
if(low[y]>dfn[x]){
f[x][1]+=f[y][0];
f[x][0]+=max(f[y][0],f[y][1]);
}
}
for(i=tai[x];i;i=mp[i].fro){
y=mp[i].to;
if(fa[y]!=x&&dfn[y]>dfn[x]){
dp(x,y);
}
}
}
inline void be(int x,int y){
mp[++cnt].to=y;
mp[cnt].fro=tai[x];
tai[x]=cnt;
}
inline void bde(int x,int y){
be(x,y);
be(y,x);
}
int main(){
int i,j,x,y;
scanf("%d%d",&n,&m);
for(i=1;i<=m;i++){
scanf("%d%d",&x,&y);
bde(x,y);
}
dfs(1);
printf("%d\n",max(f[1][0],f[1][1]));
return 0;
}

/*

*/