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;
}