狀態壓縮(TSP問題)

POJ5067

題意

​ 有一個n*m的矩陣上存在一些石頭,接下來要從起點出發把這些石頭拿上並回到起點,問最小需要多少步。

思路

​ 用sum表示除了(1,1)這個點外其他有石頭的點的個數,用x數組和y數組表示這些點的座標的位置,在這裏我用二進制下的每個位置表示點的位置,將最後的位置來表示起點,0 - n-1分別表示其他點的位置,然後開始暴力枚舉二進制下的每個數字,注意初始化dp[0] [n]=0.

代碼

#include<bits/stdc++.h>
using namespace std;
int x[20],y[20];
int ma[60][60];
int dp[1<<12][15];
int dis(int i,int j){
	return abs(x[i]-x[j])+abs(y[i]-y[j]);
}
int main(){
	int n,m;
	while(~scanf("%d%d",&n,&m)){
		int sum=0;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++){
				scanf("%d",&ma[i][j]);
				if(i==1&&j==1) continue;
				if(ma[i][j]){
					x[sum]=i;
					y[sum]=j;
					sum++;
				}
			}
		x[sum]=1,y[sum]=1;
		n=sum;
		memset(dp,0x3f,sizeof(dp));
		dp[0][n]=0;
		for(int st=0;st<(1<<n);st++){
			for(int i=0;i<=n;i++){
				dp[st][n]=min(dp[st][n],dp[st][i]+dis(i,n));
				for(int j=0;j<n;j++){
					dp[st|1<<j][j]=min(dp[st|1<<j][j],dp[st][i]+dis(i,j));
				}
			}
		}
		printf("%d\n",dp[(1<<n)-1][n]);
	}
}

HDU3311

題意

​ 給n個地點,然後輸入一個(n+1)*(n+1)的矩陣,第i行第j列表示第i個地點到第j個地點的時間,要從第0個地點用最短的時間走遍1-n個地點,並回到原點。

思路

​ 先用Floyd最短路徑算法求出任意倆點之間的最短時間,接下來就是一個TSP板子,記住最後再讓他回到原點就行了。

代碼

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std ;
int d[20][20];
int dp[1<<12][20];
int main(){
    int n;
    while(~scanf("%d",&n)){
        if(n==0) break;
        for(int i=0;i<=n;i++){
            for(int j=0;j<=n;j++){
                scanf("%d",&d[i][j]);
            }
        }
        for(int k=0;k<=n;k++)
            for(int i=0;i<=n;i++)
                for(int j=0;j<=n;j++)
                    d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
        memset(dp,0x3f,sizeof(dp));
        dp[1][0]=0;
        n++;
        for(int i=1;i<1<<n;i++)
            for(int j=0;j<n;j++) if(i>>j&1)
                for(int k=0;k<n;k++) if((i^1<<j)>>k&1)
                    dp[i][j]=min(dp[i][j],dp[i^1<<j][k]+d[k][j]);

        int ans=0x3f3f3f3f;
        for(int i=1;i<n;i++) ans=min(ans,dp[(1<<n)-1][i]+d[i][0]);
        printf("%d\n",ans); 

    }
}

HDU4856

題意

​ 有一個n*n的地圖,.表示可行走的路,#表示障礙,接下來有m行,x1,y2,x2,y2分別表示隧道的入口和出口,最後求走完這些隧道並且只能走一次需要多長時間。

思路

​ 首先bfs所有隧道從入口到另一個隧道出口的時間,如果不能到達將它置爲-1,然後就是初始化從每個起點出發(也就是for(int i=0;i<m;i++) dp[1<<i][i]=0;)剩下的就越是單純的TSP模板了。

代碼

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int INF = 0x3f3f3f3f;
struct point{
    int x1,y1,x2,y2;
}p[17];
int n,m;
char ma[17][17];
int mp[17][17],visit[17][17],dp[1<<15][17];
int step[4][2] = {{1, 0}, {-1, 0}, {0, 1}, {0, -1}};
bool check(int x,int y){
    return x<=n&&x>=1&&y<=n&&y>=1&&ma[x][y]!='#';
} 
int bfs(point x,point y){
    queue<pair<int,int> > q;
    memset(visit,-1,sizeof(visit));
    q.push(make_pair(x.x2,x.y2));
    visit[x.x2][x.y2]=0;
    while(q.size()){
        pair<int,int> now=q.front();
        q.pop();
        if(now.first==y.x1&&now.second==y.y1) return visit[y.x1][y.y1];
        for(int i=0;i<=3;i++){
            int tx=now.first+step[i][0],ty=now.second+step[i][1];
            if(check(tx,ty)&&visit[tx][ty]==-1) {
                visit[tx][ty]=visit[now.first][now.second]+1;
                q.push(make_pair(tx,ty));
            }
        }
    }
    return -1;
}
int main(){
    while(~scanf("%d%d",&n,&m)){
        for(int i=1;i<=n;i++) scanf("%s",ma[i]+1);
        for(int i=0;i<m;i++){
            scanf("%d%d%d%d",&p[i].x1,&p[i].y1,&p[i].x2,&p[i].y2);
        }
        for(int i=0;i<m;i++){
            for(int j=0;j<m;j++){
                if(i==j) mp[i][j]=0;
                else mp[i][j]=bfs(p[i],p[j]);
            }
        }
        for(int i=1;i<1<<m;i++){
            for(int j=0;j<m;j++){
                dp[i][j]=0x3f3f3f3f;
            }
        }
        for(int i=0;i<m;i++) dp[1<<i][i]=0;
        for(int i=0;i<1<<m;i++){
            for(int j=0;j<m;j++) if(i>>j&1)
                for(int k=0;k<m;k++) if(((i^1<<j)>>k&1)&&mp[k][j]!=-1)
                    dp[i][j]=min(dp[i][j],dp[i^1<<j][k]+mp[k][j]);
        }
        int ans=0x3f3f3f3f;
        for(int i=0;i<m;i++)  ans=min(ans,dp[(1<<m)-1][i]);
        if(ans==0x3f3f3f3f) printf("-1\n");
        else printf("%d\n",ans);
    }
    return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章