BZOJ_3993_[SDOI2015]星际战争_二分+网络流
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。
Output
一行,一个实数,表示X军团要摧毁Y军团的所有巨型机器人最少需要的时间。输出结果与标准答案的绝对误差不超过10-3即视为正确。
Sample Input
3 10
4 6
0 1
1 1
Sample Output
HINT
【样例说明1】
- #include <stdio.h>
- #include <string.h>
- #include <algorithm>
- #include <stdlib.h>
- using namespace std;
- typedef double f2;
- #define N 150
- #define M 300050
- #define S (n+m+1)
- #define T (n+m+2)
- #define inf 10000000
- f2 flow[M];
- int head[N],to[M],nxt[M],cnt,n,m,a[N],b[N];
- int dep[N],Q[N],l,r,sum,map[N][N];
- f2 fabs(f2 x){return x>0?x:-x;}
- inline void add(int u,int v,f2 w) {
- to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; flow[cnt]=w;
- to[++cnt]=u; nxt[cnt]=head[v]; head[v]=cnt; flow[cnt]=0;
- }
- bool bfs() {
- memset(dep,0,sizeof(dep));
- dep[S]=1;int i;l=r=0;Q[r++]=S;
- while(l<r) {
- int x=Q[l++];
- for(i=head[x];i;i=nxt[i]) if(!dep[to[i]]&&flow[i]) {
- dep[to[i]]=dep[x]+1; if(to[i]==T) return 1; Q[r++]=to[i];
- }
- }
- return 0;
- }
- f2 dfs(int x,f2 mf) {
- if(x==T) return mf;
- int i;
- f2 nf=0;
- for(i=head[x];i;i=nxt[i]) {
- if(dep[to[i]]==dep[x]+1&&flow[i]) {
- f2 tmp=dfs(to[i],min(mf-nf,flow[i]));
- nf+=tmp;
- flow[i]-=tmp;
- flow[i^1]+=tmp;
- if(nf==mf) break;
- }
- }
- dep[x]=0;
- return nf;
- }
- bool check(f2 x) {
- memset(head,0,sizeof(head)); cnt=1;
- int i,j;
- for(i=1;i<=n;i++) add(i,T,a[i]);
- for(i=1;i<=m;i++) add(S,i+n,b[i]*x);
- for(i=1;i<=m;i++) {
- for(j=1;j<=n;j++) if(map[i][j]) {
- add(i+n,j,inf);
- }
- }
- f2 ans=0,f;
- while(bfs()) while(1) {
- f2 f=dfs(S,inf);
- if(fabs(f)<1e-6) break;
- ans+=f;
- }
- return fabs(ans-sum)<1e-6;
- }
- int main() {
- scanf("%d%d",&n,&m);
- int i,j;
- for(i=1;i<=n;i++) scanf("%d",&a[i]),sum+=a[i];
- for(i=1;i<=m;i++) scanf("%d",&b[i]);
- for(i=1;i<=m;i++) {
- for(j=1;j<=n;j++) {
- scanf("%d",&map[i][j]);
- }
- }
- f2 L=0,R=1e5;
- for(i=1;i<=70;i++) {
- f2 mid=(L+R)/2;
- if(check(mid)) R=mid;
- else L=mid;
- }
- printf("%.4lf\n",L);
- }
BZOJ_3993_[SDOI2015]星际战争_二分+网络流的更多相关文章
- BZOJ3993 [SDOI2015]星际战争 【二分 + 网络流】
题目 3333年,在银河系的某星球上,X军团和Y军团正在激烈地作战.在战斗的某一阶段,Y军团一共派遣了N个巨型机器人进攻X军团的阵地,其中第i个巨型机器人的装甲值为Ai.当一个巨型机器人的装甲值减少到 ...
- bzoj3993: [SDOI2015]星际战争(二分+最大流)
题目描述 3333年,在银河系的某星球上,X军团和Y军团正在激烈地作战. 在战斗的某一阶段,Y军团一共派遣了N个巨型机器人进攻X军团的阵地,其中第i个巨型机器人的装甲值为Ai.当一个巨型机器人的装甲值 ...
- bzoj3993: [SDOI2015]星际战争(网络流)
3993: [SDOI2015]星际战争 题目:传送门 题解: 洛谷AC了,但是因为bzoj的spj有问题所以暂时没A 一道老题目了,二分时间然后网络流判断. 每次st-->武器连时间*攻击力 ...
- BZOJ 3993: [SDOI2015]星际战争 [二分答案 二分图]
3993: [SDOI2015]星际战争 题意:略 R1D2T1考了裸二分答案+二分图最大匹配... #include <iostream> #include <cstdio> ...
- 【BZOJ3993】[SDOI2015]星际战争 二分+最大流
[BZOJ3993][SDOI2015]星际战争 Description 3333年,在银河系的某星球上,X军团和Y军团正在激烈地作战.在战斗的某一阶段,Y军团一共派遣了N个巨型机器人进攻X军团的阵地 ...
- 洛谷P3324 [SDOI2015]星际战争
题目:洛谷P3324 [SDOI2015]星际战争 思路: 类似<导弹防御塔>,因为题目保证有解,花费时间小于最终答案时一定无法消灭所有敌人,只要花费时间大于等于最终答案都可以消灭所有敌人 ...
- 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: ...
- 网络流(二分):BZOJ 3993: [SDOI2015]星际战争
Description 3333年,在银河系的某星球上,X军团和Y军团正在激烈 地作战.在战斗的某一阶段,Y军团一共派遣了N个巨型机器人进攻X军团的阵地,其中第i个巨型机器人的装甲值为Ai.当一个巨型 ...
随机推荐
- javah tool for Android Native Application
javah可以在Eclipse中配置成为External Tools,选择External Tools Configurations,配置如下,经过测试通过. Location: ${system_p ...
- 框架学习:ibatis框架的结构和分析
由于最近一段时间比较忙碌,<框架学习>系列的文章一直在搁浅着,最近开始继续这个系列的文章更新. 在上篇文章中我们说到了hibernate框架,它是一种基于JDBC的主流持久化框架,是一个优 ...
- .net Entity Framework初识1
利用EF可以直接操纵数据库,在一些简单的项目里甚至完全不用写sql. 一 code first 1.在web.config中设置连接字符串 这一步可以省略.如果跳过这一步,程序会默认生成一个可用的连接 ...
- 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 ...
- ffmpeg 时间戳处理
视频的显示和存放原理 对于一个电影,帧是这样来显示的:I B B P.现在我们需要在显示B帧之前知道P帧中的信息.因此,帧可能会按照这样的方式来存储:IPBB.这就是为什么我们会有一个解码时间戳和一个 ...
- oracle数据库的备份与还原(本地及远程操作)
数据的导出 exp qh/qh@qh file='d:\backup\qh\qh20060526.dmp' grants=y full=n 1 将数据库TEST完全导出,用户名system 密 ...
- 自定义用户认证(继承django的)
1.在app下创建一个自己用户认证文件,文件名随意,记得为.py文件 2.编辑该userauth.py文件 #!/usr/bin/env python #coding:utf-8 from djang ...
- 面向对象的WebAPI框架XXL-HEX
<面向对象的WebAPI框架XXL-HEX> 一.简介 1.1 概述 XXL-HEX 是一个简单易用的WebAPI框架, 拥有 "面向对象.数据加密.跨语言" 的 ...
- GitHub学习笔记:远程端的操控
对于远端,当你新建一个项目的时候,需要在网页处新建,在新建项目的页面,会有一段提示你上传本地项目到此远端方法的代码,直接拷贝粘贴到git shell就可以解决问题,不再详述. 当你把代码上传到一个已经 ...
- Java面向对象(二、继承)
Java 继承 继承的概念 继承是java面向对象编程技术的一块基石,因为它允许创建分等级层次的类. 继承就是子类继承父类的特征和行为,使得子类对象(实例)具有父类的实例域和方法,或子类从父类继承方法 ...