BZOJ_3993_[SDOI2015]星际战争_二分+网络流

Description

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

Input

第一行,两个整数,N、M。

第二行,N个整数,A1、A2…AN。
第三行,M个整数,B1、B2…BM。
接下来的M行,每行N个整数,这些整数均为0或者1。这部分中的第i行的第j个整数为0表示第i个激光武器不可以攻击第j个巨型机器人,为1表示第i个激光武器可以攻击第j个巨型机器人。

Output

一行,一个实数,表示X军团要摧毁Y军团的所有巨型机器人最少需要的时间。输出结果与标准答案的绝对误差不超过10-3即视为正确。

Sample Input

2 2
3 10
4 6
0 1
1 1

Sample Output

1.300000

HINT

【样例说明1】

战斗开始后的前0.5秒,激光武器1攻击2号巨型机器人,激光武器2攻击1号巨型机器人。1号巨型机器人被完全摧毁,2号巨型机器人还剩余8的装甲值;
接下来的0.8秒,激光武器1、2同时攻击2号巨型机器人。2号巨型机器人被完全摧毁。
对于全部的数据,1<=N, M<=50,1<=Ai<=105,1<=Bi<=1000,输入数据保证X军团一定能摧毁Y军团的所有巨型机器人

 首先二分答案转化成判定问题。
设二分的答案为x,S向每个武器连x*b[i]的边,每个武器向机器人连inf的边,机器人向T连a[i]的边。
求最大流和sum{a[i]}比较。
注意流量需要用浮点数。
 
代码:
  1. #include <stdio.h>
  2. #include <string.h>
  3. #include <algorithm>
  4. #include <stdlib.h>
  5. using namespace std;
  6. typedef double f2;
  7. #define N 150
  8. #define M 300050
  9. #define S (n+m+1)
  10. #define T (n+m+2)
  11. #define inf 10000000
  12. f2 flow[M];
  13. int head[N],to[M],nxt[M],cnt,n,m,a[N],b[N];
  14. int dep[N],Q[N],l,r,sum,map[N][N];
  15. f2 fabs(f2 x){return x>0?x:-x;}
  16. inline void add(int u,int v,f2 w) {
  17. to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; flow[cnt]=w;
  18. to[++cnt]=u; nxt[cnt]=head[v]; head[v]=cnt; flow[cnt]=0;
  19. }
  20. bool bfs() {
  21. memset(dep,0,sizeof(dep));
  22. dep[S]=1;int i;l=r=0;Q[r++]=S;
  23. while(l<r) {
  24. int x=Q[l++];
  25. for(i=head[x];i;i=nxt[i]) if(!dep[to[i]]&&flow[i]) {
  26. dep[to[i]]=dep[x]+1; if(to[i]==T) return 1; Q[r++]=to[i];
  27. }
  28. }
  29. return 0;
  30. }
  31. f2 dfs(int x,f2 mf) {
  32. if(x==T) return mf;
  33. int i;
  34. f2 nf=0;
  35. for(i=head[x];i;i=nxt[i]) {
  36. if(dep[to[i]]==dep[x]+1&&flow[i]) {
  37. f2 tmp=dfs(to[i],min(mf-nf,flow[i]));
  38. nf+=tmp;
  39. flow[i]-=tmp;
  40. flow[i^1]+=tmp;
  41. if(nf==mf) break;
  42. }
  43. }
  44. dep[x]=0;
  45. return nf;
  46. }
  47. bool check(f2 x) {
  48. memset(head,0,sizeof(head)); cnt=1;
  49. int i,j;
  50. for(i=1;i<=n;i++) add(i,T,a[i]);
  51. for(i=1;i<=m;i++) add(S,i+n,b[i]*x);
  52. for(i=1;i<=m;i++) {
  53. for(j=1;j<=n;j++) if(map[i][j]) {
  54. add(i+n,j,inf);
  55. }
  56. }
  57. f2 ans=0,f;
  58. while(bfs()) while(1) {
  59. f2 f=dfs(S,inf);
  60. if(fabs(f)<1e-6) break;
  61. ans+=f;
  62. }
  63. return fabs(ans-sum)<1e-6;
  64. }
  65. int main() {
  66. scanf("%d%d",&n,&m);
  67. int i,j;
  68. for(i=1;i<=n;i++) scanf("%d",&a[i]),sum+=a[i];
  69. for(i=1;i<=m;i++) scanf("%d",&b[i]);
  70. for(i=1;i<=m;i++) {
  71. for(j=1;j<=n;j++) {
  72. scanf("%d",&map[i][j]);
  73. }
  74. }
  75. f2 L=0,R=1e5;
  76. for(i=1;i<=70;i++) {
  77. f2 mid=(L+R)/2;
  78. if(check(mid)) R=mid;
  79. else L=mid;
  80. }
  81. printf("%.4lf\n",L);
  82. }

BZOJ_3993_[SDOI2015]星际战争_二分+网络流的更多相关文章

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

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

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

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

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

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

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

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

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

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

  6. 洛谷P3324 [SDOI2015]星际战争

    题目:洛谷P3324 [SDOI2015]星际战争 思路: 类似<导弹防御塔>,因为题目保证有解,花费时间小于最终答案时一定无法消灭所有敌人,只要花费时间大于等于最终答案都可以消灭所有敌人 ...

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

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

  8. 3993: [SDOI2015]星际战争

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

  9. 网络流(二分):BZOJ 3993: [SDOI2015]星际战争

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

随机推荐

  1. javah tool for Android Native Application

    javah可以在Eclipse中配置成为External Tools,选择External Tools Configurations,配置如下,经过测试通过. Location: ${system_p ...

  2. 框架学习:ibatis框架的结构和分析

    由于最近一段时间比较忙碌,<框架学习>系列的文章一直在搁浅着,最近开始继续这个系列的文章更新. 在上篇文章中我们说到了hibernate框架,它是一种基于JDBC的主流持久化框架,是一个优 ...

  3. .net Entity Framework初识1

    利用EF可以直接操纵数据库,在一些简单的项目里甚至完全不用写sql. 一 code first 1.在web.config中设置连接字符串 这一步可以省略.如果跳过这一步,程序会默认生成一个可用的连接 ...

  4. openssh升级的坑爹之路

    安装Zlib http://zlib.net/zlib-1.2.8.tar.gz tar -zxvf zlib-1.2.8.tar.gz cd zlib-1.2.8 ./configure --sha ...

  5. ffmpeg 时间戳处理

    视频的显示和存放原理 对于一个电影,帧是这样来显示的:I B B P.现在我们需要在显示B帧之前知道P帧中的信息.因此,帧可能会按照这样的方式来存储:IPBB.这就是为什么我们会有一个解码时间戳和一个 ...

  6. oracle数据库的备份与还原(本地及远程操作)

    数据的导出 exp qh/qh@qh  file='d:\backup\qh\qh20060526.dmp'  grants=y  full=n  1 将数据库TEST完全导出,用户名system 密 ...

  7. 自定义用户认证(继承django的)

    1.在app下创建一个自己用户认证文件,文件名随意,记得为.py文件 2.编辑该userauth.py文件 #!/usr/bin/env python #coding:utf-8 from djang ...

  8. 面向对象的WebAPI框架XXL-HEX

    <面向对象的WebAPI框架XXL-HEX>    一.简介 1.1 概述 XXL-HEX 是一个简单易用的WebAPI框架, 拥有 "面向对象.数据加密.跨语言" 的 ...

  9. GitHub学习笔记:远程端的操控

    对于远端,当你新建一个项目的时候,需要在网页处新建,在新建项目的页面,会有一段提示你上传本地项目到此远端方法的代码,直接拷贝粘贴到git shell就可以解决问题,不再详述. 当你把代码上传到一个已经 ...

  10. Java面向对象(二、继承)

    Java 继承 继承的概念 继承是java面向对象编程技术的一块基石,因为它允许创建分等级层次的类. 继承就是子类继承父类的特征和行为,使得子类对象(实例)具有父类的实例域和方法,或子类从父类继承方法 ...