【SDOI2015】【BZOJ3993】星际战争
2016-01-25 11:57
375 查看
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军团的所有巨型机器人
Source
Round 1 感谢yts1999上传
二分费用最大流判定
当时年轻的我知道是二分费用最大流判定但是不会实数二分最后跪了..当年我怎么这么傻逼(现在一样傻逼..
数组大小我随手开的..因为一开始RE了几发后来索性开的特别大..
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军团的所有巨型机器人
Source
Round 1 感谢yts1999上传
二分费用最大流判定
当时年轻的我知道是二分费用最大流判定但是不会实数二分最后跪了..当年我怎么这么傻逼(现在一样傻逼..
数组大小我随手开的..因为一开始RE了几发后来索性开的特别大..
#include<iostream> #include<cstdio> #include<cstring> #include<cmath> #include<algorithm> #define eps 1e-7 #define GET (ch>='0'&&ch<='9') #define MAXINT 0x7fffffff #define MAXN 1100 #define MAXM 10010 using namespace std; int n,m,top,s,t; int dis[MAXN],num[MAXN]; int a[MAXN],b[MAXN],map[MAXN][MAXN]; double sum; struct edge { int to; double c,f; edge *next,*rev; }e[MAXM<<1],*prev[MAXN]; void insert(int u,int v,double c) { e[++top].to=v;e[top].next=prev[u];prev[u]=&e[top];e[top].c=c; e[++top].to=u;e[top].next=prev[v];prev[v]=&e[top];e[top].c=0; prev[u]->rev=prev[v];prev[v]->rev=prev[u]; } void in(int &x) { char ch=getchar();x=0; while (!GET) ch=getchar(); while (GET) x=x*10+ch-'0',ch=getchar(); } void init() { int que[MAXN],head=0,tail=0; for (int i=0;i<=n+m+1;i++) dis[i]=MAXINT,num[i]=0; dis[t]=0;que[tail++]=t;num[0]=1; while (head!=tail) { int x=que[head++]; for (edge *i=prev[x];i;i=i->next) if (((i->rev)&&(i->rev->c<=eps))||(dis[i->to]<MAXINT)) {} else dis[i->to]=dis[x]+1,num[dis[i->to]]++,que[tail++]=i->to; } } double maxflow() { int st=s;double ret=0; edge *E[MAXN],*rep[MAXN]; for (int i=0;i<=n+m+1;i++) E[i]=prev[i]; while (dis[s]<n+m+2) { if (st==t) { double minn=1e9; for (int i=s;i!=t;i=E[i]->to) minn=min(minn,E[i]->c); for (int i=s;i!=t;i=E[i]->to) E[i]->c-=minn,E[i]->f+=minn,E[i]->rev->c+=minn,E[i]->rev->f-=minn; ret+=minn;st=s; } edge *i; for (i=E[st];i;i=i->next) if (i->c>0&&dis[st]==dis[i->to]+1) break; if (i) E[st]=i,rep[i->to]=i->rev,st=i->to; else { if ((--num[dis[st]])==0) break; E[st]=prev[st];int minn=n+m+2; for (edge *I=prev[st];I;I=I->next) if (I->c>eps) minn=min(minn,dis[I->to]); dis[st]=minn+1;++num[dis[st]]; if (st!=s) st=rep[st]->to; } } return ret; } void rebuild(double mid) { memset(e,0,sizeof(e));memset(prev,0,sizeof(prev));top=0; for (int i=1;i<=m;i++) insert(s,i,mid*(double)b[i]); for (int i=1;i<=m;i++) for (int j=1;j<=n;j++) if (map[i][j]) insert(i,j+m,1e9); for (int i=1;i<=n;i++) insert(i+m,t,a[i]); } int main() { in(n);in(m);s=0;t=n+m+1; for (int i=1;i<=n;i++) in(a[i]),sum+=(double)a[i]; for (int i=1;i<=m;i++) in(b[i]); for (int i=1;i<=m;i++) for (int j=1;j<=n;j++) in(map[i][j]); double l=1,r=1e6,mid=(l+r)/2; while (r-l>eps) { mid=(l+r)/2; rebuild(mid);init(); double t=maxflow();//printf("%.6f\n",t); if (fabs(t-sum)<eps) r=mid; else l=mid; } printf("%.6f\n",l); }
相关文章推荐
- Poj2638 网络流+最短路+二分答案
- BZOJ3275 Number (最小割)
- [笔记] 网络流-最大流 POJ-1273\HDU-4240
- 上下界网络流初探
- Edmonds-Karp 最大流 hdu 1532 Drained Ditches
- 网络流_poj1273
- POJ 1273 Drainage Ditches 最大流 dinic
- POJ1273-Drainage Ditches
- 【网络流】复杂的大门
- ACM/ICPC World Finals 2013 C Surely You Congest
- 网络流算法整理
- [BZOJ1797][AHOI2009][最大流][强连通分量]Mincut最小割
- [BZOJ2324][ZJOI2011][最小费用最大流]营救皮卡丘
- [BZOJ1834][ZJOI2010][最大流][最小费用最大流]网络扩容
- Topcoder SRM642 TaroCutting
- POJ2391解题报告
- uva 11248
- 【网络流之最小割模型】poj3469 BZOJ3144 UVA1212
- [网络流24题 #2]太空飞行计划问题
- [网络流24题 #18]分配问题