回路计数问题代码详解

#include <bits/stdc++.h>
#define int long long
using namespace std;
// dp[i][j]：表示访问状态为i、当前位于j的方案数(i为一个21位二进制数，表示每个位置是否被访问过)
int dp[1 << 22][22] , g[22][22];
signed main()
{
    int n = 1 << 21;
		// 计算i，j之间有无通路
		for(int i = 1 ; i <= 21 ; i ++) for(int j = 1 ; j <= 21 ; j ++){
        if(__gcd(i , j) == 1) g[i - 1][j - 1] = g[j - 1][i - 1] = 1;
    }
		// 从0号楼出发，设置方案数为1
    dp[1][0] = 1;
    // 遍历21栋楼的所有访问状态，且这个遍历顺序必须从小到大，因为递推关系是小的推出大的
    for(int i = 1 ; i < n ; i ++){
				// 遍历出发点
        for(int j = 0 ; j < 21 ; j ++){
						// 保证不自相矛盾：即出发点是访问过的
            if(!(i >> j & 1)) continue ;
						// 遍历到达点
            for(int k = 0 ; k < 21 ; k ++) {
								// 1.保证出发点可到到达点	2.保证不自相矛盾：即到达点是未访问过的
                if(!g[j][k] || (i >> k & 1)) continue ;
                dp[i + (1 << k)][k] += dp[i][j];
            }
        }
    }
		// 整体来看这个遍历过程：先考虑少数楼层的转移方案数，由小得大，推出整体的转移方案数
    int res = 0;
    // 统计每栋楼都访问一遍之后位于每栋楼的方案数，求和即所求答案
    for(int i = 0 ; i < 21 ; i ++) {
        res += dp[n - 1][i];
    }
    cout << res << '\n';
    return 0;
}