hdu 4026 2011上海赛区网络赛F TSP ****

没看过TSP,先mark

 //4838039 2011-10-27 23:04:15 Accepted 4026 2343MS 31044K 3143 B C++ Geners
 //状态压缩DP的TSP问题
 //优先级位运算小于判等 ， 还有各种细节各种出错
 #include <cstdio>
 #include <cstring>
 #include <stdlib.h>
 #define mabs(a) (a>0?a:-(a))

 using namespace std;

 typedef long long ll;
 const int maxn=<<;
 int gcd(int a, int b)
 {
     return b?gcd(b, a%b):a;
 }

 int g[][];
 int map[][], cnt;
 int n, m;
 ll dp[maxn][];
 struct Node {
     int x,y;
 }node[];

 void init ()
 {
     for (int i= ; i< ; ++i)
         for (int j= ; j< ; ++j)
             g[i][j]=gcd(i, j);
 }

 bool valid(int sta, int k, int pre)
 {
     int x=node[k].x,xx=node[pre].x;
     int y=node[k].y,yy=node[pre].y;

     if(map[x][y]> || map[xx][yy]>)
         return false;//如果2个点不都是0。
     int d=g[mabs(x-xx)][mabs(y-yy)];
     for (int i= ; i<=d ; ++i)
     {
         int tmpx=xx+i*(x-xx)/d;
         int tmpy=yy+i*(y-yy)/d;
         if(map[tmpx][tmpy] == )
             return false;//中间点forbidden
         if(map[tmpx][tmpy] <= )
         {
             if((sta&(<<(-map[tmpx][tmpy]))) == )
             {
                 return false;
             }//中间有未遍历的0。
         }
     }
     return true;
 }

 void DP()
 {
     memset (dp, , sizeof(dp));
     int limit=<<cnt;
     for (int i= ; i<limit ; ++i)
     {
         for (int j= ; j<cnt ; ++j)
         {
             dp[<<j][j]=1ll;//初始以每个点开始的情况
             if(i&(<<j))
             for (int k= ; k<cnt ; ++k)
             {
                 if(k==j)continue;
                 if(i^(<<j)& (<<k) )   //k是从此时出发的点，j是要去的点
                 if(valid(i,j,k))
                 {
                         dp[i][j]+=dp[i^(<<j)][k];
                 }
             }
         }
     }
     ll ans=0ll;
     for (int i= ; i<cnt ; ++i)
         ans+=dp[limit-][i];
     printf("%I64d\n", ans);

 }

 int main ()
 {
     init ();
     while (~scanf("%d%d", &n, &m))
     {
         cnt=;
         for (int i= ; i<n ; ++i)
         {
             for (int j= ; j<m ; ++j)
             {
                 scanf("%d",&map[i][j]);
                 if(map[i][j] == )
                 {
                     node[cnt].x=i;
                     node[cnt].y=j;
                     map[i][j]=-cnt;//映射到小于等于0
                     cnt++;
                 }
                 //printf(" %d \n",map[i][j]);
             }
         }
         DP();
     }
     return ;
 }