插头DP。。
求用若干个回路走遍整张图的方案数。。然而格子是六边形的TAT
想了好久...按照列来转移好写一点。。其实和正方形格子差不多？
膜了标程才知道，如果回路数没有限制的话。。可以只记录有没有插头，而不用去记录插头属于哪个联通块。。
（要是去记录插头属于哪个联通块就会爆long long了TAT
代码依然很慢。。。
1 #include<cstdio>
2 #include<iostream>
3 #include<cstring>
4 #include<algorithm>
5 #define ll long long
6 using namespace std;
7 const int maxzt=1023333,modd=6023,inf=1002333333;
8 struct zs1{
9 struct zs{
10 int too,pre;
11 }e[maxzt];int tot,last[modd];
12 int zt[maxzt];
13 ll f[maxzt];
14 inline int get(int v){
15 int i,x=v%modd;
16 for(i=last[x];i&&e[i].too!=v;i=e[i].pre);
17 if(i)return i;
18 e[++tot].too=v,e[tot].pre=last[x],last[x]=tot;
19 f[tot]=0,zt[tot]=v;
20 return tot;
21 }
22 }hm[2];
23 bool can[13][13];
24 int mp[21];
25 int i,j,k,n,m,ans;
26
27 inline void upd(int &a,int b){if(b<a)a=b;}
28
29 inline int encode(){
30 int x=0;
31 for(int i=0;i<n<<1;i++)x=x<<1|mp[i];return x;
32 }
33 inline void decode(int x){
34 for(int i=n*2-1;i>=0;i--)mp[i]=x&1,x>>=1;
35 }
36 inline void clr(bool now){
37 memset(hm[now].last,0,modd<<2);
38 hm[now].tot=0;
39 }
40 inline void dp_blank(int x,int y,bool pre){
41 int up,zx,zs,i,zt;ll f;bool now=pre^1;int posup,poszs,poszx,ys,yx;
42 if(y&1)ys=x;else ys=x-1;yx=ys+1;
43 clr(now);//printf(" %d,%d\n",x,y);
44 for(i=1;i<=hm[pre].tot;i++){
45 zt=hm[pre].zt[i],f=hm[pre].f[i];
46 if(x==1)zt>>=1;
47 decode(zt);
48 // for(int j=0;j<n<<1;j++)printf(" %d",mp[j]);printf(" %lld\n",f);
49 if(y&1)posup=(x<<1)-2,poszs=(x<<1)-1,poszx=x<<1;
50 else posup=max(0,(x<<1)-3),poszs=(x<<1)-2,poszx=(x<<1)-1;
51 up=mp[posup],zs=mp[poszs],zx=mp[poszx];
52 // printf(" %d %d %d %d %d %d\n",up,zs,zx,posup,poszs,poszx);
53 if(up+zs+zx==0){
54 if(can[x+1][y]&&can[ys][y+1])
55 mp[poszs]=0,mp[poszx]=mp[posup]=1,
56 hm[now].f[ hm[now].get(encode()) ]+=f;
57 if(can[x+1][y]&&can[yx][y+1])
58 mp[posup]=0,mp[poszx]=mp[poszs]=1,
59 hm[now].f[ hm[now].get(encode()) ]+=f;
60 if(can[ys][y+1]&&can[yx][y+1])
61 mp[poszx]=0,mp[posup]=mp[poszs]=1,
62 hm[now].f[ hm[now].get(encode()) ]+=f;
63 }
64 if(up+zs+zx==1){
65 if(can[x+1][y])
66 mp[poszs]=mp[posup]=0,mp[poszx]=1,
67 hm[now].f[ hm[now].get(encode()) ]+=f;
68 if(can[ys][y+1])
69 mp[poszs]=mp[poszx]=0,mp[posup]=1,
70 hm[now].f[ hm[now].get(encode()) ]+=f;
71 if(can[yx][y+1])//puts("yx"),
72 mp[posup]=mp[poszx]=0,mp[poszs]=1,
73 hm[now].f[ hm[now].get(encode()) ]+=f;
74 }
75 if(up+zs+zx==2)
76 mp[poszx]=mp[poszs]=mp[posup]=0,
77 hm[now].f[ hm[now].get(encode()) ]+=f;
78 }
79 }
80 inline void dp_bar(int x,int y,bool pre){
81 int i,up,zs,zx,zt;ll f;bool now=pre^1;int posup,poszs,poszx;
82
83 clr(now);//printf(" # %d,%d\n",x,y);
84 for(i=1;i<=hm[pre].tot;i++){
85 zt=hm[pre].zt[i],f=hm[pre].f[i];
86 if(x==1)zt>>=1;
87 decode(zt);
88 // for(int j=0;j<n<<1;j++)printf(" %d",mp[j]);printf(" %lld\n",f);
89 if(y&1)posup=(x<<1)-2,poszs=(x<<1)-1,poszx=x<<1;
90 else posup=max(0,(x<<1)-3),poszs=(x<<1)-2,poszx=(x<<1)-1;
91 up=mp[posup],zs=mp[poszs],zx=mp[poszx];
92 if(up+zs+zx==0)
93 hm[now].f[ hm[now].get(encode()) ]+=f;
94 }
95 }
96 char s[23];
97 int main(){
98 while(scanf("%d",&n)!=EOF){
99 scanf("%d",&m);
100 memset(can,0,sizeof(can));
101 for(i=1;i<=n;i++)for(j=1;j<=8;j++)can[i][j]=1;
102 for(i=1;i<=m;i++)scanf("%s",s),can[s[0]-‘A‘+1][s[1]-‘A‘+1]=0;
103
104 bool now=1,pre=0;
105 clr(pre),hm[pre].f[ hm[pre].get(0) ]=1;
106 for(j=1;j<=8;j++)for(i=1;i<=n;swap(now,pre),i++)
107 if(can[i][j])dp_blank(i,j,pre);else dp_bar(i,j,pre);
108 ll ans=0;
109 for(i=1;i<=hm[pre].tot;i++)ans+=hm[pre].f[i];
110 printf("%lld\n",ans);
111 }
112 return 0;
113 }
发现取模的数字会对效率有很大影响。。。。平时似乎不会啊= =
模数取10w出头就好了。。
时间： 2024-08-10 23:14:20