OpenJudge 百练 4124:海贼王之伟大航路

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

路飞他们伟大航路行程的起点是罗格镇，终点是拉夫德鲁（那里藏匿着“唯一的大秘宝”——ONE PIECE）。而航程中间，则是各式各样的岛屿。

因为伟大航路上的气候十分异常，所以来往任意两个岛屿之间的时间差别很大，从A岛到B岛可能需要1天，而从B岛到A岛则可能需要1年。当然，任意两个岛之间的航行时间虽然差别很大，但都是已知的。

现在假设路飞一行从罗格镇（起点）出发，遍历伟大航路中间所有的岛屿（但是已经经过的岛屿不能再次经过），最后到达拉夫德鲁（终点）。假设他们在岛上不作任何的停留，请问，他们最少需要花费多少时间才能到达终点？

输入
输入数据包含多行。
第一行包含一个整数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。
输出
输出为一个整数，代表路飞一行从起点遍历所有中间岛屿（不重复）之后到达终点所需要的最少的时间。
样例输入

样例输入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

样例输出

样例输出1：
100
样例输出2：
137

提示

提示：
对于样例输入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!次运算，一定会超时。

思路

单纯的用回溯枚举一定会超时
可以使用状态压缩DP记录曾经到达过的岛屿所经历的时间
利用记忆化搜索 进行最优性剪枝

使用状态压缩DP涉及到位运算

位运算基础

位运算

按位与运算符（&）
参加运算的两个数据，按二进制位进行“与”运算。

运算规则：0&0=0; 0&1=0; 1&0=0; 1&1=1;

  即：两位同时为“1”，结果才为“1”，否则为0

例如：3&5 即 0000 0011& 0000 0101 = 00000001 因此，3&5的值得1。

按位或运算符（|）
参加运算的两个对象，按二进制位进行“或”运算。

运算规则：0|0=0； 0|1=1； 1|0=1； 1|1=1；

 即 ：参加运算的两个对象只要有一个为1，其值为1。

例如:3|5　即 00000011 | 0000 0101 = 00000111 因此，3|5的值得7。

异或运算符（^）
参加运算的两个数据，按二进制位进行“异或”运算。

运算规则：0^0=0； 0^1=1； 1^0=1； 1^1=0；

即：参加运算的两个对象，如果两个相应位为“异”（值不同），则该位结果为1，否则为0。

例如：1011 ^ 0110 = 1101

相当于二进制位的无进位相加

左移运算符（<<）
将一个运算对象的各二进制位全部左移若干位（左边的二进制位丢弃，右边补0）。

例：a = a<< 2将a的二进制位左移2位，右补0，

左移1位后a = a *2;

若左移时舍弃的高位不包含1，则每左移一位，相当于该数乘以2。

右移运算符（>>）
将一个运算对象的各二进制位全部右移若干位（右边的二进制位丢弃，符号位不变，左边补0）。

例：a = a>> 2将a的二进制位右移2位，左补0，

右移1位后a = a /2;

若右移时舍弃的高位不包含1，则每右移一位，相当于该数除以2。

由于每个岛屿只有到达和未到达两种状态 可以对应二进制的1和0
所以使用一个多位的二进制数表示到达过的岛屿所形成的集合
例如：
0001 表示到达过岛屿1
0011 表示到达过岛屿1和2
…
集合

int toltalSet;/*已经经过的岛屿所形成的集合*/

更新集合

toltalSet = toltalSet | (1<<(i-1));

1<<(i-1)表示到达岛屿i
和 totalSet进行与运算 可以将totalSet更新

DFS
使用minT记录每个状态的最优解 大于等于这个值就不搜了

int minT[MaxSet][MaxN];

最优性剪枝

if(totalTime+road[pos][i] >= minT[toltalSet|(1<<(i-1))][i])/*最优性剪枝*/
			continue;

完整代码

#include<bits/stdc++.h>
using namespace std;
const int MaxSet = 1 << 17;
const int MaxN = 25;
int N;/*岛屿个数*/
int road[35][35];/*邻接矩阵*/
int minTime = 1 << 30;/*记录最短的时间*/
int minT[MaxSet][MaxN]; /*minT[s][i]表示已经经过集合s 且最后一个岛屿是i*/
int totalTime;//当前经过的时间
int toltalSet;/*已经经过的岛屿所形成的集合*/
void dfs(int pos, int cnt) { /*上一次到达的岛屿 已经到达的岛屿个数*/
	if(totalTime>=minTime)/*最优性剪枝*/
		return ;
	if (cnt == N - 1) {/*已经经过N-1个岛屿*/
		minTime = min(minTime, totalTime + road[pos][N]);/*更新*/
		return ;
	}
	for (int i = 2; i < N; ++i) {
		if ((toltalSet & (1 << (i - 1))) != 0)/*判断这个岛屿是否已经来过*/
			continue;
		if(totalTime+road[pos][i] >= minT[toltalSet|(1<<(i-1))][i])/*最优性剪枝*/
			continue;
		minT[toltalSet|(1<<(i-1))][i] = totalTime+road[pos][i];/*更新*/
		totalTime += road[pos][i];/*更新花费的时间*/
		toltalSet = toltalSet | (1<<(i-1));/*更新经过到达过的岛屿所形成的集合*/
		dfs(i,cnt+1);
		toltalSet = toltalSet & (~(1<<(i-1)));/*回溯 消除到达岛屿i的影响*/
		totalTime -= road[pos][i];
	}
}

int main() {

	cin >> N;
	for (int i = 1; i <= N; ++i)
		for (int j = 1; j <= N; ++j)
			cin >> road[i][j];/*读入*/
	for (int i = 0; i < MaxSet; ++i)
		for (int j = 0; j < MaxN; ++j)
			minT[i][j] = 1 << 30;/*初始化*/
	toltalSet = 1;/*初始化 到达过岛屿1*/
	dfs(1,1);/*深搜*/
	cout << minTime << endl;
	return 0;
}