uva10857(状态压缩DP)
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)的更多相关文章
- hoj2662 状态压缩dp
Pieces Assignment My Tags (Edit) Source : zhouguyue Time limit : 1 sec Memory limit : 64 M S ...
- POJ 3254 Corn Fields(状态压缩DP)
Corn Fields Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 4739 Accepted: 2506 Descr ...
- [知识点]状态压缩DP
// 此博文为迁移而来,写于2015年7月15日,不代表本人现在的观点与看法.原始地址:http://blog.sina.com.cn/s/blog_6022c4720102w6jf.html 1.前 ...
- HDU-4529 郑厂长系列故事——N骑士问题 状态压缩DP
题意:给定一个合法的八皇后棋盘,现在给定1-10个骑士,问这些骑士不能够相互攻击的拜访方式有多少种. 分析:一开始想着搜索写,发现该题和八皇后不同,八皇后每一行只能够摆放一个棋子,因此搜索收敛的很快, ...
- DP大作战—状态压缩dp
题目描述 阿姆斯特朗回旋加速式阿姆斯特朗炮是一种非常厉害的武器,这种武器可以毁灭自身同行同列两个单位范围内的所有其他单位(其实就是十字型),听起来比红警里面的法国巨炮可是厉害多了.现在,零崎要在地图上 ...
- 状态压缩dp问题
问题:Ignatius has just come back school from the 30th ACM/ICPC. Now he has a lot of homework to do. Ev ...
- BZOJ-1226 学校食堂Dining 状态压缩DP
1226: [SDOI2009]学校食堂Dining Time Limit: 10 Sec Memory Limit: 259 MB Submit: 588 Solved: 360 [Submit][ ...
- Marriage Ceremonies(状态压缩dp)
Marriage Ceremonies Time Limit:2000MS Memory Limit:32768KB 64bit IO Format:%lld & %llu ...
- HDU 1074 (状态压缩DP)
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=1074 题目大意:有N个作业(N<=15),每个作业需耗时,有一个截止期限.超期多少天就要扣多少 ...
随机推荐
- 【Spiral Matrix II】cpp
题目: Given an integer n, generate a square matrix filled with elements from 1 to n2 in spiral order. ...
- [网站公告]又拍云API故障造成图片无法上传
大家好,今天早上8:30左右发现又拍云API出现故障,造成图片无法上传,调用图片上传API时出现错误:“The operation has timed out”. 该故障给大家带来了麻烦,望大家谅解! ...
- Java EE - JSP 小结
Table of Contents 前言 JSP 与 Servlet JSP 初始化参数 脚本元素 page 指令 禁用脚本元素 EL 表达式 EL 函数 taglib 指令 标记 TLD 文件的位置 ...
- ASP.NET Core [3]:进入HttpContext的世界(笔记)
原文链接:http://www.cnblogs.com/RainingNight/p/httpcontext-in-asp-net-core.html HttpContext是ASP.NET中的核心对 ...
- 贪吃蛇StringBuilder 和 定时器
ConsoleKeyInfo info = Console.ReadKey(); while (true) { if (info.Key == ConsoleKey.UpArrow)//只能输入一次但 ...
- jquery select chosen 动态绑定值
$("#ddlMstData").find("option[value=" + data.MstKey + "]").attr(" ...
- crond守护进程
Linux系统任务计划/etc/crontab cron的主配置文件,可以定义PATHcron格式如下:# .----------------分钟 (0 - 59)# | .------------- ...
- easyui中datagrid空数据集不刷新的解决方式
datagrid空间可以异步请求json数据,并将新数据覆盖原有数据,重绘数据表. 但是当回来空数据集的时候,js会产生这样一条错误: TypeError: rows is null for(var ...
- easyui中的依赖关系
参考自:http://www.easyui.info/archives/765.html 在使用easyui的过程中发现各个组件直接存在依赖关系,也就是上层的复杂组件依赖于一个或者多个简单组件,复杂组 ...
- BZOJ5297 [Cqoi2018]社交网络 【矩阵树定理】
题目链接 BZOJ5297 题解 最近这玩意这么那么火 这题要用到有向图的矩阵树定理 主对角线上对应入度 剩余位置如果有边则为\(-1\),不然为\(0\) \(M_{i,i}\)即为以\(i\)为根 ...