hdu 5450 Traversal(状压dp)

题意比较迷，第一句话的意思是那k个数不存在不用管，出题人的英语.....

但是这个题还是比较难的。

每个点只能到差值1, p, p+2的点，要求每次走数量大于等于3的环，不能重复走完所有点。走完所有点形成环其实相当于，每个点出度入度两次，要超过3个点意味着不能两两成环，更不能自环。

由于最多只能到达p+2之外的点，所以每个点一定要在p+2的距离只能形成一个环。可以维护当前点i到之前i-p-2之间点的度，考虑向度小于2的点连边然后转移，对于之前p+2个点的度小于2的状态是无效的，因为i+1不可能连向i-p-2。

然后复杂度就是O(n*3^n)。

#include <bits/stdc++.h>
using namespace std;
const int mod=1e4+7;
bool book[105];
int dp[3][531541];
int cur[15], nex[15], N, n, p;
int dx[5];
int encode(int *a, int n)
{
    int res=0;
    for(int i=0; i<n; i++)res=res*3+a[i];
    return res;
}
void decode(int *a, int s, int n)
{
    for(int i=n-1; i>=0; i--)
    {
        a[i]=s%3;
        s/=3;
    }
}
void move(int *a, int n)
{
    for(int i=n; i>=1; i--)
    {
        a[i]=a[i-1];
    }
    a[0]=0;
}

void init()
{
    int x, i, k;
    memset(book, 0, sizeof book);
    scanf("%d%d", &N, &p);
    scanf("%d", &k);
    dx[0]=1, dx[1]=p, dx[2]=p+2;
    for(i=0; i<k; i++)
    {
        scanf("%d", &x);
        book[x]=1;
    }
}
void add(int &u, int v){u=(u+v)%mod;}

int solve()
{
   int i, j, e, tar;
   n=p+2;
   for(i=0; i<n; i++)cur[i]=2;
   e=encode(cur, n);
   memset(dp, 0, sizeof dp);
   dp[0][e]=1;
   for(i=0; i<N; i++)
   {
        for(j=0; j<=e; j++)
        {
            if(dp[0][j]==0)continue;
            decode(cur, j, n);
            for(int l=0; l<n; l++)nex[l]=cur[l];
            move(nex, n);

                if(book[i+1]==1)
                {
                    nex[0]=2;
                    tar=encode(nex, n);
                    if(nex[n]==2)add(dp[1][tar], dp[0][j]);
                    continue;
                }

            nex[0]=0;
            tar=encode(nex, n);
            if(nex[n]==2)add(dp[1][tar], dp[0][j]);

            nex[0]=1;
            for(int l=0; l<3; l++)
            {
                int &v=nex[dx[l]];
                if(v<2)
                {
                    v++;
                    tar=encode(nex, n);        
                    if(nex[n]==2)add(dp[1][tar], dp[0][j]);
                    v--;
                }
            }

            nex[0]=2;
            for(int l=0; l<3; l++)
            {
                for(int o=l+1; o<3; o++)
                {
                    int &u=nex[dx[l]], &v=nex[dx[o]];
                    if(u<2 && v<2)
                    {
                    u++, v++;
                    tar=encode(nex, n);
                    if(nex[n]==2)add(dp[1][tar], dp[0][j]);
                    u--, v--;
                    }
                }
            }

        }
        for(j=0; j<=e; j++){dp[0][j]=dp[1][j];dp[1][j]=0;}
   }
   return dp[0][e];
}

int main()
{
    int e=1, t, i, j;
    cin>>t;
    while(t--)
    {
        init();
        printf("Case #%d: %d\n", e++, solve());
    }
}