状压dp-POJ - 3311-旅行商问题

题目链接：https://vjudge.net/contest/307763#problem/B
题目大意：

poj这题是变形：
一个送外卖的人，要将外卖全部送去所有地点再回到店里，求最短路。店在0点。
因为这里可以重复经过每个地点。所以首先用floyd求出最短路，然后再进行DP即可

我看题解基本上都是递推dp，我写的记忆化搜索。因为没有删除别人的变量定义，导致数组开小，又debug了2个小时。

递推的思路也不难
dp[s][j]表示经过s集合中的点，j是最后一个经过的点的最短路径。
刷表法：枚举i, j。j必须在集合中存在，i必须不存在。
dp[ s^(1<<i )][i+1]=min(dp[s^(1<<i)][i+1], dp[s][j+1]+dis[j+1][i+1]);//从j到i。
注意；下标不是0开始的时候，要转化(i+1)

//dfs
#include <bits/stdc++.h>
#define INF 100000000
using namespace std;
int dis[12][12];
int dp[12][(1<<12)];
int n;

int dfs(int i, int s)
{
    if(dp[i][s]!=-1)
    {
        return dp[i][s];
    }
    if(s==0)
    {
        return dp[i][0]=dis[i][0];
    }
    dp[i][s]=INF;
    for(int j=0;j<n;j++)
    {
        if(s&(1<<j))//去城市j+1
        {
            dp[i][s]=min(dp[i][s], dis[i][j+1]+dfs(j+1, s^(1<<j)));
        }
    }

    return dp[i][s];
}

int main()
{
	while(scanf("%d",&n) && n)
	{
		for(int i = 0;i <= n;++i)
        {
            for(int j = 0;j <= n;++j)
            {
                scanf("%d",&dis[i][j]);
            }
        }
        memset(dp, -1, sizeof(dp));

        //最短路
		for(int k = 0;k <= n;++k)
        {
            for(int i = 0;i <= n;++i)
            {
                for(int j = 0;j <= n;++j)
                {
                    if(dis[i][k] + dis[k][j] < dis[i][j])
                    {
                        dis[i][j] = dis[i][k] + dis[k][j];
                    }
                }
            }
        }

        cout<<dfs(0, (1<<n)-1)<<endl;

	}
	return 0;
}

//递推
#include <bits/stdc++.h>
#define INF 100000000
using namespace std;
int dis[12][12];
int dp[(1<<12)][12];
int n;

int main()
{
	while(scanf("%d",&n) && n)
	{
		for(int i = 0;i <= n;++i)
        {
            for(int j = 0;j <= n;++j)
            {
                scanf("%d",&dis[i][j]);
            }
        }
        memset(dp, -1, sizeof(dp));

        //最短路
		for(int k = 0;k <= n;++k)
        {
            for(int i = 0;i <= n;++i)
            {
                for(int j = 0;j <= n;++j)
                {
                    if(dis[i][k] + dis[k][j] < dis[i][j])
                    {
                        dis[i][j] = dis[i][k] + dis[k][j];
                    }
                }
            }
        }

        memset(dp, 3, sizeof(dp));
        for(int s=0;s<(1<<n);s++)
        {
            dp[0][0]=0;
            for(int i=0;i<n;i++)
            {
                if(!(s&(1<<i)))//i+1没有去过
                {
                    if(s==0)
                    {
                        dp[s^(1<<i)][i+1]=dis[0][i+1];
                    }
                    for(int j=0;j<n;j++)
                    {
                        if(s&(1<<j))//j+1去过
                        {
                            dp[s^(1<<i)][i+1]=min(dp[s^(1<<i)][i+1], dp[s][j+1]+dis[j+1][i+1]);//刷表法
                        }
                    }
                }
            }
        }
        int Min=INF;

        for(int i=1;i<=n;i++)
        {
            Min=min(dp[(1<<n)-1][i]+dis[i][0], Min);
        }

        cout<<Min<<endl;
	}
	return 0;
}