【NOIP模拟】管道

题目：管道

解析：

       状压+记忆化搜索

 

代码：

#include <bits/stdc++.h>
using namespace std;
 
const int RLEN=1<<18|1;
inline char nc() {
	static char ibuf[RLEN],*ib,*ob;
	(ib==ob) && (ob==(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
	return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
	char ch=nc(); int i=0,f=1;
	while(!isdigit(ch)) {if(ch=='-')f=-1; ch=nc();}
	while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=nc();}
	return i*f;
}

const int N=20, mod=998244353;
 
int n,m;
vector <int> g[N];
 
int f[1<<N][N];
 
 
inline int dfs(int mask,int v) {
    if(f[mask][v]) return f[mask][v];
    f[mask][v]=mask;
    for(int i=0;i<g[v].size();++i) {
        int t=g[v][i];
        if(!(mask&(1<<t))) f[mask][v]|=dfs(mask|(1<<t),t);
    }
    return f[mask][v];
}
 
long long dp[1<<N][N];
 
inline int DFS(int mask,int v) {
    if(f[mask][v]==mask) return 1;
    if(dp[mask][v]!=-1) return dp[mask][v];
    dp[mask][v]=0;
    for(int i=0;i<g[v].size();++i) {
        int t=g[v][i];
        if(mask&(1<<t)) continue;
        int x=DFS(mask|(1<<t),t);
        int y=DFS(f[mask|(1<<t)][t],v);
        dp[mask][v]+=(long long)x*y%mod;
    }
    return dp[mask][v]%=mod;
}
 
int main() {
	n=rd(), m=rd();
	for(int i=1;i<=m;i++) {
		int x=rd()-1, y=rd()-1;
		g[x].push_back(y);
		g[y].push_back(x);
	}
    for(int i=0;i<(1<<n);++i) 
        for(int j=0;j<n;++j) 
            if(i&(1<<j)) f[i][j]=dfs(i,j);
    memset(dp,-1,sizeof(dp));
    long long ans=0;
    for(int i=0;i<n;++i) 
        ans=ans+DFS(1<<i,i);
    cout<<ans%mod<<'\n';
}