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.当一个巨型 ...
随机推荐
- Gulp基础知识
首先,我们需要了解Gulp能做些什么? 编译 sass sass是什么?(使CSS可以用编程的方式写,加快我们开发的速度) ...
- JAVA 综合面试题
JAVA 综合面试题 2007-08-12 目录 TOC \o "1-3" \h \z \u Java面试题整理 9 Java面向对象 9 1. super()与this()的区别 ...
- 并发编程(四):atomic
本篇博客我们主要讲述J.U.C包下的atomic包,在上篇博客"并发模拟"的最后,我们模拟高并发的情形时出现了线程安全问题,怎么解决呢?其实解决的办法有很多中,如直接在add()方 ...
- Java反射之修改常量值
1. 通过反射修改常量的值 package com.blueStarWei.invoke; import java.lang.reflect.Field; public class ModifyFin ...
- Day20 Django的使用_基础
老师网址: https://www.cnblogs.com/yuanchenqi/articles/7652353.html 1,复习上级课,一对一,一对多,多对多的使用 models.py: cla ...
- springMVC中添加restful 风格
RESTful架构:是一种设计的风格,并不是标准,只是提供了一组设计原则和约束条件,也是目前比较流行的一种互联网软件架构.它结构清晰.符合标准.易于理解.扩展方便,所以正得到越来越多网站的采用. 关于 ...
- 算法竞赛之递归——输出1-n的所有排列
本文是博主原创文章,未经允许不得转载.我的csdn博客也同步发布了此文, 链接 https://blog.csdn.net/umbrellalalalala/article/details/79792 ...
- ArcCore重构-目标文件结构化
基于官方arc-stable-9c57d86f66be,AUTOSAR版本3.1.5 基本问题 3. 编译系统中所有代码文件通过搜索路径(VPATH)中搜索,存在名称污染问题,需加入路径信息: ...
- Spring+Redis的部署与Redis缓存使用示例
由于项目的业务需要,这两天折腾了一下Spring-redis配置,有了前面用Spring托管hibernate的经验,这次可以说是顺风顺水,大概说一下流程. ubuntu 安装 redis sudo ...
- Total Command使用笔记
一.快键键(基于水晶2右)以下数字小键盘无效 Tab 左右窗口切换 Ctrl+d 进入工作目录ctrl+d+数字 指定目录alt+←/→ 后退/前进目录ctrl+\ 跳转到根目录Ctrl+b,不分层级 ...