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】星际战争的更多相关文章

  1. BZOJ 3993 [SDOI 2015] 星际战争 解题报告

    首先我们可以二分答案. 假设当前二分出来的答案是 $Ans$ ,那么我们考虑用网络流检验: 设武器为 $X$,第 $i$ 个武器的攻击力为 $B_i$: 设机器人为 $Y$,第 $i$ 个机器人的装甲 ...

  2. [SDOI 2015] 星际战争

    [题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=3993 [算法] 首先发现问题具有单调性 , 不妨二分答案mid 考虑网络流 : 将源 ...

  3. [BZOJ 3992] [SDOI 2015] 序列统计(DP+原根+NTT)

    [BZOJ 3992] [SDOI 2015] 序列统计(DP+原根+NTT) 题面 小C有一个集合S,里面的元素都是小于质数M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数列,数 ...

  4. BZOJ 3990 [SDOI 2015] 排序 解题报告

    这个题哎呀...细节超级多... 首先,我猜了一个结论.如果有一种排序方案是可行的,假设这个方案是 $S$ . 那么我们把 $S$ 给任意重新排列之后,也必然可以构造出一组合法方案来. 于是我们就可以 ...

  5. BZOJ 3992 [SDOI 2015] 序列统计 解题报告

    这个题最暴力的搞法就是这样的: 设 $Dp[i][j]$ 为前 $i$ 个数乘积为 $j$ 的方案数. 转移的话就不多说了哈... 当前复杂度 $O(nm^2)$ 注意到,$M$ 是个质数,就说明 $ ...

  6. [BZOJ 3992] [SDOI 2015] 序列统计

    Description 传送门 Solution [一] 设 \(f[i][j]\) 表示前 \(i\) 个数的乘积在模 \(p\) 意义下等于 \(j\) 的方案数,有 \[ f[i][j]=\su ...

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

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

  8. BZOJ 3993 [SDOI2015]星际战争 | 网络流 二分答案

    链接 BZOJ 3993 题解 这道题挺棵的-- 二分答案t,然后源点向武器连t * b[i], 武器向能攻击的敌人连1, 敌人向汇点连a[i],如果最大流等于所有敌人的a[i]之和则可行. #inc ...

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

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

  10. 3993: [SDOI2015]星际战争

    3993: [SDOI2015]星际战争 Time Limit: 10 Sec  Memory Limit: 128 MBSec  Special JudgeSubmit: 1244  Solved: ...

随机推荐

  1. Redis命令拾遗二(散列类型)

    本文版权归博客园和作者吴双共同所有,欢迎转载,转载和爬虫请注明原文地址 :博客园蜗牛NoSql系列地址  http://www.cnblogs.com/tdws/tag/NoSql/ Redis命令拾 ...

  2. alienware Win8 系统安装

    原作者网名 alienware-小来: 我的外星人 老是装系统出错.我觉得写的不错.把原作者的东西拿过来.. 对于win7系统的用户来说想要体验下win8.1系统,或者是原来win8.1系统加装固态后 ...

  3. servlet开发中遇到的问题集合

    问题1: servlet插入数据库时中文会乱码. 解决方法:在数据库连接地址最后增加两个转码参数(?useUnicode=true&characterEncoding=utf8) url=jd ...

  4. Java多线程--让主线程等待子线程执行完毕

    使用Java多线程编程时经常遇到主线程需要等待子线程执行完成以后才能继续执行,那么接下来介绍一种简单的方式使主线程等待. java.util.concurrent.CountDownLatch 使用c ...

  5. BroadcastReceiver几种常见监听

    1.BroadcastReceiver监听拨号 <intent-filter android:priority="1000" > <action android: ...

  6. 【blade利刃出鞘】一起进入移动端webapp开发吧

    前言 在移动浪潮袭来的时候,小钗有幸进入框架组做webapp框架开发,过程中遇到了移动端的各种坑,也产生了各种激情,就我们公司的发展历程来说 第一阶段:使用传统方式开发移动站点,少量引入HTML5元素 ...

  7. 自定义PopupWindow

    PopupWindow,一个弹出窗口控件,可以用来显示任意View,而且会浮动在当前activity的顶部 自定义PopupWindow. 1.extends PopupWindow 2.构造方法中可 ...

  8. css权威指南--笔记

    第1章 css和文档 1,元素:替换元素(img input),非替换元素(大多数span). 2,link:rel(代表关系:stylesheet,候选样式表:alternate styleshee ...

  9. gulp入门小记

    由于我所在的项目组一直在用gulp构建工具,而我只是在前人搭好的环境下每次运行gulp packJs来打包js,对里面的东西全然不知,刚好最近有些时间就想自己从学学将gulp怎么用于构建前端项目中,这 ...

  10. 01背包问题python 2.7实现

    版权声明:本文为博主原创文章,转载请注明转自 http://www.cnblogs.com/kdxb/p/6140625.html #!/usr/bin/env python # -*- coding ...