【SDOI2015】星际战争 网络流

题目描述

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

数据范围

1<=N, M<=50，1<=Ai<=105，1<=Bi<=1000

样例输入

5 4
30 24 11 48 27
4 3 6 1
0 0 0 1 1
1 0 0 1 0
0 1 1 0 1
0 1 0 0 1

样例输出

11.142858

解题思路

二分答案，网络流判定，如果满流既可以摧毁所有机器人。

代码

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#define eps 1e-5
#define Maxn 500005//数组都是随便开的。
#define Maxe 500005
using namespace std;
inline int Getint(){int x=0,f=1;char ch=getchar();while('0'>ch||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while('0'<=ch&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;}
int S,T,N,cnt=0,dis[Maxn],GAP[Maxn],h[Maxn];
struct node{int to,next;double v;int pair;}e[Maxe];
void AddEdge(int x,int y,double v,int pair){e[cnt]=(node){y,h[x],v,pair};h[x]=cnt;}
void AddEdge(int x,int y,double v){AddEdge(x,y,v,++cnt+1);AddEdge(y,x,0,++cnt-1);}
double SAP(int x,double Maxflow){
    if(x==T)return Maxflow;
    double tmp=Maxflow;
    for(int p=h[x];p;p=e[p].next){
        int y=e[p].to;
        double flow=min(tmp,e[p].v);
        if(flow&&dis[x]==dis[y]+1){
            double ret=SAP(y,flow);
            tmp-=ret;
            e[p].v-=ret;
            e[e[p].pair].v+=ret;
            if(!tmp||dis[S]==N)return Maxflow-tmp;
        }
    }
    if(--GAP[dis[x]]==0)dis[S]=N;
    else GAP[++dis[x]]++;
    return Maxflow-tmp;
}
double SAP(){
    memset(GAP,0,sizeof(GAP));
    memset(dis,0,sizeof(dis));
    GAP[0]=N;
    double Ans=0;
    while(dis[S]<N)Ans+=SAP(S,1<<30);
    return Ans;
}
int a[55],b[55],n,m,Sum=0;
bool Map[55][55];
bool Check(double x){
    memset(h,0,sizeof(h));
    memset(e,0,sizeof(e));
    cnt=0;
    for(int i=1;i<=n;i++)AddEdge(i,T,a[i]);
    for(int i=1;i<=m;i++)AddEdge(S,i+n,x*b[i]);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            if(Map[i][j])
                AddEdge(i+n,j,1<<30);
    if(fabs(SAP()-Sum)<=eps)return true;
    return false;
}
int main(){
    n=Getint(),m=Getint();
    for(int i=1;i<=n;i++)Sum+=a[i]=Getint();
    for(int i=1;i<=m;i++)b[i]=Getint();
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            Map[i][j]=Getint();
    S=0,T=n+m+1,N=T+1;
    double L=0,r=1e5;
    while(L+eps<=r){
        double mid=(L+r)/2;
        if(Check(mid))r=mid;
        else L=mid;
    }
    printf("%.6lf",(L+r)/2);
    return 0;
}

转载于:https://www.cnblogs.com/Cedric341561/p/6811017.html