【折半 + 高维前缀和】The 2019 Asia Yinchuan First Round Online Programming? - K - Vertex Covers

题目链接 https://nanti.jisuanke.com/t/41295

---

题意

给出一个n个点的图(n<=36)，定义一个点覆盖的价值是所有点权值的积，问所有点覆盖的价值的和是多少。

---

题解

看错题好难受。
折半，把n拆成左右两半。枚举左边点的情况，对于合法情况（包含所有 端点都在左边的边）处理出来，并用高维前缀和做一个统计。
右边也把合法情况都处理出来，对于每种合法情况，如果此时存在一些连向两边的边，而这些边的右端点不在我们枚举的情况内，左端点就必须选择。这样对于每种右边的合法情况，能够确定左边必须选的点集，我们可以通过刚刚求的高维前缀和，来得到包含这些必选点集的值。
复杂度大概是$T{2}^{\frac{n}{2}}m$

---

#include <bits/stdc++.h>
using namespace std;
const int N=1e6+7;
int dp[N];
int t,n,m,q,cs;
struct Edge{
    int v,nxt;
}e[3000];
int p[50],edn;
void add(int u,int v){
    e[++edn]=(Edge){v,p[u]};p[u]=edn;
    e[++edn]=(Edge){u,p[v]};p[v]=edn;
}
int val[50];
int main()
{
    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d",&n,&m,&q);
        for(int i=1;i<=n;i++){
            scanf("%d",&val[i]);
            p[i]=-1;
        }
        edn=-1;
        for(int i=1,u,v;i<=m;i++){
            scanf("%d%d",&u,&v);
            add(u,v);
        }
        int n1=n/2;
        int n2=n-n1;
        for(int s=0;s<(1<<n1);s++){
            bool ok=true;
            dp[s]=1;
            for(int u=1;u<=n1;u++){
                for(int i=p[u];~i;i=e[i].nxt){
                    int v=e[i].v;
                    if(v>n1) continue;
                    if((s&(1<<(u-1)))||(s&(1<<(v-1)))) continue;
                    ok=false;break;
                }
                if(!ok){dp[s]=0;break;}
                if(s&(1<<(u-1))) dp[s]=1ll*dp[s]*val[u]%q;
            }
        }
        for(int i=0;i<n1;i++){
            for(int j=(1<<n1)-1;j>=0;j--){
                if(!(j&(1<<i))) dp[j]=(dp[j]+dp[j|(1<<i)])%q;
            }
        }
        int ans=0;
        for(int s=0;s<(1<<n2);s++){
            bool ok=true;
            int tmp=1;
            int link=0;
            for(int u=n1+1;u<=n1+n2;u++){
                for(int i=p[u];~i;i=e[i].nxt){
                    int v=e[i].v;
                    if(v<=n1){
                        if(!(s&(1<<(u-1-n1)))) link|=1<<(v-1);
                        continue;
                    }
                    if((s&(1<<(u-1-n1)))||(s&(1<<(v-1-n1)))) continue;
                    ok=false;break;
                }
                if(!ok){tmp=0;break;}
                if(s&(1<<(u-1-n1))) tmp=1ll*tmp*val[u]%q;
            }
            ans=(ans+1ll*tmp*dp[link]%q)%q;
        }
        printf("Case #%d: %d\n",++cs,ans);
    }
}