bzoj3993 [SDOI2015]星际战争

Description

---

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，输入数据保证X军团一定能摧毁Y军团的所有巨型机器人

Solution

---

滞涨选手总是弄不清n和m的区别（比如我）

容易发现一个激光造成的伤害=单位时间伤害*时间，那么考虑二分答案时间t跑最大流就行
一个技巧就是不用小数直接乘10000倍这样
滞涨选手总是弄不清n和m的区别（比如我）

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))

typedef long long LL;
const LL INF=50000000000000LL;
const int N=205;
const int E=20005;

struct edge {int x,y; LL w; int next;} e[E];

int queue[N],head,tail;
int dis[N],can[N][N];
int ls[N],edCnt;
LL a[N],b[N],tot;

void addEdge(int x,int y,LL w) {
    e[++edCnt]=(edge) {x,y,w,ls[x]}; ls[x]=edCnt;
    e[++edCnt]=(edge) {y,x,0,ls[y]}; ls[y]=edCnt;
}

bool bfs(int st,int ed) {
    head=1; tail=0; queue[++tail]=st;
    fill(dis,-1); dis[st]=1;
    while (head<=tail) {
        int now=queue[head++];
        for (int i=ls[now];i;i=e[i].next) {
            if (e[i].w>0&&dis[e[i].y]==-1) {
                dis[e[i].y]=dis[now]+1;
                if (e[i].y==ed) return true;
                queue[++tail]=e[i].y;
            }
        }
    }
    return false;
}

LL find(int now,int ed,LL mn) {
    if (now==ed||!mn) return mn;
    LL ret=0;
    for (int i=ls[now];i;i=e[i].next) {
        if (e[i].w>0&&dis[now]+1==dis[e[i].y]) {
            LL d=find(e[i].y,ed,std:: min(mn-ret,e[i].w));
            e[i].w-=d; e[i^1].w+=d;
            ret+=d;
            if (ret==mn) break;
        }
    }
    return ret;
}

LL dinic(int st,int ed) {
    LL ret=0;
    while (bfs(st,ed)) {
        ret+=find(st,ed,INF);
    }
    return ret;
}

bool check(int n,int m,LL lim) {
    fill(ls,0); edCnt=1;
    rep(i,1,m) addEdge(0,i,(LL)b[i]*lim);
    rep(i,1,n) addEdge(i+m,n+m+1,a[i]*10000);
    rep(i,1,m) rep(j,1,n) {
        if (can[i][j]) addEdge(i,j+m,INF);
    }
    LL mxFlow=dinic(0,n+m+1);
    if (mxFlow>=tot) return true;
    return false;
}

void solve(int n,int m) {
    LL l=0,r=INF;
    while (l<=r) {
        LL mid=(l+r)/2;
        if (check(n,m,mid)) r=mid-1;
        else l=mid+1;
    }
    double prt=(double)(r)/10000.0;
    printf("%.4lf\n", prt);
}

int main(void) {
    int n,m; scanf("%d%d",&n,&m);
    rep(i,1,n) scanf("%d",&a[i]);
    rep(i,1,m) scanf("%d",&b[i]);
    rep(i,1,n) tot+=a[i]*10000;
    rep(i,1,m) rep(j,1,n) scanf("%d",&can[i][j]);
    solve(n,m);
    return 0;
}