f[i][j]表示走到j點時，狀态為i的方案數量
M=1<<21 為22位的二進制數，M-1為21位1，表示最終狀态
轉移方程為：
（轉移方程的限制條件見代碼注釋）
答案：881012367360
#include<iostream>
#include<algorithm>
#include<map>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;
const int N=21,M=1<<N;
//22位二進制 M-1 = 21個1
ll f[M][N];
bool g[N][N];
int gcd(int a,int b)
{
return b?gcd(b,a%b):a;
}
ll dp()
{
f[1][0]=1;//最初從0點開始，狀态為00000....1
for(int i=1;i<=M-1;i++)//狀态
for(int j=0;j<=20;j++)//點
if(i>>j&1) //找到狀态i中走過的點j
for(int k=0;k<=20;k++)//對于j點，找能轉移到j的點k
//在未走j之前，狀态為 i-(1<<j) 判斷是否走過k點
if( (i-(1<<j)) >>k&1 && g[k][j])
f[i][j]+=f[(i-(1<<j))][k];
ll res=0;
for(int i=1;i<=20;i++)//倒數第二步在2~21點的情況(對應1~20)
res+=f[M-1][i];
return res;
}
int main()
{
for(int i=1;i<=21;i++)
for(int j=1;j<=21;j++)
if(gcd(i,j)==1)//将1~21映射到0~20
g[i-1][j-1]=true;
cout<<dp();
}