【HNOI2004】邮递员 |
Time Limit:10000MS Memory Limit:65536K
Total Submit:47 Accepted:21
Case Time Limit:1000MS
Description
Smith在P市的邮政局工作,他每天的工作是从邮局出发,到自己所管辖的所有邮筒取信件,然后带回邮局。
他所管辖的邮筒非常巧地排成了一个m*n的点阵(点阵中的间距都是相等的)。左上角的邮筒恰好在邮局的门口。
Smith是一个非常标新立异的人,他希望每天都能走不同的路线,但是同时,他又不希望路线的长度增加,他想知道他有多少条不同的路线可走。
【任务描述】
你的程序需要根据给定的输入,给出符合题意的输出:
输入包括点阵的m和n的值;
你需要根据给出的输入,计算出Smith可选的不同路线的总条数;
Input
输入文件postman.in只有一行。包括两个整数m, n(1 <= m <= 10, 1 <= n <= 20),表示了Smith管辖内的邮筒排成的点阵。
Output
只有一行,只有一个整数,表示Smith可选的不同路线的条数。
Sample Input
2 2
说明:该输入表示,Smith管辖了2*2的一个邮筒点阵。
Sample Output
2
Hint
Source
xinyue
题目:http://mail.bashu.cn:8080/bs_oj/showproblem?problem_id=2502
分析:这题可以用插头DP一条回路来做,如果该图存在回路的话,直接求一条回路的数量再乘以2,因为反向走也算一种,如果不存在回路的话,一开始我以为直接求一条路径的数量,这样太复杂了,不过题目要求使得路径最短。。。这样的话就好办多了
好吧,这种做法没有证明,不知道对不对,只做参考,被误解的同学抱歉了T_T
对于一个图,如果不存在回路的话必定能使得终点在三个位置,路径才会最短,如下图红色点
其实这三个点为终点的路径总数是一样的,那么我们就可以只求一个再乘3,就用最上面那个吧,然后由于起点和终点都确定了,所以有是一条回路问题了,不过貌似数据要用大数(高精度),我的代码就没写了,所以最后一组数据过不了。。。
代码:
#include<cstdio>
#include<cstring>
using namespace std;
const int mm=15511;
struct data
{
int s[mm],h[mm],p[mm],t;
long long d[mm];
inline int hash(int x,long long dd)
{
int i,c=x%mm;
for(i=h[c];i>=0;i=p[i])
if(s[i]==x)
{
d[i]+=dd;
return i;
}
s[t]=x,p[t]=h[c],h[c]=t,d[t]=dd;
return t++;
}
void clear()
{
t=0,memset(h,-1,sizeof(h));
}
}f[2];
int i,j,g1,g2,u,k,n,m,a,b,a1,a2,b1,b2;
bool g[33][13];
inline int eat(int s,int a,int b,int c,bool f)
{
int n=1,x;
while(n)
{
if(f)a<<=2,b<<=2,c<<=2;
else a>>=2,b>>=2,c>>=2;
x=s&c;
if(x==a)++n;
if(x==b)--n;
}
return (s^b)|a;
}
inline bool ok(int c)
{
if(c==1)return g[i+1][j];
if(c==2)return g[i][j+1];
if(c==3)return g[i+1][j]&&g[i][j+1];
return 1;
}
inline void work(int S,long long d)
{
int x,y;
x=a&S,y=b&S;
if(x==0&&y==0)
{
if(ok(3))f[g2].hash(S|a1|b2,d);
}
else if(x==0&&y==b1)
{
if(ok(2))f[g2].hash(S,d);
if(ok(1))f[g2].hash((S^b1)|a1,d);
}
else if(x==a1&&y==0)
{
if(ok(1))f[g2].hash(S,d);
if(ok(2))f[g2].hash((S^a1)|b1,d);
}
else if(x==a1&&y==b1)f[g2].hash(eat(S^a1^b1,b1,b2,b,1),d);
else if(x==0&&y==b2)
{
if(ok(2))f[g2].hash(S,d);
if(ok(1))f[g2].hash((S^b2)|a2,d);
}
else if(x==a2&&y==0)
{
if(ok(1))f[g2].hash(S,d);
if(ok(2))f[g2].hash((S^a2)|b2,d);
}
else if(x==a2&&y==b2)f[g2].hash(eat(S^a2^b2,a2,a1,a,0),d);
else if(x==a2&&y==b1)f[g2].hash(S^a2^b1,d);
}
inline long long DP()
{
long long ss;
if(n==1||m==1)return 1;
f[0].clear();
if((n*m)&1)f[0].hash(33,1),ss=3;
else f[0].hash(0,1),ss=2;
for(g1=1,g2=i=0;i<=n;++i)
{
for(k=0;k<f[g2].t;++k)f[g2].s[k]<<=2;
a=3,b=3<<2,a1=1,a2=2,b1=1<<2,b2=2<<2;
for(j=0;j<m;a<<=2,b<<=2,a1<<=2,a2<<=2,b1<<=2,b2<<=2,++j)
if(g[i][j])for(g1=!g1,g2=!g2,f[g2].clear(),k=0;k<f[g1].t;++k)
{
if(i==n-1&&j==m-1)
{
if((f[g1].s[k]&(a1|b2))==(a1|b2))return f[g1].d[k]*ss;else continue;
}
work(f[g1].s[k],f[g1].d[k]);
}
}
return 0;
}
int main()
{
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
memset(g,0,sizeof(g));
scanf("%d%d",&n,&m);
if(n<m)
{
int tt=n;
n=m;
m=tt;
}
for(i=0;i<n;++i)
for(j=0;j<m;++j)
g[i][j]=1;
printf("%I64d\n",DP());
return 0;
}