[WC2018]州区划分

[WC2018]州区划分

注意审题：

1.有序选择

2.若干个州

3.贡献是州满意度的乘积

枚举最后一个州是哪一个，合法时候贡献sum[s]^p，否则贡献0

存在欧拉回路：每个点都是偶度数，且图连通（dfs验证）

然后愉快子集卷积即可。

PS:
FMT辣鸡，

FWT可以节省一倍的常数！这样才能通过此题

#include<bits/stdc++.h>
#define reg register int
#define il inline
#define fi first
#define se second
#define mk(a,b) make_pair(a,b)
#define numb (ch^'0')
using namespace std;
typedef long long ll;
template<class T>il void rd(T &x){
    char ch;x=;bool fl=false;
    while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
    for(x=numb;isdigit(ch=getchar());x=x*+numb);
    (fl==true)&&(x=-x);
}
template<class T>il void output(T x){if(x/)output(x/);putchar(x%+'');}
template<class T>il void ot(T x){if(x<) putchar('-'),x=-x;output(x);putchar(' ');}
template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');}

namespace Miracle{
const int N=;
const int mod=;
int n,m,p;
int qm(int x,int y){
    int ret=;
    while(y){
        if(y&) ret=(ll)ret*x%mod;
        x=(ll)x*x%mod;
        y>>=;
    }
    return ret;
}
int ad(int x,int y){
    return x+y>=mod?x+y-mod:x+y;
}
int f[N][<<];
int g[N][<<];
int w[N],sum[<<];
int sz[<<];
int to[N];
void FWT(int *f,int n){
    for(reg p=;p<=n;p<<=){
        for(reg l=;l<n;l+=p){
            for(reg k=l;k<l+p/;++k){
                f[k+p/]=ad(f[k],f[k+p/]);
            }
        }
    }
}
void IFWT(int *f,int n){
    for(reg p=;p<=n;p<<=){
        for(reg l=;l<n;l+=p){
            for(reg k=l;k<l+p/;++k){
                f[k+p/]=ad(f[k+p/],mod-f[k]);
            }
        }
    }
}
int val(int v){
    // cout<<" bv "<<v<<endl;
    return p==?:(p==?v:(ll)v*v%mod);
}
int tot;
int vis[N];
int inv[];
void dfs(int x,int s){
    // cout<<" x "<<x<<" s "<<s<<endl;
    vis[x]=;
    ++tot;
    for(reg i=;i<n;++i){
        if((!vis[i])&&(s&(<<i))&&(to[x]&(<<i))){
            dfs(i,s);
        }
    }
}
int main(){
    rd(n);rd(m);rd(p);
    int x,y;
    inv[]=;
    for(reg i=;i<=;++i){
        inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod;
    }
    for(reg i=;i<=m;++i){
        rd(x);rd(y);--x;--y;
        to[x]|=(<<y);to[y]|=(<<x);
    }

    for(reg i=;i<n;++i){
        rd(w[i]);
    }
    for(reg i=;i<(<<n);++i){
        sz[i]=sz[i>>]+(i&);
        for(reg j=;j<n;++j){
            if(i&(<<j)) sum[i]+=w[j];
        }
    }
    // cout<<" pre 1"<<endl;
    for(reg s=;s<(<<n);++s){
        int ji=;
        tot=;
        memset(vis,,sizeof vis);

        dfs((__builtin_ctz(s&(-s))),s);
        // cout<<" s "<<s<<" "<<__builtin_ctz(s&(-s))<<" tot "<<tot<<endl;
        // cout<<" tot "<<tot<<" sz "<<sz[s]<<endl;
        if(tot!=sz[s]) ji=;
        for(reg i=;i<n;++i){
            if(s&(<<i)){
                ji+=(sz[to[i]&s]&);
            }
        }
        // cout<<" ss "<<s<<" ji "<<ji<<endl;
        if(ji>){
            // cout<<" val "<<val(sum[s])<<endl;
            g[sz[s]][s]=val(sum[s]);
        }
    }
    for(reg i=;i<(<<n);++i){
        f[][i]=;
    }
    // return 0;
    // cout<<" built "<<endl;
    for(reg i=;i<=n;++i){
        FWT(g[i],(<<n));
        // prt(g[i],0,(1<<n)-1);

        for(reg j=;j<=i;++j){
            for(reg s=;s<(<<n);++s){
                f[i][s]=ad(f[i][s],(ll)g[j][s]*f[i-j][s]%mod);
            }
        }
        IFWT(f[i],(<<n));
        // prt(f[i],0,(1<<n)-1);
        for(reg s=;s<(<<n);++s){
            // if(sz[s]==i)
            f[i][s]=(ll)val(inv[sum[s]])*f[i][s]%mod;
            // else f[i][s]=0;
        }
        if(i!=n) FWT(f[i],(<<n));
    }
    printf("%d",f[n][(<<n)-]);
    return ;
}

}
signed main(){
    Miracle::main();
    return ;
}

/*
   Author: *Miracle*
   Date: 2019/4/13 19:58:12
*/