Formula 1

题意

在\(n*m\)的矩阵中，有些格子有树，没有树的格子不能到达，找一条回路，吃完所有的树，求有多少种方法。

解法

因为只要一条回路，所以我们必须维护插头的连通性。
具体的可以参照 这位大佬的博客

代码

注意开long long。

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <cctype>
#define del(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef long long ll;
template <typename T>
inline void read(T &x) {x=0;char c=getchar();T k=1;while(!isdigit(c)) {if(c=='-') k=-1;c=getchar();}while(isdigit(c)) {x=x*10+c-'0';c=getchar();}x*=k;
}const int maxn=15;
const int maxhash=100000;
char G[maxn][maxn];
int _hash[maxhash];
ll sta[2][600000],sum[2][600000];
int cur,n,m,en,em;
int tot[2];
int jz[maxn];void _init() {read(n),read(m);for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) {cin>>G[i][j];if(G[i][j]=='.') en=i,em=j;}for(int i=1;i<=m;i++) jz[i]=i<<1;
}void hash_insert(ll s,ll data) {int pos=s%maxhash;while(_hash[pos]) {if(sta[cur][_hash[pos]]==s) {sum[cur][_hash[pos]]+=data;return;}if(++pos==maxhash) pos=0;}++tot[cur];_hash[pos]=tot[cur];sta[cur][tot[cur]]=s;sum[cur][tot[cur]]=data;
}
ll ans;
void work() {tot[0]=1;sum[0][1]=1;for(int i=1;i<=n;i++) {for(int k=1;k<=tot[cur];k++)sta[cur][k]=sta[cur][k]<<2;for(int j=1;j<=m;j++) {cur^=1;tot[cur]=0;del(_hash,0);del(sta[cur],0);del(sum[cur],0);for(int k=1;k<=tot[1-cur];k++) {ll s=sta[1-cur][k],data=sum[1-cur][k];int x=(s>>jz[j-1])%4;int y=(s>>jz[j])%4;ll temp;if(G[i][j]!='.') {if(x==0&&y==0) hash_insert(s,data);}else {if(x==0&&y==0) {if(G[i][j+1]=='.'&&G[i+1][j]=='.') {temp=s+1*(1<<jz[j-1])+2*(1<<jz[j]);                     hash_insert(temp,data);}continue;}if(x==0&&y>0) {if(G[i][j+1]=='.')hash_insert(s,data);if(G[i+1][j]=='.') {temp=s-y*(1<<jz[j])+y*(1<<jz[j-1]);                         hash_insert(temp,data);}continue;}if(x>0&&y==0) {if(G[i+1][j]=='.')hash_insert(s,data);if(G[i][j+1]=='.') {temp=s-x*(1<<jz[j-1])+x*(1<<jz[j]);                     hash_insert(temp,data);}continue;}if(x==1&&y==1) {int f=1;for(int v=j+1;v<=m;v++) {int fff=(s>>jz[v])%4;if(fff==1) f++;if(fff==2) f--;if(!f) {temp=s-2*(1<<jz[v])+1*(1<<jz[v]);break;}}temp=temp-1*(1<<jz[j-1])-1*(1<<jz[j]);                  hash_insert(temp,data);continue;}if(x==2&&y==2) {int f=1;for(int v=j-2;v>=1;v--) {int fff=(s>>jz[v])%4;if(fff==1) f--;if(fff==2) f++;if(!f) {temp=s-1*(1<<jz[v])+2*(1<<jz[v]);break;}}temp=temp-2*(1<<jz[j-1])-2*(1<<jz[j]);      hash_insert(temp,data);continue;}if(x==2&&y==1) {temp=s-2*(1<<jz[j-1])-1*(1<<jz[j]);                 hash_insert(temp,data);continue;}if(x==1&&y==2) {if(i==en&&j==em) {ans+=data;}}}}}}
}int main() {_init();work();printf("%lld\n",ans);return 0;
}