2018北大暑校acm算法训练课程 海贼王之伟大航路 状压dp

题干信息:

“我是要成为海贼王的男人！”，路飞一边喊着这样的口号，一边和他的伙伴们一起踏上了伟大航路的艰险历程。


路飞他们伟大航路行程的起点是罗格镇，终点是拉夫德鲁（那里藏匿着“唯一的大秘宝”——ONE PIECE）。而航程中间，则是各式各样的岛屿。
因为伟大航路上的气候十分异常，所以来往任意两个岛屿之间的时间差别很大，从A岛到B岛可能需要1天，而从B岛到A岛则可能需要1年。当然，任意两个岛之间的航行时间虽然差别很大，但都是已知的。
现在假设路飞一行从罗格镇（起点）出发，遍历伟大航路中间所有的岛屿（但是已经经过的岛屿不能再次经过），最后到达拉夫德鲁（终点）。假设他们在岛上不作任何的停留，请问，他们最少需要花费多少时间才能到达终点？
Input
输入数据包含多行。 
第一行包含一个整数N(2 < N ≤ 16)，代表伟大航路上一共有N个岛屿（包含起点的罗格镇和终点的拉夫德鲁）。其中，起点的编号为1，终点的编号为N。 
之后的N行每一行包含N个整数，其中，第i(1 ≤ i ≤ N)行的第j(1 ≤ j ≤ N)个整数代表从第i个岛屿出发到第j个岛屿需要的时间t(0 < t < 10000)。第i行第i个整数为0。
Output
输出为一个整数，代表路飞一行从起点遍历所有中间岛屿（不重复）之后到达终点所需要的最少的时间。
Sample Input
样例输入1：
4
0 10 20 999
5 0 90 30
99 50 0 10
999 1 2 0

样例输入2：
5
0 18 13 98 8
89 0 45 78 43 
22 38 0 96 12
68 19 29 0 52
95 83 21 24 0
Sample Output
样例输出1：
100

样例输出2：
137
Hint
提示： 
对于样例输入1：路飞选择从起点岛屿1出发，依次经过岛屿3，岛屿2，最后到达终点岛屿4。花费时间为20+50+30=100。 
对于样例输入2：可能的路径及总时间为： 
1,2,3,4,5: 18+45+96+52=211 
1,2,4,3,5: 18+78+29+12=137 
1,3,2,4,5: 13+38+78+52=181 
1,3,4,2,5: 13+96+19+43=171 
1,4,2,3,5: 98+19+45+12=174 
1,4,3,2,5: 98+29+38+43=208 
所以最短的时间花费为137 
单纯的枚举在N=16时需要14!次运算，一定会超时。

这个题暴力是不行的 阶乘的复杂度太高,不管怎么剪枝都不行

思考一下我们是否存在重复计算

暴力的想法是什么?

枚举所有情况,维护到达最终点的最大值
考虑重复情况

重复情况
选择1-2-3-4 值为120
选择1-3-2-4 值为110
这是集合为1234 终点为4的情况,以上两种情况中间的交换对后续的搜索没有任何影响,但是最小值却发生了改变
在暴力的情况下,我们无法对该情况进行剪枝,于是我们想用一种结构来进行剪枝

问题答案并不询问顺序

但是由于各点路径长度不一样,所以尽管点集不一样,到达该节点的终点不同,所获得的最短路也不同

故我们维护的应该是一个点集为s,终点为j的二维集合dp[s][j]

状态转移方程:
Dp[s][j]=min{dp[pre_s][k]+w[k][j]};
Pre_s等于s集合中抛去点j
(注意,尽管集合一样,但是走的顺序不同会导致最短路径的不同,这一点体现在k的变化,即包含1234终点为4的最短路和包含1234终点为3的最短路是分开存储的,这里需要枚举k来计算)
枚举每个k,w[k][j]是k到j的边权值

使用dp[s][j]表示经过集合s中的每个点恰好一次,且最后走的一个点是j的最佳路径长度

最终答案为min(dp[all][j]) (0<=j < N)

如何表示已经走过的点的集合?
状态压缩

点的位数为16位,故我们完全可以使用一个变量表示点集

表示某一个点走过就是(1<< j )
去掉某个点j就是s=s-(1<< j )
Dp[s][j]=min(dp[s-1<< j ] + w[k][j])//这里之前写错了

1-2-3-4(1-2-4-8)
故表示1-2-3-4的集合为(1+2+4+8)=15

这里使用的状态压缩为二进制压缩
还可以使用康托展开等等压缩手段
感兴趣的可以了解一下:
hdu1043/poj1077 一个多组输入,一个单组输入(康托展开)
Poj 3279 二进制枚举题

这里的代码是修改自戳我

如果侵犯了的话我删了重写==
主要是自己实现的差点超了时==

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<stack>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<sstream>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int MAX_N=100000;
const int INF = 0x3f3f3f3f;
#define eps 1e-7
//#define Kongxiangzhouye
//https://www.cnblogs.com/gideonzsd/p/6422208.html
//代码优化自上述博客链接,源码请戳链接查看
inline int read(){
    int a=0;
    int fh=1;
    char c=getchar();
    while(!(c<='9'&&c>='0')){
        if(c=='-')
            fh*=-1;
        c=getchar();
    }
    while(c<='9'&&c>='0'){
        a=a*10+c-'0';
        c=getchar();
    }
    return a*fh;
}
int step;//当前的最短路长度 
int minstep;
int w[16][16];
bool vis[15];
int poww[20];
int n;
int dp[15][1<<15];
int po;//po为位运算的点集 
void init(){
    po=0;
    step=0;
    minstep=INF;
    memset(dp,INF,sizeof(dp));
    memset(w,0,sizeof(w));
    memset(vis,0,sizeof(vis));
}
void dfs(int cur,int qc){
    if(cur+2==n){//点集中没有包含1 故包含n-2即为完全走过所有道路 
        step+=w[qc][n-1];//加上最后要走的路 
        minstep=min(minstep,step);
        step-=w[qc][n-1];
        return;
    }
    for(int i=1;i<n-1;i++){//不包含起点和终点的循环 
        if(!vis[i]){//判断该点是否已经被走过 
            if(step>=minstep)//如果当前步数已经大于最终的步数则直接结束(原来暴力的剪枝) 
                continue;
            if((dp[i][po+poww[i-1]]<=step+w[qc][i]))//如果已走集合且终点为i-1的最短路已经小于当前情况 
                continue;
            po+=poww[i-1];
            step+=w[qc][i];
            dp[i][po]=step;
            vis[i]=true;
            dfs(cur+1,i);
            po-=poww[i-1];
            vis[i]=false;
            step-=w[qc][i];
        }
    }
}
int main(){
    ios::sync_with_stdio(false);
    #ifdef Kongxiangzhouye
        freopen("in.txt","r",stdin);
        freopen("out.txt","w",stdout);
    #endif
    for(int i=0; i<14; i++)//预处理位运算的结果以减少重复计算 
        poww[i]=1<<i;
    while (scanf("%d",&n)!=EOF){
        init();
        for(int i=0; i<n; i++){
            for(int j=0; j<n; j++){
                w[i][j]=read();
            }
        }
        dfs(0,0);
        cout<<minstep<<endl;
    }
    return 0;
}