【HNOI2004】邮递员 （插头DP一条回路）

【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;
}