hdu 3001 Travelling TSP变形 三进制状压dp

// hdu 3001 TSP问题的变形
// 这次到每个点最多两次,所以可以用三进制的类推
// dp[S][u]表示当前在u点访问状态为S时所得到的最小的开销
// 采用刷表法,即用当前的状态推出它所能转移的状态
// dp[S][u] 可以到达的状态为dp[S+state[v]][v](dist[u][v]!=inf)
// dp[S+state[v]][v] = max(dp[S+state[v]][v],dp[S][u]+dist[u][v]);
// 其中每个点最多访问2次
// 技巧就是先把所有的状态表示出来visit[S][u];
// 表示当前在u点访问状态为S时u的访问状态。
// 其中我们要求的就是最后每个点都访问过的状态S时的dp[S][v]{v<n};
//
// 这道题开始完全想套一个裸的TSP模板过的,结果样例都过不了
// 哎,发现自己一直很水,继续加油吧。。。
// 而且,犯了一个非常严重的错误,导致调了2个小时的bug
// 哎。。。加油,不过最终还是ac啦
//const int maxn = 15;
//int dp[1<<maxn][maxn];
//int mp[maxn][maxn];
//int n,m;
//void init(){
//	int a,b,c;
//	memset(mp,0,sizeof(mp));
//	memset(dp,0x3f,sizeof(dp));
//	for (int i=0;i<m;i++){
//		scanf("%d%d%d",&a,&b,&c);
//		mp[a][b] = mp[b][a] = c;
//	}
//	for (int i=0;i<(1<<n);i++)
//		dp[i][0] = mp[0][i];
//}
//
//void print(){
//	for (int i=0;i<(1<<n);i++){
//		for (int j=0;j<n;j++)
//			printf("%d ",dp[i][j]);
//		puts("");
//	}
//}
//
//void solve(){
//	for (int S=0;S<(1<<n);S++){
//		for (int u=1;u<=n;u++)
//			for (int v=1;v<=n;v++){
//				if ((S>>v)&1){
//					dp[S][u] = min(dp[S][u],dp[S^(1<<v)][v]+mp[v][u]);
//				}
//			}
//	}
//	print();
//	int ans = 0;
//	for (int i=0;i<n;i++)
//		ans = max(ans,dp[(1<<n)-1][i]);
//	printf("%d\n",ans);
//}
//
//int main() {
//    freopen("G:\\Code\\1.txt","r",stdin);
//	while(scanf("%d%d",&n,&m)!=EOF){
//		init();
//		solve();
//	}
//	return 0;
//}

#include <algorithm>
#include <bitset>
#include <cassert>
#include <cctype>
#include <cfloat>
#include <climits>
#include <cmath>
#include <complex>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <deque>
#include <functional>
#include <iostream>
#include <list>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <stack>
#include <vector>
#define ceil(a,b) (((a)+(b)-1)/(b))
#define endl '\n'
#define gcd __gcd
#define highBit(x) (1ULL<<(63-__builtin_clzll(x)))
#define popCount __builtin_popcountll
typedef long long ll;
using namespace std;
const int MOD = 1000000007;
const long double PI = acos(-1.L);

template<class T> inline T lcm(const T& a, const T& b) { return a/gcd(a, b)*b; }
template<class T> inline T lowBit(const T& x) { return x&-x; }
template<class T> inline T maximize(T& a, const T& b) { return a=a<b?b:a; }
template<class T> inline T minimize(T& a, const T& b) { return a=a<b?a:b; }

const int maxn = 12;
const int maxm = 60000;
int visit[maxm][maxn];
int dp[maxm][maxn];
int state[maxn];
int mp[maxn][maxn];
const int inf = 0x3f3f3f3f;
int n,m;

void print(){
	for (int i=0;i<n;i++){
		for (int j=0;j<n;j++)
			printf("%d ",mp[i][j]);
		puts("");
	}
}

void init(){
	memset(mp,inf,sizeof(mp));
	memset(dp,inf,sizeof(dp));
	for (int i=0;i<n;i++){
		dp[state[i]][i] = 0;
	}
	int a,b,c;
	for (int i=0;i<m;i++){
		scanf("%d%d%d",&a,&b,&c);
		a--,b--;
		mp[a][b] = mp[b][a] = min(mp[a][b],c);
	}
	//for (int i=0;i<n;i++)
	//	mp[i][i] = 0;
//	print();
}
void solve(){
	int ans = inf;
	for (int S=0;S<state[n];S++){
		bool flag = false;
		for (int u=0;u<n;u++){
			if (visit[S][u] == 0)	flag = true;
			if (dp[S][u] == inf)	continue;
			for (int v=0;v<n;v++) if (v!=u){
				if (visit[S][v]>=2)	continue;
				if (mp[u][v] == inf) continue;
				dp[S+state[v]][v] = min(dp[S+state[v]][v],dp[S][u]+mp[u][v]);
			}
		}
		if (!flag){
			for (int u=0;u<n;u++)
				ans = min(ans,dp[S][u]);
		}
	}
	if (ans == inf)
		ans = -1;
	printf("%d\n",ans);

}

void get(){
	state[0]=1;
	for (int i=1;i<=10;i++)
		state[i] = state[i-1] * 3;
	for (int S=0;S<=state[10];S++){
		int x = S;
		for (int i=0;i<=10;i++){
			visit[S][i] = x%3;
			x/=3;
		}
	}
}

int main() {
	//freopen("G:\\Code\\1.txt","r",stdin);
	get();
	while(scanf("%d%d",&n,&m)!=EOF){
		init();
		solve();
	}
	return 0;
}

时间: 2024-10-13 14:52:26

hdu 3001 Travelling TSP变形 三进制状压dp的相关文章

HDU 3001 三进制状压DP

N个城市,M条道路,每条道路有其经过的代价,每一个城市最多能够到达两次,求走全然部城市最小代价,起点随意. 三进制状压.存储每一个状态下每一个城市经过的次数. 转移方程: dp[i+b[k]][k]=Min(dp[i+b[k]][k],dp[i][j]+dis[j][k]); #include "stdio.h" #include "string.h" const int inf=0x3f3f3f3f; int b[15],mark[60010][15],dp[60

hdu 3001 Travelling 经过所有点(最多两次)的最短路径 三进制状压dp

题目链接 题意 给定一个\(N\)个点的无向图,求从任意一个点出发,经过所有点的最短路径长度(每个点至多可以经过两次). 思路 状态表示.转移及大体思路 与 poj 3311 Hie with the Pie 经过所有点(可重)的最短路径 floyd + 状压dp 相同. 但,因为是每个点 至多可以经过两次,所以应该用 三进制 来表示状态. 因为三进制不能直接通过移位来表示,所以要 预处理 出每个数字\(state\)的三进制表示中每一位\(i\)上的值\(dig[state][i]\). 注意

HDU 3001 三进制 状压dp

Travelling Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 3789    Accepted Submission(s): 1182 Problem Description After coding so many days,Mr Acmer wants to have a good rest.So travelling is

hdu 3001 Travelling(状态压缩 三进制)

Travelling Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 6600    Accepted Submission(s): 2144 Problem Description After coding so many days,Mr Acmer wants to have a good rest.So travelling is

Travelling (三进制+状压dp)

题目链接 1 #include <bits/stdc++.h> 2 using namespace std; 3 typedef long long ll; 4 inline ll read(){ 5 int x=0,f=1;char ch=getchar(); 6 while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();} 7 while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch

三进制状压 HDOJ 3001 Travelling

题目传送门 题意:从某个点出发,所有点都走过且最多走两次,问最小花费 分析:数据量这么小应该是状压题,旅行商TSP的变形.dp[st][i]表示状态st,在i点时的最小花费,用三进制状压.以后任意进制状压都会了. #include <bits/stdc++.h> using namespace std; const int INF = 0x3f3f3f3f; const int N = 11; int dp[60000][10]; int bit[N]; int d[N][N]; int p[

hdu 3341 Lost&#39;s revenge(AC自动机+变进制状压DP)

题目链接:hdu 3341 Lost's revenge 题目大意:给定一些需要匹配的串,然后在给定一个目标串,现在可以通过交换目标串中任意两个位置的字符,要求最 后生成的串匹配尽量多的匹配串,可以重复匹配. 解题思路:这题很明显是AC自动机+DP,但是dp的状态需要开40?40?40?40(记录每种字符的个数),空间承受 不了,但是其实因为目标串的长度有限,为40:所以状态更本不需要那么多,最多只有10?10?10?10,但是通过 40进制的hash转换肯定是不行,可以根据目标串中4种字符的个

三进制状压(涂抹果酱)

题意:https://www.acwing.com/problem/content/1067/ Tyvj 两周年庆典要到了,Sam 想为 Tyvj 做一个大蛋糕. 蛋糕俯视图是一个 N×M 的矩形,它被划分成 N×M 个边长为 1×1×1 的小正方形区域(可以把蛋糕当成 N 行 M 列的矩阵). 蛋糕很快做好了,但光秃秃的蛋糕肯定不好看! 所以,Sam 要在蛋糕的上表面涂抹果酱. 果酱有三种,分别是红果酱.绿果酱.蓝果酱,三种果酱的编号分别为 1,2,3. 为了保证蛋糕的视觉效果,Admin 下

hdu 5135 Little Zu Chongzhi&#39;s Triangles 状压DP

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5135 从大到小贪心是不对的! 比如 4 30 40 50 89 这组数据 结果应该是600 如果取后三条边只会是一个细长细长的三角形根本没多大面积 只不过因为数据水然后贪心碰巧能过而已 看到网上题解很多人用贪心我好郁闷... 状压DP 看到边的数目很少就应该想到这个了 先枚举可能出现的所有三角形 记录他们使用了那些边还有他们的面积 然后再枚举状态之中 套一层遍历所有三角形的循环 如果当前状态和三角形