按照他们的代码稍微修改并加上自我理解的注释的版本
首先, 对于sta->x1->x2->x3->…->end->sta的成环形式,将sta作为终点,即dp[s][sta]表示以状态s中最小点为起点,sta为终点的路径数,如果还存在sta到起点的路径则加上环,我们再枚举j来表示下一个点(此处的j应大于起点,若是等于则是对应成环的情况),也就是枚举i所连的非该集合的点,要注意j不在当前点集,那么状态转移方程式就可以出来了:
f[S|(1<<j)][j]+=f[S][i](j ∉S and i,j相连),相当于固定起点(防止同一个环被记录多次,比如1->2->1,因为1最小所以是起点,而如果不固定起点2->1->2也会被算一次)寻找新终点,把sta连到j处。
注意点集中只有两点的情况会被算上(有条边1到2,则2也能到1,因为不存在重复边所以必然是一条边用了两次,这种与环无关的情况会被记录),所以要减掉边数m 。
同时每个环会被算两次(顺时针和逆时针),最后还要除以二。
#include <bits/stdc++.h>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
#define mem(x, a) memset(x, a, sizeof(x));
#define N 19
int n, m;
bool Map[N + 5][N + 5]; //finds(S):找S二进制中最低非零位的位置
ll f[(1 << N) + 5][N + 5];
int finds(int x)
{
for (int i = 0; i < n; i++)
if (x & (1 << i))
return i;
}
int main()
{ //f[S][i]:点集S,i作为路径的终点,以finds(i)作为起点
ll ans = 0;
cin >> n >> m;
for (int i = 1; i <= m; ++i)
{ //这里以所有点以0开始会好弄一些
int x, y;
cin >> x >> y;
x -= 1;
y -= 1;
Map[x][y] = Map[y][x] = 1;
}
for (int i = 0; i < n; i++)
f[1 << i][i] = 1; //初始化
for (int S = 1; S <= (1 << n) - 1; S++)
{ //枚举点集
for (int i = 0; i < n; i++)
{ //刷表只能由自己转移到别人
if (!f[S][i])
continue; //自己状态必须合法访问过
int p = finds(S); //p为当前集合S位置最小的一点,即起点
for (int j = p; j < n; j++)
{ //固定起点寻找终点
if (!Map[i][j])
continue; //必须与之前的终点连边
if (S & (1 << j))
{ //j访问到当前点集
if (j == p) //j正好为起点则出现了圈
ans += f[S][i]; //累加答案
}
else //状态转移
f[S | (1 << j)][j] += f[S][i];
}
}
} //输出答案
printf("%lld\n", (ans - m) / 2ll);
return 0;
}