poj 1739 Tony's Tour 插头dp
第一道插头dp，还是非常艰难的，做法的话cdq论文中已经讲得很清楚了，用的是最小表示法，但是这个题目因为障碍的存在导致细节处理上面更加麻烦。
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=9,maxm=3e3+9;;
char a[maxn][maxn];
bool first[9][9];
int n,m;
int cag(int txt)
{
int ss[11];
int lon=-1;
while(txt)
{
ss[++lon]=txt%8;
txt/=8;
}
short fla[11],top=0;
memset(fla,0,sizeof(fla));
for(int i=0;i<=lon;i++)
if(ss[i])
{
if(!fla[ss[i]])
{
fla[ss[i]]=++top;
ss[i]=top;
}
else
ss[i]=fla[ss[i]];
}
txt=0;
int tt=1;
for(int i=0;i<=lon;i++)
{
txt+=ss[i]*tt;
tt*=8;
}
return txt;
}
struct
{
int hash[maxm],lon;
struct
{
int key,next;
long long sum;
}data[maxm];
void clear()
{
memset(hash,-1,sizeof(hash));
lon=0;
}
void push(int key,long long sum)
{
key=cag(key);
// cout<<key<<endl;
int tmp=key%maxm;
for(int k=hash[tmp];k!=-1;k=data[k].next)
{
if(data[k].key==key)
{
data[k].sum+=sum;
return ;
}
}
data[++lon].key=key;
data[lon].sum=sum;
data[lon].next=hash[tmp];
hash[tmp]=lon;
}
}dp[2];
int getbit(int tmp,int k)
{
tmp=(tmp&(7<<k*3));
return tmp>>k*3;
}
int cal(int tmp)
{
int ans=0;
while(tmp)
{
ans=max(ans,tmp%8);
tmp/=8;
}
return ans+1;
}
int ss[11];
int mer(int tmp,int t,int s)
{
int txt=tmp,lon=-1;
while(txt)
{
ss[++lon]=txt%8;
txt/=8;
}
// cout<<endl;
// printf("%d %d\n",t,s);
// for(int i=0;i<=lon;i++)
// printf("%d",ss[i]);
// cout<<endl;
int a=ss[t],b=ss[s];
int c=max(a,b);
int d=min(a,b);
ss[t]=0;
ss[s]=0;
for(int i=0;i<=lon;i++)
if(ss[i]==c)
ss[i]=d;
// for(int i=0;i<=lon;i++)
// printf("%d",ss[i]);
// cout<<endl;
txt=0;
int tt=1;
for(int i=0;i<=lon;i++)
{
txt+=ss[i]*tt;
tt*=8;
}
return txt;
}
int ff(int tmp)
{
if(tmp==0) return tmp;
int txt=tmp,lon=-1;
while(txt)
{
ss[++lon]=txt%8;
txt/=8;
}
int i,j;
for(i=0;i<=lon;i++)
if(ss[i])
break;
for(j=lon;j>=0;j--)
if(ss[i])
break;
ss[i]=ss[j];
txt=0;
int tt=1;
for(int i=0;i<=lon;i++)
{
txt+=ss[i]*tt;
tt*=8;
}
return txt;
}
void solve()
{
int now=0,pre=1;
dp[now].clear();
dp[pre].clear();
for(int i=n;i>=1;i--)
for(int j=1;j<=m;j++)
{
now=now^1;
pre=pre^1;
dp[now].clear();
// cout<<dp[pre].lon<<endl;
if(a[i][j]=='#')
{
for(int k=1;k<=dp[pre].lon;k++)
{
int tmp=dp[pre].data[k].key;
long long sum=dp[pre].data[k].sum;
int up=getbit(tmp,j);
int r=getbit(tmp,0);
if(!up&&!r)
dp[now].push(tmp,sum);
}
continue;
}
if(i==n&&j==1)
{
dp[pre].push(1,1);
}
for(int k=1;k<=dp[pre].lon;k++)
{
int tmp=dp[pre].data[k].key;
long long sum=dp[pre].data[k].sum;
int up=getbit(tmp,j);
int r=getbit(tmp,0);
// for(int p=0;p<=m;p++)
// printf("%d",getbit(tmp,p));
// printf("\n");
// if(i==1&&j==m)
// printf("%d %d %d %d %d %d\n\n",i,j,up,r,tmp,sum);
if(tmp==0) continue;
if(up&&r)
{
if(up!=r)
dp[now].push(mer(tmp,j,0),sum);
else if(first[i][j])
dp[now].push(mer(tmp,j,0),sum);
}
else if(up&&!r)
{
if(j!=m)
dp[now].push(tmp-(up<<j*3)+up,sum);
// else if(j==m&&i==n)
// dp[now].push(tmp-(up<<j*3),sum);
if(j!=m||i!=n)
if(i!=1)
dp[now].push(tmp,sum);
}
else if(!up&&r)
{
if(j!=m||i!=n)
if(i!=1)
dp[now].push(tmp+(r<<j*3)-r,sum);
if(j!=m)
dp[now].push(tmp,sum);
else if(j==m&&i==n)
{
int txt=ff(tmp-r);
dp[now].push(txt,sum);
}
}
else if(!up&&!r)
{
if(j!=m&&i!=1)
dp[now].push(tmp+((cal(tmp))<<j*3)+cal(tmp),sum);
else if(j==m&&i==n)
{
int txt=ff(tmp+((cal(tmp))<<j*3));
dp[now].push(txt,sum);
}
}
}
}
long long ans=0;
for(int i=1;i<=dp[now].lon;i++)
if(dp[now].data[i].key==0)
ans+=dp[now].data[i].sum;
cout<<ans<<endl;
}
int main()
{
// freopen("in.txt","r",stdin);
while(scanf("%d %d",&n,&m)!=EOF)
{
if(n==0&&m==0) break;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
scanf("%c",&a[i][j]);
if(a[i][j]==' '||a[i][j]=='\n')
j--;
}
memset(first,0,sizeof(first));
for(int i=1;i<=n;i++)
for(int j=m;j>=1;j--)
if(a[i][j]=='.')
{
first[i][j]=1;
i=n+1;
break;
}
solve();
}
return 0;
}