题目链接:郊区春游
思路:首先先用floyd求一下各个点之间的最短路,那么接下来就是找出一条最短路将R个点串起来。最简单的就是枚举每一种排列然后求最小值,那么时间复杂度就是15!,也就是1307674368000,这样肯定超时,那么我们就是用状压的方法来做,dp [i] [j]:表示以i为当前的路径状态,以j为终点的最小距离。(i在二进制表示下1表示该位置已经走过,0表示该位置未走过)通过枚举每一种状态(及从1到1<<r-1),然后找出当前走过的点,再去找出没有走过的点,更新一下就行了。
状态转移方程:
dp[i+(1<<(k-1))][k]=min(dp[i+(1<<(k-1))][k],dp[i][j]+dist[a[j]][a[k]]);
核心代码:
for(int i=1;i<=len;i++){
for(int j=1;j<=r;j++){
if(i>>(j-1)&1){ //找出当前已经走过的点
// cout<<"j:"<<j<<endl;
for(int k=1;k<=r;k++){
if(i>>(k-1)&1) continue; //找出没有走过的点
dp[i+(1<<(k-1))][k]=min(dp[i+(1<<(k-1))][k],dp[i][j]+dist[a[j]][a[k]]); //用以j为终点,去更新以k为终点。
}
}
}
}
代码:
#include<iostream>
#define IOS ios::sync_with_stdio(false);
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<deque>
using namespace std;
typedef long long ll;
typedef double dd;
typedef pair<int, int> PII;
const int N=210;
int dx[4]={-1,0,0,1};
int dy[4]={0,1,-1,0};
int lowbit(int x){
return x&-x;
}
int n,m,r;
int a[20];
int dp[1<<16][20];
int dist[N][N];
void fload(){
for(int k=1;k<=n;k++){
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
dist[i][j]=min(dist[i][j],dist[i][k]+dist[k][j]);
}
}
}
}
int main(){
cin>>n>>m>>r;
memset(dist,0x3f,sizeof dist);
for(int i=1;i<=r;i++) cin>>a[i];
for(int i=1;i<=m;i++){
int x,y,w;
cin>>x>>y>>w;
dist[x][y]=min(dist[x][y],w);
dist[y][x]=min(dist[y][x],w);
}
for(int i=1;i<=n;i++) dist[i][i]=0;
fload();
int len=(1<<r)-1;
memset(dp,0x3f,sizeof dp);
for(int i=1;i<=r;i++){
dp[1<<(i-1)][i]=0;
}
for(int i=1;i<=len;i++){
for(int j=1;j<=r;j++){
if(i>>(j-1)&1){
// cout<<"j:"<<j<<endl;
for(int k=1;k<=r;k++){
if(i>>(k-1)&1) continue;
dp[i+(1<<(k-1))][k]=min(dp[i+(1<<(k-1))][k],dp[i][j]+dist[a[j]][a[k]]);
}
}
}
}
int ans=1e9;
for(int i=1;i<=r;i++) ans=min(ans,dp[len][i]);
cout<<ans;
return 0;
}