[题目链接]

https://www.lydsy.com/JudgeOnline/problem.php?id=3993

[算法]

首先发现问题具有单调性 , 不妨二分答案mid

考虑网络流 :

将源点向每个"激光武器”连一条流量为mid * Bi的边

将每个“激光武器”向每个其可以攻击的“机器人”连一条流量为正无穷的边

将每个“机器人”向汇点连一条流量为Ai的边

判断是否满流即可

时间复杂度 : O(dinic(N + M , M ^ 2) * logV)

[代码]

为避免精度误差 , 可以在整数域上进行二分 , 最后以浮点数形式输出

#include<bits/stdc++.h>
using namespace std;
#define N 510
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
const ll inf = 1e15;struct edge
{int to;ll w;int nxt;
} e[N * N * 5];int n , m , tot , S , T;
int dep[N] , head[N] , g[N][N];
ll a[N] , b[N];template <typename T> inline void chkmax(T &x,T y) { x = max(x,y); }
template <typename T> inline void chkmin(T &x,T y) { x = min(x,y); }
template <typename T> inline void read(T &x)
{T f = 1; x = 0;char c = getchar();for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;for (; isdigit(c); c = getchar()) x = (x << 3) + (x << 1) + c - '0';x *= f;
}
inline void addedge(int u , int v , ll w)
{++tot;e[tot] = (edge){v , w , head[u]};head[u] = tot;++tot;e[tot] = (edge){u , 0 , head[v]};head[v] = tot;
}
inline bool bfs()
{queue< int > q;for (int i = 1; i <= T; ++i)dep[i] = -1;q.push(S);dep[S] = 1;while (!q.empty()){int cur = q.front();q.pop();for (int i = head[cur]; i; i = e[i].nxt){int v = e[i].to;ll w = e[i].w;if (w > 0 && dep[v] == -1){dep[v] = dep[cur] + 1;q.push(v);if (v == T) return true;}}}return false;
}
inline ll dinic(int u , ll flow)
{ll k , rest = flow;if (u == T)return flow;for (int i = head[u]; i && rest; i = e[i].nxt){int v = e[i].to;ll w = e[i].w;if (dep[v] == dep[u] + 1 && w){k = dinic(v , min(w , rest));e[i].w -= k;e[i ^ 1].w += k;if (!k) dep[v] = 0;rest -= k;}}return flow - rest;
}
inline bool check(ll mid)
{S = n + m + 1 , T = S + 1;for (int i = 1; i <= T; ++i) head[i] = 0;for (int i = 1; i <= tot; ++i) e[i].nxt = 0;tot = 1;for (int i = 1; i <= m; ++i) addedge(S , i , b[i] * mid);for (int i = 1; i <= m; ++i){for (int j = 1; j <= n; ++j){if (g[i][j])addedge(i , j + m , inf);}}ll sum = 0;for (int i = 1; i <= n; ++i) {addedge(i + m , T , a[i]);sum += a[i];}ll res = 0;while (bfs()){while (double flow = dinic(S , inf)) res += flow;}return res == sum;
}int main()
{scanf("%d%d" , &n , &m); for (int i = 1; i <= n; ++i) {scanf("%lld" , &a[i]);a[i] *= 1000;}for (int i = 1; i <= m; ++i) scanf("%lld" , &b[i]);for (int i = 1; i <= m; ++i){for (int j = 1; j <= n; ++j){scanf("%d" , &g[i][j]);}}ll l = 0 , r = inf , ans = 0;while (l <= r){int mid = (l + r) >> 1;if (check(mid)){ans = mid;r = mid - 1;} else l = mid + 1;}printf("%.6lf\n" , (double)(ans / 1000.0));return 0;}

转载于:https://www.cnblogs.com/evenbao/p/10549250.html

[SDOI 2015] 星际战争相关推荐

  1. BZOJ3993: [SDOI2015]星际战争

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

  2. 星际战争2服务器维护,星际战争2初次反击虚空技能洗练技巧

    <星际战争2初次反击>关于虚空的重要性不必多说了,这里主要讲一下虚空如何洗练技能比较好,下面一起来看看吧. 虚空在pvp站场上的重要性不亚于顶级的武器.好的虚空是你吊打高玩的必备之物.首先 ...

  3. LRU算法,走迷宫,数根,星际战争

    问答题1 问答题1:一进程刚获得三个主存块的使用权,若该进程访问页面的次序是{1321215123},采用 LRU 算法,缺页数是几次? LRU 算法是先置换出不常用的页面,如果一个页面刚刚被访问过, ...

  4. XMU2018摸底测试 星际战争I(Easy)

    星际战争I(Easy) Description(描述) 现在是星际战争的第一个年头. 绝地武士躲藏在一颗大小为N \times MN×M 的矩形二维星球上. 贸易联盟想要消灭他们,但是经历过之前的战斗 ...

  5. 【Luogu-P3324 [SDOI2015] / DSY-1993】星际战争

    传送门:P3324 [SDOI2015]星际战争 二分 + 最大流 Solution\mathfrak{Solution}Solution 1 不可否认,看到题面的第一反应是构造二分图,让武器在一边, ...

  6. XMU2018摸底测试 星际战争I(Normal)

    星际战争I(Normal) Description 现在是星际战争的第一个年头. 绝地武士躲藏在一颗大小为N \times MN×M 的矩形二维星球上. 贸易联盟想要消灭他们,但是经历过之前的战斗之后 ...

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

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

  8. 「BZOJ 3994」「SDOI 2015」约数个数和「莫比乌斯反演」

    题意 设\(d(x)\)为\(x\)的约数个数,求\(\sum_{i=1}^{n}\sum_{j=1}^{m}d(ij)\). 题解 首先证个公式: \[d(ij) = \sum_{x|i}\sum_ ...

  9. bzoj3993 [SDOI2015]星际战争

    Description \(3333\) 年,在银河系的某星球上, X军团和Y军团正在激烈地作战.在战斗的某一阶段,Y军团一共派遣了 \(N\) 个巨型机器人进攻X军团的阵地,其中第i个巨型机器人的装 ...

  10. [SDOI2015]星际战争

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

最新文章

  1. 讯飞语音的中的bug用户校验失败
  2. IDEA 中比较骚后缀补全技巧!你可能没用过
  3. pinctrl框架【转】
  4. Invalid bound statement (not found): com.xsw.dao.CategoryDao.getCategoryById] with root cause
  5. android -------- NDK 入门指南
  6. 文献记录(part80)--基于平均互信息的最优社区发现方法
  7. HDU 4635(强连通分量分解
  8. NPOI 导出Excel 数据方式
  9. CSS相关知识点:6种清除浮动和BFC
  10. mongoexport导出mongodb数据库中的数据
  11. 迭代器 in C++
  12. java字符转为数字_java判断字符串是否可以转为数字
  13. 2020.7.20-每日进步
  14. 关于机器学习的十个实例
  15. 序列号大全,装机必备
  16. 如何做好科研实验记录
  17. WIN10系统盘安装提示MBR分区问题解决步骤
  18. windows无法自动修复此计算机 鼠标,win10出现你的账户已被停用。开机按F8进不去。自动修复鼠标和键盘都用...
  19. 防弹玻璃为啥会被钢球砸碎?这就是一道高中物理题!
  20. 第10节_我的日记本开发手记(10)——使用自定义图标字体

热门文章

  1. 2022推荐大宗商品行业研究报告投资策略产业发展前景市场行情分析(附件中为网盘链接,报告持续更新)
  2. python数据可视化——曲线图
  3. ruby + cucumber环境搭建
  4. 第3节:alphapose环境配置
  5. 远程控制android盒子,电视盒子ADB教程 通过ADB远程安装应用方法(2)
  6. git远程代码回滚_git 远程分支回滚
  7. vue批量打包下载图片
  8. PyG 中Message Passing机制详解
  9. 网页设计如何排成一列_网页排版设计中对齐技巧能让网站产生高端感(下)
  10. Online Convex Making Gradient Descent Optimal for Strongly Convex Stochastic Optimization