首页 > 代码库 > poj 3311 状压DP
poj 3311 状压DP
经典TSP变形
学到:1、floyd O(n^3)处理任意两点的最短路
2、集合的位表示,我会在最后的总结出写出。注意写代码之前一定设计好位的状态,本题中,第0位到第n位分别代表第i个城市,1是已经走过,0没走过
那么DP方程 :dp[s][i]--当前在城市i,状态为s(s存储的是走过了那些城市)
3、最后要求形成回路,那么就是min(dp[1<<(n+1)-1][i],dp[0][i])
#include <cstdio> #include <cstring> #include <algorithm> #include <string> #include <iostream> #include <cmath> #include <map> #include <queue> using namespace std; #define ls(rt) rt*2 #define rs(rt) rt*2+1 #define ll long long #define rep(i,s,e) for(int i=s;i<e;i++) #define repe(i,s,e) for(int i=s;i<=e;i++) #define CL(a,b) memset(a,b,sizeof(a)) #define IN(s) freopen(s,"r",stdin) const int MAXN = 12; int dis[MAXN][MAXN]; int dp[1<<MAXN][MAXN]; const int INF = 1e9+10; int n; void floyd() { rep(k,0,n+1) rep(i,0,n+1) rep(j,0,n+1) dis[i][j]=min(dis[i][k]+dis[k][j],dis[i][j]); } int main() { //IN("poj3311.txt"); int len; while(~scanf("%d",&n) && n) { rep(i,0,n+1) rep(j,0,n+1) dp[i][j]=dis[i][j]=INF; rep(i,0,n+1) rep(j,0,n+1) { scanf("%d",&len); dis[i][j]=min(dis[i][j],len); } floyd();//求出任意两点的距离 int S=1<<(n+1); rep(i,0,S) rep(j,0,n+1) { dp[i][j]=INF; } for(int s=0;s<S;s++)//枚举所有的状态 rep(i,0,n+1) { if(s&(1<<(i))) { if(s==(1<<i))dp[s][i]=dis[0][i]; else rep(j,0,n+1) if(s&(1<<j) && i!=j) { dp[s][i]=min(dp[s^(1<<i)][j]+dis[j][i],dp[s][i]); } } } int ans=INF; for(int i=0;i<n+1;i++) ans=min(ans,dp[(S-1)][i]+dis[i][0]); printf("%d\n",ans); } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。