题解-AGC016F

早期 AGC 的难度确实虚高（但 F 也不是我这种菜鸡能够挑战）。

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int n,m,pw[225],g[20][20],cnt[1<<16][20],f[1<<16];
int lowbit(int s){return s&-s;}
int main(){
	scanf("%d%d",&n,&m);
	pw[0]=1;
	for(int i=1,u,v;i<=m;i++){
		pw[i]=(pw[i-1]+pw[i-1])%mod;
		scanf("%d%d",&u,&v);
		u--;v--;
		g[u][v]=1;
	}
	for(int s=1;s<(1<<n);s++){
		for(int i=0;i<n;i++){
			cnt[s][i]=cnt[s-lowbit(s)][i]+g[i][__builtin_ctz(lowbit(s))];
		}
	}
	for (int s=0;s<(1<<n);s++)if((s&3)==3){
		f[s]=1;
		for (int t=s&(s-1);t;t=(t-1)&s)if((t&3)==3||(t&3)==0){
			int now=1;
			if ((t&3)==3) {
				for(int i=0;i<n;i++)if((s>>i)&1){
					if((t>>i)&1)now=1ll*now*(pw[cnt[s^t][i]]-1)%mod;
					else now=1ll*now*pw[cnt[t][i]]%mod;
				}
				now=1ll*now*f[t]%mod;
			} else {
				for(int i=0;i<n;i++)if((s>>i)&1){
					if((t>>i)&1)now=1ll*now*(pw[cnt[s^t][i]]-1)%mod*pw[cnt[t][i]]%mod;
					else now=1ll*now*pw[cnt[t][i]]%mod;
				}
			}
			f[s]=(f[s]+now)%mod;
		}
	}
	printf("%d\n",(pw[m]-f[(1<<n)-1]+mod)%mod);
	return 0;
}

标签：cnt,pw,int,题解,AGC016F,1ll,now,mod
来源： https://www.cnblogs.com/zcr-blog/p/16343796.html