【状压+容斥】BZOJ4455 [ZJOI2016] 小星星

【题目】
lydsy
一幅$n$个点的树和一棵$n$个点的树，问有多少种点对应方式使得树上存在一条边则图上也存在一条边。
$n\le 17$

【解题思路】
直接求一一对应并不好算，我们可以求任意对应再减去不合法。考虑没有一一对应的关系是，我们可以先暴力枚举以下哪些点可以和树上点对应，其他点则不能对应（相当于枚举一个映射）。

令$f_{i,j}$表示$i$子树对应图上点是$j$的方案数。

那么每次暴力枚举和当前点匹配的点，再找到这个点在图中所有儿子，用子树进行转移，一次复杂度$O(n^3)$。

因此总的复杂度就是$O(n^3\cdot 2^n)$

【参考代码】

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N=20;
int n,m,cnt,a[N];
ll ans,f[N][N];

int read()
{
    int ret=0;char c=getchar();
    while(!isdigit(c)) c=getchar();
    while(isdigit(c)) ret=ret*10+(c^48),c=getchar();
    return ret;
}

vector<int>G[N];
void add1(int u,int v){G[u].push_back(v);G[v].push_back(u);}

int head[N],tot;
struct Tway{int v,nex;}e[N<<1];
void add2(int u,int v)
{
    e[++tot]=(Tway){v,head[u]};head[u]=tot;
    e[++tot]=(Tway){u,head[v]};head[v]=tot;
}

void dfs(int x,int fa)
{
    for(int i=head[x];i;i=e[i].nex) if(e[i].v^fa) dfs(e[i].v,x);
    for(int i=1;i<=cnt;++i)
    {
        f[x][a[i]]=1;
        for(int j=head[x];j;j=e[j].nex)
        {
            int v=e[j].v;
            if(v==fa) continue;
            ll sum=0;
            for(int k=0;k<(int)G[a[i]].size();++k) sum+=f[v][G[a[i]][k]];
            f[x][a[i]]*=sum;
            if(!f[x][a[i]]) break;
        }
    }
}

int main()
{
#ifdef Durant_Lee
    freopen("BZOJ4455.in","r",stdin);
    freopen("BZOJ4455.out","w",stdout);
#endif
    n=read();m=read();
    for(int i=1;i<=m;++i) add1(read(),read());
    for(int i=1;i<n;++i) add2(read(),read());
    for(int s=1;s<(1<<n);++s)
    {
        cnt=0;memset(f,0,sizeof(f));
        for(int i=0;i<n;++i) if(s&(1<<i)) a[++cnt]=i+1;
        dfs(1,0);ll sum=0;
        for(int i=1;i<=cnt;++i) sum+=f[1][a[i]];
        ans+=((n-cnt)&1?-1:1)*sum;
    }
    printf("%lld\n",ans);
    return 0;
}