uva10857

题意

兔子希望在平面上 n 个点上放蛋,每个点最多放一个蛋,初始兔子在 (0, 0) 点,这里有无数个蛋,兔子可以回到这个点取蛋,兔子的速度为 \(v * 2^{-i}\)(i 为携带蛋的数量)。对于每个点 (x, y) ,在 \(720+x/2000\) 分钟后不能在该点放蛋了。问最多放几个蛋。

分析

状态压缩DP。

dp[i][j][k] 表示到 i 点时,所有点的状态为 j ,身上携带 k 个蛋的最短时间。

可以使用 BFS 进行记忆化搜索。但是直接这样搜 \(O(2^{17}*18^3)\) 复杂度太高。

发现兔子速度的计算公式,\(v * 2^{-i}\),当 i 很大时,速度会非常小。

考虑一种特殊情况,只有两个点 A(0, a), B(0, a + b) ,如果每次只携带一个蛋,那么分别在 A B 放 2 个蛋所需时间为 \(a * 2 + a + a * 2 + b * 2 = 5a + 2b\) ,一次性携带两个蛋所需时间为 \(4a + 2b\),携带两个蛋更优;

假如有三个点 A(0, a) B(0, a + b) C(0, a + b + c),一次性携带三个蛋所需时间为 \(8a + 4b + 2c\),而选择先带一个蛋放到 A,再带两个蛋,所需时间为 \(2 * a + a + 4 * a + 4 * b + 2 * c = 7a + 4b + 2c\),也就是说分开来带更优,对于更多的点也可以证明。

所以 k 的最大值为 2,也就是说从初始点最多只能带一个或两个蛋出来。

code

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 18;
const double INF = 0x3f3f3f3f;
const double eps = 1e-9;
int n;
double v, dis[MAXN][MAXN];
double dist(double x1, double y1, double x2, double y2) {
return sqrt((x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1));
}
int x[MAXN], y[MAXN];
double dp[MAXN][1 << MAXN][3];
double maxtime[MAXN];
struct P {
int pos, state, cnt;
};
void bfs() {
queue<P> q;
q.push(P{0, 0, 2});
q.push(P{0, 0, 1});
dp[0][0][0] = 0;
dp[0][0][1] = 0;
dp[0][0][2] = 0;
while(!q.empty()) {
P p = q.front(); q.pop();
for(int i = 1; i <= n; i++) {
if(!((p.state >> i) & 1)) {
P np;
np.state = p.state + (1 << i);
if(p.cnt) {
np.cnt = p.cnt - 1;
np.pos = i;
double ctime = dis[p.pos][i] / v * (1 << p.cnt);
if(dp[p.pos][p.state][p.cnt] + ctime < maxtime[i]) {
if(fabs(dp[np.pos][np.state][np.cnt] - INF) < eps) { // eps
q.push(P{np.pos, np.state, np.cnt});
}
dp[np.pos][np.state][np.cnt] = min(dp[np.pos][np.state][np.cnt], dp[p.pos][p.state][p.cnt] + ctime);
}
} else {
double ctime = dis[p.pos][0] / v; // 从当前点回到 0 点
double ctime1 = dis[0][i] / v * 2, ctime2 = dis[0][i] / v * 4; // 从 0 点拿 1 个或 2 个到当前点的时间
np.pos = i;
np.cnt = 0;
if(dp[p.pos][p.state][p.cnt] + ctime + ctime1 < maxtime[i]) {
if(fabs(dp[np.pos][np.state][np.cnt] - INF) < eps) {
q.push(P{np.pos, np.state, np.cnt});
}
dp[np.pos][np.state][np.cnt] = min(dp[np.pos][np.state][np.cnt], dp[p.pos][p.state][p.cnt] + ctime + ctime1);
}
np.cnt = 1;
if(dp[p.pos][p.state][p.cnt] + ctime + ctime2 < maxtime[i]) {
if(fabs(dp[np.pos][np.state][np.cnt] - INF) < eps) {
q.push(P{np.pos, np.state, np.cnt});
}
dp[np.pos][np.state][np.cnt] = min(dp[np.pos][np.state][np.cnt], dp[p.pos][p.state][p.cnt] + ctime + ctime2);
}
}
}
}
}
}
int bits1[1 << MAXN];
int main() {
while(cin >> n >> v && (n + v)) {
x[0] = 0; y[0] = 0;
for(int i = 0; i <= n; i++) {
for(int j = 0; j < 3; j++) {
for(int k = 0; k < (1 << (1 + n)); k++) {
dp[i][k][j] = INF;
}
}
}
for(int i = 1; i <= n; i++) {
cin >> x[i] >> y[i];
}
for(int i = 0; i <= n; i++) {
maxtime[i] = (720 + (1.0 * x[i] / 2000)) * 60;
for(int j = 0; j <= n; j++) {
dis[i][j] = dist(x[i], y[i], x[j], y[j]);
}
}
bfs();
int ans = 0;
for(int i = 0; i <= n; i++) {
for(int k = 0; k < 3; k++) {
for(int j = 0; j < (1 << (n + 1)); j++) {
bits1[j] = bits1[j >> 1] + (j & 1);
if(dp[i][j][k] < INF) {
ans = max(ans, bits1[j]);
}
}
}
}
cout << ans << endl;
}
return 0;
}

uva10857(状态压缩DP)的更多相关文章

  1. hoj2662 状态压缩dp

    Pieces Assignment My Tags   (Edit)   Source : zhouguyue   Time limit : 1 sec   Memory limit : 64 M S ...

  2. POJ 3254 Corn Fields(状态压缩DP)

    Corn Fields Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 4739   Accepted: 2506 Descr ...

  3. [知识点]状态压缩DP

    // 此博文为迁移而来,写于2015年7月15日,不代表本人现在的观点与看法.原始地址:http://blog.sina.com.cn/s/blog_6022c4720102w6jf.html 1.前 ...

  4. HDU-4529 郑厂长系列故事——N骑士问题 状态压缩DP

    题意:给定一个合法的八皇后棋盘,现在给定1-10个骑士,问这些骑士不能够相互攻击的拜访方式有多少种. 分析:一开始想着搜索写,发现该题和八皇后不同,八皇后每一行只能够摆放一个棋子,因此搜索收敛的很快, ...

  5. DP大作战—状态压缩dp

    题目描述 阿姆斯特朗回旋加速式阿姆斯特朗炮是一种非常厉害的武器,这种武器可以毁灭自身同行同列两个单位范围内的所有其他单位(其实就是十字型),听起来比红警里面的法国巨炮可是厉害多了.现在,零崎要在地图上 ...

  6. 状态压缩dp问题

    问题:Ignatius has just come back school from the 30th ACM/ICPC. Now he has a lot of homework to do. Ev ...

  7. BZOJ-1226 学校食堂Dining 状态压缩DP

    1226: [SDOI2009]学校食堂Dining Time Limit: 10 Sec Memory Limit: 259 MB Submit: 588 Solved: 360 [Submit][ ...

  8. Marriage Ceremonies(状态压缩dp)

     Marriage Ceremonies Time Limit:2000MS     Memory Limit:32768KB     64bit IO Format:%lld & %llu ...

  9. HDU 1074 (状态压缩DP)

    题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=1074 题目大意:有N个作业(N<=15),每个作业需耗时,有一个截止期限.超期多少天就要扣多少 ...

随机推荐

  1. CentOS搭建pptpd服务笔记

    pptpd.rpm 包下载.http://poptop.sourceforge.net/yum/stable/packages/ 参考资料:http://www.oschina.net/questio ...

  2. 关于在虚拟机上安装iOS所遇到的问题

    先配一张成功的图片: 要想实现上图所示的,就要安装 unlocker206, 安装完unlocker206 不用改里面的内容!!!看到**经验上说还要改内容,结果我试了很多次,根本不行,其实问题不在改 ...

  3. selenium初识(一)

    Selenium是一个开源的便携式的自动化软件测试工具,用于测试web应用程序.有能力在不同浏览器和操作系统运行.它是一套工具,帮助我们有效地给予web应用程序的自动化. Selenium分为以下几个 ...

  4. Python下安装protobuf

    1. 下载安装包 2. 解压缩 tar –xzvf protobuf-2.6.1.tar.gz 3. 安装protoc 在python中使用protocbuf需要Protocal Buffer 编译器 ...

  5. hexo 添加标签

    --- title: title #文章標題 date: 2016-06-01 23:47:44 #文章生成時間 categories: "Hexo教程" #文章分類目錄 可以省略 ...

  6. [codeforces] 97B Superset || 平面分治

    原题 给出一个平面的一些点,让你添加点,使得所有点对满足以下三个要求中的一个: 1.在一个水平面上 2.在一个竖直线上 3.以这两个点为对角的矩形内包含有其他点 输出一种可行解 因为只需要可行解,且只 ...

  7. ZCC loves cube(cube)

    题目描述 调戏完了狗,ZCC开始玩起了积木.ZCC的面前有一块n*n的棋盘,他要用这些1*1*1的积木在棋盘上搭出一个宏伟的建筑.积木有三种颜色,ZCC认为一个建筑要被称为宏伟的应该满足能从正面看到的 ...

  8. 【HDU 5961 传递】

    Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Submission ...

  9. PHP AES128加密解密

    <?php /** * Class AES */ class AES { public static function encrypt($input, $key) { $size = mcryp ...

  10. 在AppCode中的razor调用HtmlHelper方法和UrlHelper方法

    原文发布时间为:2011-05-17 -- 来源于本人的百度文章 [由搬家工具导入] 可以写一个帮助类,如下 using System.Web.WebPages;using System.Web.Mv ...