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
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军团的所有巨型机器人
一眼二分+最大流
莫非网络流代码抄多了自带升级效果？？？
ST---->i 能够造成的伤害
i---->j 能打的机器人，INF
j----> ed 护甲值
然后枚举时间就OK了
丧心病狂这道题竟然卡我精度
代码如下：
#include<cmath> #include<cstdio> #include<cstring> #include<cstdlib> #include<algorithm> #define eps 1e-7 using namespace std; double INF=(1<<29); struct dicnic{ int x,y,next,other; double c; }a[2100000];int len,last[210000]; void ins(int x,int y,double c) { int k1,k2; k1=++len; a[len].x=x;a[len].y=y;a[len].c=c; a[len].next=last[x];last[x]=len; k2=++len; a[len].x=y;a[len].y=x;a[len].c=0; a[len].next=last[y];last[y]=len; a[k1].other=k2; a[k2].other=k1; } int st,ed,head,tail; int h[210000],list[210000]; bool bt_h() { memset(h,0,sizeof(h));h[st]=1; list[1]=st;head=1;tail=2; while(head!=tail) { int x=list[head]; for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(h[y]==0&&a[k].c>0) { h[y]=h[x]+1; list[tail++]=y; } } head++; } if(h[ed]>0)return true; return false; } double findflow(int x,double f) { if(x==ed)return f; double s=0,t; for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(h[y]==h[x]+1&&a[k].c>0&&s<f) { s+=(t=findflow(y,min(a[k].c,f-s))); a[k].c-=t;a[a[k].other].c+=t; } } if(s==0)h[x]=0; return s; } int n,m; double b[51]; double c[51]; int f[51][51]; bool check(double x) { double sum=0.0; len=0;memset(last,0,sizeof(last)); for(int i=1;i<=m;i++)ins(st,i,x*c[i]); for(int i=1;i<=m;i++) for(int j=1;j<=n;j++) if(f[i][j]==1) ins(i,j+m,INF); for(int i=1;i<=n;i++)ins(i+m,ed,b[i]),sum+=b[i]; double ans=0.0; while(bt_h()==true)ans+=findflow(st,INF); if(abs(ans-sum)<=eps)return true; return false; } int main() { scanf("%d%d",&n,&m);st=0,ed=n+m+1; len=0;memset(last,0,sizeof(last)); for(int i=1;i<=n;i++)scanf("%lf",&b[i]); for(int i=1;i<=m;i++)scanf("%lf",&c[i]); for(int i=1;i<=m;i++) for(int j=1;j<=n;j++) scanf("%d",&f[i][j]); double l=0.0,r=10000000.0;double ans; while(l<=r) { double mid=(l+r)/2; if(check(mid)==true){r=mid-eps;ans=mid;} else l=mid+eps; } printf("%.6lf\n",ans); return 0; }
by_lmy
