【BZOJ 3993】【SDOI 2015】星际战争
http://www.lydsy.com/JudgeOnline/problem.php?id=3993
调了好长时间啊
这道题设时间为time,那么对于m个武器从S向这m个点连容量为time*Bi的边,代表能造成的总伤害。
对于每个武器向每个能打到的机器人连容量为无穷的边。
对每个机器人向T连自己的装甲量Ai为容量的边。
又因为每个机器人的装甲必须被打完,所以设机器人到T的边的下界为Ai(上下界相等)。
从T向S连容量无穷大的边,构成无源汇上下界网络流,所以构造附加网络然后二分时间time就可以了。
注意附加网络有一些边是走不通的,所以把那些边删掉也可以
(删掉后和原网络没什么区别啊Σ( ° △ °|||)︴那我还讲这么偏QAQ)。
Dinic模板打残致命伤!!!以后改数组边界的时候注意要看程序中调用数组边界的常量是否也需要改。
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define forU(i, a, b) for (int i = (a), UP = (b); i <= UP; ++i)
#define forD(i, a, b) for (int i = (a), DOWN = (b); i >= DOWN; --i)
using namespace std;
inline double fabs(double x) {return x < 0 ? -x : x;}
const int N = 55;
const int M = N * N * 2;
struct node {
int nxt, to; double c;
node(int _nxt = 0, int _to = 0, double _c = 0) : nxt(_nxt), to(_to), c(_c) {}
} E[M];
int cnt, point[N << 1], cur[N << 1];
void ins(int u, int v, double c) {
E[++cnt] = node(point[u], v, c); point[u] = cnt;
E[++cnt] = node(point[v], u, 0); point[v] = cnt;
}
bool can[N][N];
int q[N << 1], n, m, S, T, d[N << 1];
double A[N], B[N];
bool BFS() {
memset(d, 0, sizeof(int) * (T + 1));
int head = 0, tail = 1, u, v; q[1] = S; d[S] = 1;
while (head != tail) {
u = q[++head];
for (int i = point[u]; i; i = E[i].nxt)
if (fabs(E[i].c) > 1e-12 && !d[v = E[i].to]) {
d[v] = d[u] + 1;
q[++tail] = v;
}
}
return d[T];
}
double DFS(int u, double a) {
if (fabs(a) < 1e-12 || u == T) return a;
double f, flow = 0;
for (int &i = cur[u]; i; i = E[i].nxt)
if (d[E[i].to] == d[u] + 1 && fabs(f = DFS(E[i].to, min(a, E[i].c))) > 1e-12) {
a -= f; flow += f; E[i].c -= f; E[i ^ 1].c += f;
if (fabs(a) < 1e-12) break;
}
return flow;
}
double Dinic() {
double flow = 0;
while (BFS()) {
memcpy(cur, point, sizeof(int) * (T + 1));
flow += DFS(S, 100000000000000.0);
}
return flow;
}
double tot = 0;
bool can2(double tim) {
memset(point, 0, sizeof(int) * (T + 1));
cnt = 1;
forU (i, 1, m) {
ins(S, i, tim * B[i]);
forU (j, 1, n)
if (can[i][j])
ins(i, j + m, 100000000000000.0);
}
forU (i, 1, n) ins(i + m, T, A[i]);
return fabs(Dinic() - tot) < 1e-5;
}
int main() {
scanf("%d%d", &n, &m); S = n + m + 1; T = S + 1;
forU (i, 1, n) scanf("%lf", A + i), tot += A[i];
forU (i, 1, m) scanf("%lf", B + i);
int num;
double left = 0, right = 0, mid;
forU (i, 1, m)
forU (j, 1, n) {
scanf("%d", &num);
can[i][j] = num;
}
forU (i, 1, n) {
double fornow = 0;
forU (j, 1, m)
fornow = max(fornow, B[j]);
right += ceil(A[i] / fornow);
}
while (right - left > 1e-5) {
mid = (left + right) / 2.0;
if (can2(mid)) right = mid;
else left = mid;
}
printf("%lf\n", left);
return 0;
}
【BZOJ 3993】【SDOI 2015】星际战争的更多相关文章
- BZOJ 3993 [SDOI 2015] 星际战争 解题报告
首先我们可以二分答案. 假设当前二分出来的答案是 $Ans$ ,那么我们考虑用网络流检验: 设武器为 $X$,第 $i$ 个武器的攻击力为 $B_i$: 设机器人为 $Y$,第 $i$ 个机器人的装甲 ...
- [SDOI 2015] 星际战争
[题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=3993 [算法] 首先发现问题具有单调性 , 不妨二分答案mid 考虑网络流 : 将源 ...
- [BZOJ 3992] [SDOI 2015] 序列统计(DP+原根+NTT)
[BZOJ 3992] [SDOI 2015] 序列统计(DP+原根+NTT) 题面 小C有一个集合S,里面的元素都是小于质数M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数列,数 ...
- BZOJ 3990 [SDOI 2015] 排序 解题报告
这个题哎呀...细节超级多... 首先,我猜了一个结论.如果有一种排序方案是可行的,假设这个方案是 $S$ . 那么我们把 $S$ 给任意重新排列之后,也必然可以构造出一组合法方案来. 于是我们就可以 ...
- BZOJ 3992 [SDOI 2015] 序列统计 解题报告
这个题最暴力的搞法就是这样的: 设 $Dp[i][j]$ 为前 $i$ 个数乘积为 $j$ 的方案数. 转移的话就不多说了哈... 当前复杂度 $O(nm^2)$ 注意到,$M$ 是个质数,就说明 $ ...
- [BZOJ 3992] [SDOI 2015] 序列统计
Description 传送门 Solution [一] 设 \(f[i][j]\) 表示前 \(i\) 个数的乘积在模 \(p\) 意义下等于 \(j\) 的方案数,有 \[ f[i][j]=\su ...
- BZOJ 3993: [SDOI2015]星际战争 [二分答案 二分图]
3993: [SDOI2015]星际战争 题意:略 R1D2T1考了裸二分答案+二分图最大匹配... #include <iostream> #include <cstdio> ...
- BZOJ 3993 [SDOI2015]星际战争 | 网络流 二分答案
链接 BZOJ 3993 题解 这道题挺棵的-- 二分答案t,然后源点向武器连t * b[i], 武器向能攻击的敌人连1, 敌人向汇点连a[i],如果最大流等于所有敌人的a[i]之和则可行. #inc ...
- bzoj千题计划131:bzoj3993: [SDOI2015]星际战争
http://www.lydsy.com/JudgeOnline/problem.php?id=3993 二分答案 源点向武器连 mid*攻击力的边 机器人向汇点连 防御力 的边 武器i能攻击机器人j ...
- 3993: [SDOI2015]星际战争
3993: [SDOI2015]星际战争 Time Limit: 10 Sec Memory Limit: 128 MBSec Special JudgeSubmit: 1244 Solved: ...
随机推荐
- linux vi 命令大全
进入vi的命令 vi filename :打开或新建文件,并将光标置于第一行首 vi +n filename :打开文件,并将光标置于第n行首 vi + filename :打开文件,并将光标置于最后 ...
- 课堂Java小程序(加减乘除与验证码)
一.编写一个程序,用户输入两个数,求出其加减乘除,并用消息框 显示计算结果. 1.设计思想:从键盘输入两个数字和运算符,然后计算.将输入的数字及运算符由字符型转换为整型,再用if判断输入的运算符,根据 ...
- SpringMVC启动过程详解(li)
通过对SpringMVC启动过程的深入研究,期望掌握Java Web容器启动过程:掌握SpringMVC启动过程:了解SpringMVC的配置文件如何配置,为什么要这样配置:掌握SpringMVC是如 ...
- Console.In.ReadToEnd() 控制台 输入完毕
输入完数据后 按回车(另起一行) ctrl+z enter .......百度了半天 没百度到..最后还是google 强大..解决了问题 ..
- AlloyTouch与three.js 3D模型交互
如你所见,上面的cube的旋转.加速.减速停止都是通过AlloyTouch去实现的. 演示 代码 <script src="asset/three.js"></s ...
- Android Weekly Notes Issue #222
Android Weekly Issue #222 September 11th, 2016 Android Weekly Issue #222 ARTICLES & TUTORIALS Fo ...
- IOS 杂笔-14(被人遗忘的owner)
*owner在开发中现在已经很少用了 有兴趣的童鞋可以看看* 我们遇到owner通常是在类似 [[[NSBundle mainBundle] loadNibNamed:@"Food" ...
- linux php编译安装
1.下载php安装包 http://cn2.php.net/get/php-5.5.36.tar.gz/from/this/mirror 默认情况下Nginx和PHP他俩之间是一点感觉没有的.Apac ...
- SE Springer小组《Spring音乐播放器》软件需求说明3
3 需求规定 3.1对功能的规定 基本功能与相关的输入输出如下表所示.歌曲播放.停止.暂停等功能调用MCI库,数据在MCI库下如何运作与用户的直观感受无关,就不具体列出. 输入 处理 输出 用户登录信 ...
- [project euler] program 4
上一次接触 project euler 还是2011年的事情,做了前三道题,后来被第四题卡住了,前面几题的代码也没有保留下来. 今天试着暴力破解了一下,代码如下: (我大概是第 172,719 个解出 ...