bzoj3162 独钓寒江雪（树形dp+Hash）

就是求无根树的本质不同的独立集个数。如果不考虑本质不同，则有树形dp。
考虑本质不同，就要考虑无根树的同构，可以用hash来搞。具体题解见：portal

我的hash真是冲突到死【再见】，hash值最好从大往小排。。。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
#define mod 1000000007
#define ull unsigned long long
#define k1 11113
#define k2 11117
inline int read(){int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();return x*f;
}
int n,h[N],num=1,G0=0,G1=0,mn=inf,sz[N],root,a[N],f[2][N],inv[N],prod[N];
//f[0][x],以x为根的子树,不选x的本质不同独立集个数，f[1][x]--选x的个数
ull hs[N];
struct edge{int to,next;
}data[N<<1];
inline void add(int x,int y){data[++num].to=y;data[num].next=h[x];h[x]=num;
}
inline void dfs(int x,int Fa){sz[x]=1;int mx=0;for(int i=h[x];i;i=data[i].next){int y=data[i].to;if(y==Fa) continue;dfs(y,x);sz[x]+=sz[y];mx=max(mx,sz[y]);}mx=max(mx,n-sz[x]);if(mx<mn) mn=mx,G0=x,G1=0;else if(mx==mn) G1=x;
}
inline bool cmp(int a,int b){return hs[a]==hs[b]?f[0][a]<f[0][b]:hs[a]>hs[b];}
inline int C(int n,int m){int res=prod[m];for(int i=0;i<m;++i) res=(ll)res*(n-i)%mod;return res;
}
inline void dfs2(int x,int Fa){for(int i=h[x];i;i=data[i].next){int y=data[i].to;if(y==Fa) continue;dfs2(y,x);}int tot=0;for(int i=h[x];i;i=data[i].next){int y=data[i].to;if(y==Fa) continue;a[++tot]=y;}sort(a+1,a+tot+1,cmp);f[0][x]=f[1][x]=1;int j;for(int i=1;i<=tot;++i){for(j=i+1;j<=tot&&hs[a[j]]==hs[a[i]]&&f[0][a[j]]==f[0][a[i]];++j);--j;f[0][x]=(ll)f[0][x]*C(j-i+1+f[0][a[i]]+f[1][a[i]]-1,j-i+1)%mod;f[1][x]=(ll)f[1][x]*C(j-i+1+f[0][a[i]]-1,j-i+1)%mod;i=j;}hs[x]=2333;for(int i=1;i<=tot;++i) hs[x]=hs[x]*k1+hs[a[i]]*k2;
}
int main(){
//  freopen("a.in","r",stdin);n=read();inv[1]=prod[1]=1;for(int i=2;i<=n;++i) inv[i]=(ll)inv[mod%i]*(mod-mod/i)%mod;for(int i=2;i<=n;++i) prod[i]=(ll)prod[i-1]*inv[i]%mod;for(int i=1;i<n;++i){int x=read(),y=read();add(x,y);add(y,x);}dfs(1,0);if(!G1) root=G0;else{root=n+1;for(int i=h[G0];i;i=data[i].next)if(data[i].to==G1){data[i].to=data[i^1].to=root;break;}add(root,G0);add(root,G1);}dfs2(root,0);if(!G1){printf("%d\n",(f[0][root]+f[1][root])%mod);return 0;}if(hs[G0]==hs[G1]&&f[0][G0]==f[0][G1]){int x=f[0][G0];int ans=(ll)f[0][G0]*f[1][G0]%mod;ans=(ans+(ll)x*(x+1)/2)%mod;printf("%d\n",ans);return 0;}printf("%d\n",((ll)f[0][G0]*f[0][G1]+(ll)f[0][G0]*f[1][G1]+(ll)f[1][G0]*f[0][G1])%mod);return 0;
}

bzoj3162 独钓寒江雪（树形dp+Hash）相关推荐

洛谷4895 BZOJ3162 独钓寒江雪树形dp 树哈希
题目链接题意: 给定一棵无根树,求其中本质不同的独立集的个数.独立集就是一个集合中的点之间都没有边直接相连.n<=5e5n<=5e5n<=5e5,对1e9+71e9+71e9+7取 ...
POJ 3342 树形DP+Hash
这是很久很久以前做的一道题,可惜当时WA了一页以后放弃了. 今天我又重新捡了起来.(哈哈1A了) 题意: 没有上司的舞会+判重思路: hash一下+树形DP 题目中给的人名hash到数字,再进行运算 ...
bzoj 3162: 独钓寒江雪树哈希+树形dp
题意给出一棵无标号无根树,问本质不同的最大独立集数量.答案模1e9+7. n<=500000 分析对于一般的情况,我们可以先找出树的重心作为根,然后进行树形dp.这样做有什么好处呢?通过根的 ...
codeforce708C:树形dp+二次扫描
题目链接题目大意: 给你一棵树,问你对于每一个点是否可以在树上删掉一条边,再增加一条边,使它成为树的重心. tips:一个点成为重心的条件为它的每个子树大小不超过n/2. 解题思路:首先树上问题一般 ...
[BZOJ3197][Sdoi2013]assassin（树形DP+树同构+二分图最优匹配）
关于树同构,有一个神奇的性质: 一棵树的重心只有 111 个或 2" role="presentation" style="position: relative ...
BNUOJ 52305 Around the World 树形dp
题目链接: https://www.bnuoj.com/v3/problem_show.php?pid=52305 Around the World Time Limit: 20000msMemory ...
[树形dp] Jzoj P5233 概率博弈
Description 小A和小B在玩游戏.这个游戏是这样的: 有一棵n个点的以1为根的有根树,叶子有权值.假设有m个叶子,那么树上每个叶子的权值序列就是一个1->m 的排列. 一开始在1号点有 ...
fwt优化+树形DP HDU 5909
1 //fwt优化+树形DP HDU 5909 2 //见官方题解 3 // BestCoder Round #88 http://bestcoder.hdu.edu.cn/ 4 5 #include ...
BZOJ 1040 ZJOI2008 骑士树形DP
题目大意:给定一个基环树林,每一个点上有权值,要求选择一个权值和最大的点集,要求点集中的随意两个点之间不能直接相连最大点独立集--考虑到n<=100W,网络流铁定跑不了,于是我们考虑树形DP ...

bzoj3162 独钓寒江雪（树形dp+Hash）

bzoj3162 独钓寒江雪（树形dp+Hash）相关推荐

最新文章

热门文章