POJ 3311 旅行商问题 状态压缩 dp

最近做 DP 做的起劲啊。。。

旅行商问题的状态压缩 dp[i][j] 表示 遍历 i 中所有城市一遍 并且现在在 j 的最短距离是多少。

第一个点可以遍历两次所以先不设置为1，并且这个状态是合法状态，这个点曾经是最困扰我的。

#include <stdio.h>
#include <iostream>
#include <queue>
#include <algorithm>
#include <map>
#include <vector>
#include <cmath>
#include <string.h>
#include <stdlib.h>
#include <time.h>
#include <fstream>

using namespace std;

#define READ freopen("acm.in","r",stdin)
#define ll long long
#define PII pair<int,int>
#define PDI pair<double,int>
#define PDD pair<double,double>
#define MPI map<int,int>::iterator 
#define fst first
#define sec second
#define MS(x,d) memset(x,d,sizeof(x))
#define INF 0x3f3f3f3f
#define ALL(x) x.begin(),x.end()
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define MAX 500000
#define ROOT 0,n-1,1
#define PB push_back
#define FOR(a,b,c) for(int a=b;a<c;a++)
#define MOD 706679
#define keyTree (ch[ ch[root][1] ][0])

int dist[50][50];
int dp[1<<12][20];
int main()
{
    int n;
    READ;
    while(scanf("%d",&n)!=EOF&&n!=0)
    {
        n++;
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                scanf("%d",&dist[i][j]);
        for(int k=0;k<n;k++)
            for(int i=0;i<n;i++)
                for(int j=0;j<n;j++)
                {
                    if(dist[i][j]>dist[i][k]+dist[k][j])
                        dist[i][j]=dist[i][k]+dist[k][j];
                }
        for(int i=0;i<(1<<12);i++)
            for(int j=0;j<20;j++)
                dp[i][j]=INF;
        dp[0][0]=0;
        for(int i=1;i<(1<<n);i++)
        {
            for(int j=0;j<n;j++)
            {
                if(((1<<j)&i))
                {
                    for(int k=0;k<n;k++)
                    {
                        if(!(i-(1<<j)) || ((1<<k)&(i-(1<<j))) )
                        {
                            dp[i][j]=min(dp[i][j],dp[i-(1<<j)][k]+dist[k][j]);
                        }
                    }
                }
            }
        }
        cout<<dp[(1<<n)-1][0]<<endl;
    }
    return 0;
}