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;
}