状态压缩dp,经典MTSP问题,第一问比较容易求,先求出满足情况的状态放在了mov里面,然后模拟0-1背包的过程求出最少用多少judge。第二问就比较复杂了,第二问要求我们算出最短的时间,即最短总路程。这样的话就比较难求了。。。第二问看了一下网上大神的代码才懂的。。。汗。。。第二问,dt[i][j],表示当状态为i(注意:i & 1 == 1)时, 最后一个走过的节点是j时最短路径, best[i] 表示状态完成状态 i 需要的最短路径。然后具体实现还是看代码吧。。。
#include<algorithm>
#include<iostream>
#include<cstring>
#include<vector>
#include<cstdio>
#include<cmath>
#define SL strlen
#define PB push_back
#define LL long long
#define INF 0X3f3f3f3f
#define CLR(a, b) memset(a, b, sizeof(a))
using namespace std;
const int N = 17;
struct Point
{
int x, y;
}p[N];
int dp[1 << N], c[N];
int dt[1 << N][N], best[1 << N];
int g[N][N];
vector<int> mov;
int n, m, i, j, k;
int len(Point a, Point b)
{
return ceil(sqrt((double)(a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y)));
}
bool ok(int k)
{
int ret = 0; j = 0;
while(k)
{
if(k & 1)
{
ret += c[j];
}
k >>= 1;
j ++;
}
return ret <= m;
}
void input()
{
for(i = 0; i < n; i ++)
{
scanf("%d%d", &p[i].x, &p[i].y);
}
for(i = 0; i < n; i ++)
{
scanf("%d", &c[i]);
}
}
void work()
{
for(i = 0; i < n; i ++)
{
for(j = 0; j < n; j ++)
{
g[i][j] = len(p[i], p[j]);
}
}
mov.clear();
for(i = 1; i < (1 << n); i ++)
{
if(ok(i)) mov.PB(i);
}
}
bool solve()
{
CLR(dp, INF);
dp[0] = 0;
for(i = 0; i < mov.size(); i ++)
{
for(j = (1 << n) - 1; j >= 1; j --)
{
if((j & mov[i]) == mov[i])
{
dp[j] = min(dp[j], dp[j ^ mov[i]] + 1);
}
}
}
if(dp[(1 << n) - 1] == INF) { printf("-1 ");return 0; }
//else printf("%d ", dp[(1 << n) - 1]); return 1;
}
void TSP()
{
CLR(best, INF);
CLR(dt, INF);
dt[1][0] = 0;
for(i = 0; i < mov.size(); i ++)
{
for(j = 0; j < n; j ++)
{
if((1 << j) & mov[i])
{
best[mov[i]] = min(best[mov[i]], dt[mov[i]][j] + g[j][0]);
for(k = 0; k < n; k ++)
{
if(!((1 << k) & mov[i]))
{
dt[mov[i] | (1 << k)][k] = min(dt[mov[i] | (1 << k)][k], dt[mov[i]][j] + g[j][k]);
}
}
}
}
}
for(i = 1; i < (1 << n); i ++)
{
if(i & 1) for(j = i & (i - 1); j; j = i & (j - 1))
{
best[i] = min(best[i], best[j] + best[(i ^ j) | 1]);
}
}
printf("%d\n", best[(1 << n) - 1]);
}
int main()
{
//freopen("input.txt", "r", stdin);
while(scanf("%d%d", &n, &m) != EOF)
{
input();
work();
if(solve()) TSP();
else puts("-1");
}
}
分享到:
相关推荐
HDU的1250,主要是利用高精度加法,但是代码有点繁琐,效率不是很高
HDU1059的代码
杭电ACMhdu1163
hdu1001解题报告
hdu 1574 passed sorce
HDU的一题........HDU DP动态规
hdu acm 教案 搜索入门 hdu acm 教案 搜索入门
hdu2101AC代码
搜索 dfs 解题代码 hdu1241
hdu acm 教案 动态规划(1) hdu acm 教案 动态规划(1)
hdu 5007 Post Robot 字符串枚举。 暴力一下就可以了。
ACM HDU题目分类,我自己总结的大概只有十来个吧
hdu1290 解题报告 献给杭电五十周年校庆的礼物 (切西瓜问题,即平面分割空间)
HDU最全ac代码
hdu 1166线段树代码
hdu动态规划算法集锦
hdu题目分类
自己做的HDU ACM已经AC的题目
hdu-acm源代码(上百题)hdu-acm源代码、hdu-acm源代码hdu-acm源代码
HDU图论题目分类