Ural 1519 Formula 1 插头DP(单回路)

  题目链接:http://acm.timus.ru/problem.aspx?space=1&num=1519

   一条哈密顿回路路,反正是写到蛋疼了,不过终于解决了,插头DP这玩意,太容易出错了,要注意block的处理。

  1 //STATUS:C++_AC_343MS_961KB
  2 #include<stdio.h>
  3 #include<stdlib.h>
  4 #include<string.h>
  5 #include<math.h>
  6 #include<iostream>
  7 #include<string>
  8 #include<algorithm>
  9 #include<vector>
 10 #include<queue>
 11 #include<stack>
 12 #include<map>
 13 using namespace std;
 14 #define LL long long
 15 #define pii pair<int,int>
 16 #define Max(a,b) ((a)>(b)?(a):(b))
 17 #define Min(a,b) ((a)<(b)?(a):(b))
 18 #define mem(a,b) memset(a,b,sizeof(a))
 19 #define lson l,mid,rt<<1
 20 #define rson mid+1,r,rt<<1|1
 21 const int N=15,INF=0x3f3f3f3f,MOD=4001,STA=1000010;
 22 const double DNF=100000000000;
 23 
 24 int g[N][N],code[N],ma[N];
 25 int n,m,ex,ey;
 26 
 27 struct Hash{
 28     int first[MOD],next[STA],size;
 29     LL f[STA],sta[STA];
 30     void init(){
 31         size=0;
 32         mem(first,-1);
 33     }
 34     void add(LL st,LL ans){
 35         int i,u=st%MOD;
 36         for(i=first[u];i!=-1;i=next[i]){
 37             if(sta[i]==st){
 38                 f[i]+=ans;
 39                 return;
 40             }
 41         }
 42         sta[size]=st;
 43         f[size]=ans;
 44         next[size]=first[u];
 45         first[u]=size++;
 46     }
 47 }hs[2];
 48 
 49 void shift(int p)    //换行移位
 50 {
 51     int k;
 52     LL sta;
 53     for(k=0;k<hs[!p].size;k++){
 54         sta=hs[!p].sta[k]<<3;
 55         hs[p].add(sta,hs[!p].f[k]);
 56     }
 57 }
 58 
 59 LL getsta()    //最小表示法
 60 {
 61     LL i,cnt=1,sta=0;
 62     mem(ma,-1);
 63     ma[0]=0;
 64     for(i=0;i<=m;i++){
 65         if(ma[code[i]]==-1)ma[code[i]]=cnt++;
 66         code[i]=ma[code[i]];
 67         sta|=(LL)code[i]<<(3*i);
 68     }
 69     return sta;
 70 }
 71 
 72 void getcode(LL sta)
 73 {
 74     int i;
 75     for(i=0;i<=m;i++){
 76         code[i]=sta&7;
 77         sta>>=3;
 78     }
 79 }
 80 
 81 void unblock(int i,int j,int p)
 82 {
 83     int k,t;
 84     LL cnt,x,y;
 85     for(k=0;k<hs[!p].size;k++){
 86         getcode(hs[!p].sta[k]);
 87         x=code[j],y=code[j+1];
 88         cnt=hs[!p].f[k];
 89         if(x && y){          //合并连通分量
 90             code[j]=code[j+1]=0;
 91             if(x!=y){
 92                 for(t=0;t<=m;t++)
 93                     if(code[t]==y)code[t]=x;
 94                 hs[p].add(getsta(),cnt);
 95             }
 96             else if(i==ex && j==ey){     //最后一个点特殊处理
 97                 hs[p].add(getsta(),cnt);
 98             }
 99         }
100 
101         else if(x&&!y || !x&&y){    //延续连通分量
102             t=x?x:y;
103             if(g[i+1][j]){
104                 code[j]=t;code[j+1]=0;
105                 hs[p].add(getsta(),cnt);
106             }
107             if(g[i][j+1]){
108                 code[j]=0;code[j+1]=t;
109                 hs[p].add(getsta(),cnt);
110             }
111         }
112         else if(g[i+1][j] && g[i][j+1]){   //创建新连通分量
113             code[j]=code[j+1]=8;
114             hs[p].add(getsta(),cnt);
115         }
116     }
117 }
118 
119 void block(LL j,int p)
120 {
121     int k;
122     for(k=0;k<hs[!p].size;k++){
123         getcode(hs[!p].sta[k]);
124         code[j]=code[j+1]=0;
125         hs[p].add(getsta(),hs[!p].f[k]);
126     }
127 }
128 
129 LL slove()
130 {
131     int i,j,p;
132     hs[0].init();
133     hs[p=1].init();
134     hs[0].add(0,1);
135     for(i=0;i<n;i++){
136         for(j=0;j<m;j++){
137             if(g[i][j])unblock(i,j,p);
138             else block(j,p);    //p=!p优化
139             hs[p=!p].init();
140         }
141         shift(p);   //换行移位
142         hs[p=!p].init();
143     }
144     for(i=0;i<hs[!p].size;i++){
145         if(hs[!p].sta[i]==0)return hs[!p].f[i];
146     }
147 
148     return 0;
149 }
150 
151 int main()
152 {
153  //   freopen("in.txt","r",stdin);
154     int i,j;
155     LL ans;
156     char c;
157     while(~scanf("%d%d",&n,&m) && (n || m))
158     {
159         mem(g,0);
160         ex=-1;
161         for(i=0;i<n;i++){
162             for(j=0;j<m;j++){
163                 scanf(" %c",&c);
164                 g[i][j]=(c=='.');
165                 if(g[i][j])ex=i,ey=j;
166             }
167         }
168         if(ex==-1)ans=0;
169         else ans=slove();
170 
171         printf("%I64d\n",ans);
172     }
173     return 0;
174 }

 

posted @ 2013-04-01 01:10  zhsl  阅读(610)  评论(0编辑  收藏  举报