您的位置:首页 > 其它

BZOJ 4455: [Zjoi2016]小星星

2016-09-12 21:06 567 查看
liaoliaoSB的题就先做到这了~毕竟后缀数组不好。。借口

这题可是高深的容斥原理啊,你要我证我可证不清楚。

做法: 先不限制一对一,我们可以通过枚举哪几个不选,求出全部集合的并集(至少0个禁止选的-至少1个禁止选的+至少2个禁止选的..),集合i是指不选i的情况

最后补集就是答案了。。。

代码:(短+快~)

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
struct node{int y,next;}a[40]; int first[20],len;
void ins(int x,int y)
{
a[++len].y=y,a[len].next=first[x],first[x]=len;
}
bool v[20][20];
int al,p[20];
LL f[20][20];
void dfs(int x,int fa)
{
int i,j,k,y; LL s;
for(k=first[x];k;k=a[k].next)
{
y=a[k].y;
if(y==fa)continue;
dfs(y,x);
}
for(i=1;i<=al;i++)
{
f[x][i]=1;
for(k=first[x];k;k=a[k].next)
{
y=a[k].y; if(y==fa)continue;
s=0;
for(j=1;j<=al;j++)if(v[p[i]][p[j]])s+=f[y][j];
f[x][i]*=s;
}
}
}
int main()
{
int n,m,x,y,i,j;
scanf("%d%d",&n,&m);
for(i=1;i<=m;i++)
{
scanf("%d%d",&x,&y);
v[x][y]=v[y][x]=1;
}
for(i=1;i<n;i++)scanf("%d%d",&x,&y),ins(x,y),ins(y,x);
LL ans=0,s;
for(i=0;i<(1<<n)-1;i++)
{
al=0;
for(j=1;j<=n;j++)if( !(i&(1<<j-1)) )p[++al]=j;
dfs(1,0);
s=0;
for(j=1;j<=al;j++)s+=f[1][j];
if( (n-al)&1 )ans-=s; else ans+=s;
}
printf("%lld\n",ans);
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: