题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1964
题目大意:单回路最小花费。
题目思路:插头dp。
#include<stdio.h>
#include<string.h>
#define Max 1001000
#define inf 0x3f3f3f3f
#define Hash 3007
int min(int a,int b)
{
return a<b?a:b;
}
int n,m,ex,ey,cur,mp[15][15][2],a1,a2,a3,b1,b2,b3;
int stack[15],f[15];
char str[25][25];
struct node
{
int size,next[Max],p[Hash],state[Max];
int f[Max];
inline void init()
{
memset(p,-1,sizeof(p));
size=0;
}
inline void push(int st,int val)
{
// printf("akkk\n");
int i,u=st%Hash;
for(i=p[u];i!=-1;i=next[i])
{
if(state[i]==st)
{
f[i]=min(f[i],val);
return;
}
}
state[size]=st;f[size]=val;
next[size]=p[u];
p[u]=size++;
}
}dp[2];
inline void decode(int st)
{
int top=0;
for(int i=0;i<=m;i++)
{
if((st&3)==1)
stack[top++]=i;
else if((st&3)==2)
{
f[stack[top-1]]=i;
f[i]=stack[top-1];
top--;
}
st>>=2;
}
}
inline void shift()
{
for(int k=0;k<dp[cur].size;k++)
dp[cur^1].push(dp[cur].state[k]<<2,dp[cur].f[k]);
}
inline void dpblank(int i,int j)
{
int k,left,up;
for(k=0;k<dp[cur].size;k++)
{
int st=dp[cur].state[k];
left=st&a3;
up=st&b3;
if(left&&up)
{
if(left==a1&&up==b2)
{
if(i==n&&j==m)
dp[cur^1].push(st^left^up,dp[cur].f[k]);
}
else
{
if(left==a2&&up==b1)
dp[cur^1].push(st^left^up,dp[cur].f[k]);
if(left==a1&&up==b1)
{
decode(st);
dp[cur^1].push(st^left^up^(3<<(2*f[j])),dp[cur].f[k]);
}
if(left==a2&&up==b2)
{
decode(st);
dp[cur^1].push(st^left^up^(3<<(2*f[j-1])),dp[cur].f[k]);
}
}
}
else if(left)
{
if(mp[i][j][0]!=-1)
dp[cur^1].push((st^left)|(left<<2),dp[cur].f[k]+mp[i][j][0]);
if(mp[i][j][1]!=-1)
dp[cur^1].push(st,dp[cur].f[k]+mp[i][j][1]);
}
else if(up)
{
if(mp[i][j][0]!=-1)
dp[cur^1].push(st,dp[cur].f[k]+mp[i][j][0]);
if(mp[i][j][1]!=-1)
dp[cur^1].push((st^up)|(up>>2),dp[cur].f[k]+mp[i][j][1]);
}
else if(mp[i][j][0]!=-1&&mp[i][j][1]!=-1)
dp[cur^1].push(st|a1|b2,dp[cur].f[k]+mp[i][j][0]+mp[i][j][1]);
}
}
inline void solve()
{
int i,j,k;
int ans=inf;
cur=0;
dp[0].init();
dp[0].push(0,0);
for(i=1;i<=n;i++)
{
dp[cur^1].init();
shift();
cur^=1;
b1=1;b2=2;
for(j=1;j<=m;j++)
{
a1=b1;a2=b2;a3=a1|a2;
b1<<=2;b2<<=2;b3=b1|b2;
dp[cur^1].init();
dpblank(i,j);
cur^=1;
}
}
for(k=0;k<dp[cur].size;k++)
ans=min(ans,dp[cur].f[k]);
printf("%d\n",ans);
}
int main()
{
int i,j;
int t;
scanf("%d",&t);
while(t--)
{
//ex=0;
// n=2*n+1;
// m=2*m+1;
scanf("%d%d",&n,&m);
getchar();
memset(mp,-1,sizeof(mp));
for(i=1;i<=2*n+1;i++)
{
gets(str[i]+1);
}
for(i=1;i<=n;i++)
for(j=1;j<=m;j++)
{
if(j<m)
{
mp[i][j][0]=str[2*i][2*j+1]-'0';
}
if(i<n)
mp[i][j][1]=str[2*i+1][2*j]-'0';
// printf("i %d j %d %d %d\n",i,j,mp[i][j][0],mp[i][j][1]);
}
solve();
}
return 0;
}