題目傳送門

網路流的用法又漲姿勢了——網路流可以用於判斷答案的可行性。

這題我們首先考慮建圖：

1. 第i個機器人向超級匯點連一條流量為ai的邊。
2. 超級源點向第i個武器連一條流量為?的邊。（?表示流量暫時未定）
3. 如果第i個武器可以打到第j個機器人，就連一條流量為+∞的邊。

然後考慮二分一個時間mid作為答案，這時考慮第2條建圖，這時的?就等於bi×mid。

如果最大流等於所有機器人的裝甲值，那麼當前的時間mid就是可行的。

附上AC程式碼：

#include <cstdio>
#include <cctype>
#include <cstring>
 

#include <queue>
#include <cmath>
using namespace std;

const int N=110;
struct side{
    int to;double w;int nt;
}s[N*N];
int n,m,x,b[N],map[N][N],st,ed,h[N],num,sum,d[N],cur[N];
queue <int> que;

inline char nc(void){
    static char ch[100010],*p1=ch,*p2=ch;
    return p1==p2&&(p2=(p1=ch)+fread(ch,1
 
,100010,stdin),p1==p2)?EOF:*p1++;
}

inline void read(int &a){
    static char c=nc();int f=1;
    for (;!isdigit(c);c=nc()) if (c=='-') f=-1;
    for (a=0;isdigit(c);a=(a<<3)+(a<<1)+c-'0',c=nc());
    return (void)(a*=f);
}

inline void add(int x,int y,int w){
    s[num]=(side){y,w,h[x]},h[x]=num++;
    s[num]=(side){x,0
 
,h[y]},h[y]=num++;
}

inline int bfs(int st,int ed){
    while (!que.empty()) que.pop();
    memset(d,0,sizeof d),d[st]=1,que.push(st);
    while (!que.empty()){
        int p=que.front();que.pop();
        for (int i=h[p]; ~i; i=s[i].nt)
            if (s[i].w&&!d[s[i].to]) d[s[i].to]=d[p]+1,que.push(s[i].to);
        if (d[ed]) break;
    }
    return d[ed];
}

inline double so(int x,double w){
    if (x==ed) return w;
    double sum=0,f;
    for (int &i=cur[x]; ~i; i=s[i].nt)
        if (s[i].w&&d[s[i].to]==d[x]+1&&(f=so(s[i].to,min(w-sum,s[i].w)))){
            s[i].w-=f,s[i^1].w+=f,sum+=f;
            if (sum==w) return sum;
        }
    if (!sum) d[x]=0;
    return sum;
}

inline bool check(double mid){
    for (int i=0; i<num; i+=2) s[i].w+=s[i^1].w,s[i^1].w=0;
    for (int i=h[st]; ~i; i=s[i].nt) s[i].w=b[s[i].to-n]*mid;
    double ret=0;
    while (bfs(st,ed)) memcpy(cur,h,sizeof h),ret+=so(st,0x7fffffff);
    return fabs(ret-sum)<1e-6;
}

int main(void){
    read(n),read(m),st=0,ed=n+m+1,memset(h,-1,sizeof h);
    for (int i=1; i<=n; ++i) read(x),add(i,ed,x),sum+=x;
    for (int i=1; i<=m; ++i) read(b[i]),add(st,i+n,2e9);
    for (int i=1; i<=m; ++i)
        for (int j=1; j<=n; ++j)
            if (read(x),x) add(i+n,j,2e9);
    double l=0,r=2e9,ans=0,mid;
    while (r-l>1e-6)
        if (check(mid=(l+r)/2)) ans=mid,r=mid;
        else l=mid;
    return printf("%f\n",ans),0;
}