Brief Description

3333年,在银河系的某星球上,X军团和Y军团正在激烈地作战。在战斗的某一阶段,Y军团一共派遣了N个巨型机器人进攻X军团的阵地,其中第i个巨型机器人的装甲值为Ai。当一个巨型机器人的装甲值减少到0或者以下时,这个巨型机器人就被摧毁了。X军团有M个激光武器,其中第i个激光武器每秒可以削减一个巨型机器人Bi的装甲值。激光武器的攻击是连续的。这种激光武器非常奇怪,一个激光武器只能攻击一些特定的敌人。Y军团看到自己的巨型机器人被X军团一个一个消灭,他们急需下达更多的指令。为了这个目标,Y军团需要知道X军团最少需要用多长时间才能将Y军团的所有巨型机器人摧毁。但是他们不会计算这个问题,因此向你求助。

Algorithm Design

明明是非常水的套路题却卡了很久, 看来是太长时间没有做网络流的题目手生了, 我真的是思博啊, 我这样的人就应该早点滚粗吧Orz

我们考虑直接二分答案, 这样所有武器能输出的最大火力就有了, 按照火力建模, 随便搞搞就有了.

我真是思博啊.

数组开小了...然后不爽一下子开150万就过了Orz

Code

#include <algorithm>
#include <cmath>
#include <cstdio>
#include <queue>
const double eps = 1e-9;
#define ld long double
const int maxn = 110 << 1;
const int inf = 0x3f3f3f;
struct edge {
int to, next;
ld cap;
edge(int a = 0, int b = 0, ld c = 0) : to(a), next(b), cap(c) {}
} e[1500005];
int cnt = 1, n, m, s, t, v;
ld L, R;
int head[maxn], mi[maxn], A[maxn], B[maxn], tod;
ld dist[maxn];
void add_edge(int u, int v, ld c) {
e[++cnt] = edge(v, head[u], c);
head[u] = cnt;
}
void add(int u, int v, ld c) {
add_edge(u, v, c);
add_edge(v, u, 0);
}
void bfs(int s) {
std::fill(dist, dist + v, -1);
dist[s] = 0;
std::queue<int> q;
q.push(s);
while (!q.empty()) {
int u = q.front();
q.pop();
for (int i = head[u]; i; i = e[i].next) {
if (dist[e[i].to] == -1 && e[i].cap > eps) {
dist[e[i].to] = dist[u] + e[i].cap;
q.push(e[i].to);
}
}
}
}
ld dfs(int s, int t, ld flow) {
if (s == t)
return flow;
for (int i = head[s]; i; i = e[i].next) {
// printf("%.2f ", e[i].cap);
if (e[i].cap > eps && dist[e[i].to] > dist[s]) {
ld d = dfs(e[i].to, t, std::min(e[i].cap, flow));
if (d > eps) {
e[i].cap -= d;
e[i ^ 1].cap += d;
return d;
}
}
}
return 0;
}
ld dinic(int s, int t) {
ld flow = 0;
while (1) {
bfs(s);
if (dist[t] == -1)
return flow;
ld d;
while ((d = dfs(s, t, inf)))
flow += d;
}
return flow;
}
bool check(ld mid) {
for (int i = 2; i <= cnt; i += 2)
e[i].cap += e[i ^ 1].cap, e[i ^ 1].cap = 0;
for (int i = head[s]; i; i = e[i].next)
e[i].cap = mid * B[e[i].to];
ld flow = dinic(s, t);
return std::fabs(flow - (ld)tod) <= eps;
}
int main() {
#ifndef ONLINE_JUDGE
freopen("input", "r", stdin);
#endif
scanf("%d %d", &n, &m);
s = 0, t = n + m + 1, v = t + 1, L = 0, R = 0;
for (int i = 1; i <= n; i++) {
scanf("%d", &A[i]);
mi[i] = inf;
}
for (int i = 1; i <= m; i++)
scanf("%d", &B[i]);
for (int i = 1; i <= m; i++) {
for (int j = 1; j <= n; j++) {
int x;
scanf("%d", &x);
if (x == 1) {
add(i, j + m, inf);
mi[j] = std::min(mi[j], B[i]);
}
}
add(0, i, inf);
}
for (int i = 1; i <= n; i++) {
R += A[i] / mi[i];
add(i + m, t, A[i]);
tod += A[i];
}
// printf("%d\n", tod);
R = 5e6;
while (R - L > eps) {
ld mid = (L + R) / 2;
if (check(mid))
R = mid;
else
L = mid;
}
printf("%.7Lf", (L + R) / 2);
}

[bzoj3993][SDOI2015]星际战争-二分+最大流的更多相关文章

  1. 【BZOJ3993】[SDOI2015]星际战争 二分+最大流

    [BZOJ3993][SDOI2015]星际战争 Description 3333年,在银河系的某星球上,X军团和Y军团正在激烈地作战.在战斗的某一阶段,Y军团一共派遣了N个巨型机器人进攻X军团的阵地 ...

  2. bzoj千题计划131:bzoj3993: [SDOI2015]星际战争

    http://www.lydsy.com/JudgeOnline/problem.php?id=3993 二分答案 源点向武器连 mid*攻击力的边 机器人向汇点连 防御力 的边 武器i能攻击机器人j ...

  3. BZOJ 3993: [SDOI2015]星际战争 [二分答案 二分图]

    3993: [SDOI2015]星际战争 题意:略 R1D2T1考了裸二分答案+二分图最大匹配... #include <iostream> #include <cstdio> ...

  4. bzoj3993: [SDOI2015]星际战争(网络流)

    3993: [SDOI2015]星际战争 题目:传送门 题解: 洛谷AC了,但是因为bzoj的spj有问题所以暂时没A 一道老题目了,二分时间然后网络流判断. 每次st-->武器连时间*攻击力 ...

  5. 【二分答案】【最大流】bzoj3993 [Sdoi2015]星际战争

    二分Time,S->炮[i]:Time*b[i] 炮[i]->机器人[i]:INF 机器人[i]->T:a[i]. 判断是否满流. #include<cstdio> #i ...

  6. [SDOI2015][bzoj3993] 星际战争 [二分+最大流]

    题面 传送门 思路 首先,有一个非常明显的模型: 将激光武器放到一边,机器人放到另一边,从每一个激光武器向它可以攻击的机器人连边,形成一个二分图 建立附加源点$ss$和附加汇点$tt$,$ss$连所有 ...

  7. bzoj3993: [SDOI2015]星际战争(二分+最大流)

    题目描述 3333年,在银河系的某星球上,X军团和Y军团正在激烈地作战. 在战斗的某一阶段,Y军团一共派遣了N个巨型机器人进攻X军团的阵地,其中第i个巨型机器人的装甲值为Ai.当一个巨型机器人的装甲值 ...

  8. BZOJ3993 [SDOI2015]星际战争 【二分 + 网络流】

    题目 3333年,在银河系的某星球上,X军团和Y军团正在激烈地作战.在战斗的某一阶段,Y军团一共派遣了N个巨型机器人进攻X军团的阵地,其中第i个巨型机器人的装甲值为Ai.当一个巨型机器人的装甲值减少到 ...

  9. BZOJ3993 [SDOI2015]星际战争

    二分答案...然后最大流验证是否可行... 没了,好水啊QAQ /************************************************************** Prob ...

随机推荐

  1. c# 调用c++ 使用指针传递的时候

    http://www.cnblogs.com/warensoft/archive/2011/12/09/warenosoft3d.html 上面这篇文章很好解释了. 简单记录一下: 1. 声明  注意 ...

  2. [C/C++] 友元函数和友元类

    A---友元函数: class Data{ public: ... friend int f(int &m);//友元函数 ... } 友元函数是可以直接访问类的私有成员的非成员函数.它是定义 ...

  3. GSL介绍【转】

    GSL(GNU Scientific Library)是一个C写成的用于科学计算的库,下面是一些相关的包 http://www.thebigdata.cn/JiShuBoKe/5612.html

  4. Codeforces Round#516 Div.1 翻车记

    A:开场懵逼.然后发现有人1min过,于是就sort了一下,于是就过了.正经证明的话,考虑回文串两端点一定是相同的,所以最多有Σcnti*(cnti+1)/2个,cnti为第i种字母出现次数.而sor ...

  5. BZOJ4573:[ZJOI2016]大森林——题解

    http://www.lydsy.com/JudgeOnline/problem.php?id=4573 https://www.luogu.org/problemnew/show/P3348#sub ...

  6. Codeforces Round #466 (Div. 2) E. Cashback

    Codeforces Round #466 (Div. 2) E. Cashback(dp + 贪心) 题意: 给一个长度为\(n\)的序列\(a_i\),给出一个整数\(c\) 定义序列中一段长度为 ...

  7. 51nod 1275 连续子段的差异(twopointer+单调队列)

    对于每一个i找到最近的j满足最大值-最小值>K,对答案的贡献为j-i,用单调队列维护最值即可 #include<iostream> #include<cstdlib> # ...

  8. Linux系统通过AWS命令行上传文件至S3

    打开你的AWS控制台: 在IAM中创建一个新用户(比如test),创建时它会自动创建一个用户安全凭证,是由“访问密钥ID”和“私有访问密钥”组成的,请记住它并下载该凭证,后面会用到它: 选择你刚创建的 ...

  9. 将shell返回的结果保存至数组

    如下,我需要将u1和u2提取出保存至数组,方便后续的调用 root@ubuntu:~# lxc list+------+---------+------------------------------ ...

  10. HDU 5650 异或

    so easy Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Sub ...