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),每个作业需耗时,有一个截止期限.超期多少天就要扣多少 ...
随机推荐
- USACO Section1.1 Broken Necklace 解题报告
beads解题报告 —— icedream61 博客园(转载请注明出处)---------------------------------------------------------------- ...
- php静态文件返回304
有时一些静态文件(如图片)会由php输出,会发现请求都是200,静态文件每次都去服务器上请求太浪费资源了,这时如何让浏览器缓存图片呢?就需要我们在php中输出304了. 我们可以利用php中的 HTT ...
- Python 字符串格式化输出方式
字符串格式化有两种方式:百分号方式.format方式. 其中,百分号方式比较老,而format方式是比较先进的,企图替代古老的方式,目前两者共存. 1.百分号方式 格式:%[(name)][flags ...
- 查看Linux系统用户登录日志
日志: /var/log/wtmp说明: 此文件是二进制文件,查看方法如下:该日志文件永久记录每个用户登录.注销及系统的启动.停机的事件.因此随着系统正常运行时间的增加,该文件的大小也会越来越大,增加 ...
- 一个初学者的辛酸路程-jQuery
前言: 主要概要: 1.HTML+CSS补充 2.DOM事件 3.jQuery示例 内容概要: 1.布局 代码如下 <!DOCTYPE html> <html lang=" ...
- 【转载】全面解析Unity3D自动生成的脚本工程文件
我们在Unity3D开发的时候,经常会看到它会产生不少固定命名工程文件,诸如: Assembly-CSharp-vs.csproj Assembly-CSharp-firstpass-vs.csp ...
- Linq中Count()和Any()引发的效率问题
1.count和any 今天看了鹤冲天的文章:Linq:切勿使用 Count() > 0 来判断集合非空 有所收获,写下文章总结一下: 先看如下代码: static void Main(st ...
- ASP.NET Core [2]:Middleware-请求管道的构成(笔记)
原文链接:http://www.cnblogs.com/RainingNight/p/middleware-in-asp-net-core.html 中间件处理请求主要分为三个阶段:1. 中间件的注册 ...
- 构建Docker镜像两种方式的比较-Dockerfile方式和S2I方式
前言 写Dockerfile是构建Docker镜像最通常的方式,接触过Docker的童鞋多少了解一些.前段时间研究OpenShift(paas的一种),发现了另外一种构建Docker镜像的方式:S2I ...
- Linux查看端口被占用情形
查看某端口的占用情况: lsof -i:<端口号> 例如:lsof -i:8080 netstat -apn|grep <端口号> 例如: netstat -apn | gre ...