3122 奶牛代理商 VIII

时间限制: 3 s

空间限制: 256000 KB

题目等级 : 大师 Master

题目描述 Description

小徐是USACO中国区的奶牛代理商，专门出售质优价廉的“FJ"牌奶牛。

有一天，她的奶牛卖完了，她得去美国进货。

她需要去N个奶牛农场询问价格（小徐是个认真的人，买东西一定要货比三家）。

给你一个邻接矩阵，表示N个农场间的路径长度，求小徐最少走多少路。(从农场1出发，最后回到出发点买）

输入描述 Input Description

邻接矩阵

输出描述 Output Description

答案（见描述）

样例输入 Sample Input

0 1 2

3 0 10

2 0 0

样例输出 Sample Output

数据范围及提示 Data Size & Hint

N<=15，路径长度<=1000

TSP

分类标签 Tags 点此展开

状压DP好恶心啊。。

这道题的关键点有两个，

1.走过所有的点

2.最短路径

第2个最短路径比较好解决，n<=16的话，，一遍Floyd就可以

但是第一个条件，要走过所有的点。

我们可以考虑用状态压缩的方法来实现

我们可以用一个二进制的字符串表示这个点是否走过，比如说110表示已经走过了1和2这两个点，第3个点还没有经过

代码思路：

设置一个dp数组,dp[now][j]表示在now状态下，到达点j所需要的花费

首先我们需要暴力枚举i和j两点，来求最短距离

其次，我们还需要枚举一个能够包揽所有状态的变量now，来记录每一个能够到达的状态

当状态now可以到达j的话，那么说明我们可以通过这个状态到达i（i和j之间必定有路径）

最后枚举每个点，取一下最小值就可以

细节问题：

1.跑floyd的时候不要预先设定最大值，因为每两个点（不相同）之间必定有边相连

2.dp数组的第一位必须要开的足够大，最小是2^16，因为第一维记录的是状态而不是大小

3.now<=(1<<n)-1:

　　当n==3时，1<<3 == 2^3 == 8 == 1000

　　1000-1=111正好是三个点都到达的理想情况

4.now&(1<<(j-1))

　　j-1是为了不超边界且枚举出所有情况

　　首先要明确，1<<(j-1)得到的一定是一个2^x的数，转换成二进制一定是1+000.....的形式

　　那么当now&(1<<(j-1))有值时，就说状态now是可以到达j点的

5.now|(1<<(i-1))

　　在进行这个运算的时候，now一定是满足now&(1<<(j-1))！=0的（程序满足顺序执行）

　　那么通过now状态一定可以到达j点

　　而且1<<(i-1)也一定是一个2^x的数

　　所以now|(1<<(i-1))就是一个可以通过j到达i产生的新状态

　举个栗子：

　now=110010

　i=100

　那么结果就是110110

代码：

 #include<iostream>

 #include<cstdio>

 #include<cstring>

 #include<cmath>

 using namespace std;

 const int MAXN=;

 int read(int & n)

 {

     char p='+';int x=;

     while(p<''||p>'')

         p=getchar();

     while(p>=''&&p<='')

     x=x*+p-,p=getchar();

     n=x;

 }

 int dp[MAXN*][MAXN];

 int dis[MAXN][MAXN];

 int n;

 void floyed()

 {

     for(int i=;i<=n;i++)

         for(int j=;j<=n;j++)

             for(int k=;k<=n;k++)

                 if(i!=j&&j!=k&&i!=k)

                     dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);

 }

 void zhuangya()

 {

     /*for(int i=1;i<=n;i++)

     {

         for(int j=1;j<=n;j++)

         {

             cout<<dis[i][j]<<" ";

         }

         cout<<endl;

     }*/

     memset(dp,0xf,sizeof(dp));

     dp[][]=;

     for(int now=;now<=(<<n)-;now++)

         for(int i=;i<=n;i++)

             for(int j=;j<=n;j++)

                 if((now&(<<(j-)))&&i!=j)

                 {

                     dp[now|(<<(i-))][i]=

                         min

                         (

                             dp[now|(<<(i-))][i],

                             dp[now][j]+dis[j][i]

                         );

                 }

     int ans=0x7ffffff;

     for(int i=;i<=n;i++)

     {

         ans=min(ans,dp[(<<n)-][i]+dis[i][]);

     }

     cout<<ans;

 }

 int main()

 {

     read(n);

     for(int i=;i<=n;i++)

         for(int j=;j<=n;j++)

             read(dis[i][j]);

     floyed();

     zhuangya();

     return ;

 }