Luogu P3324 [SDOI2015]星际战争

二分+最大流

首先考虑二分答案

然后可以发现对于已知时间，判断是否可以将所有机器人摧毁可以用网络流

建立源点和汇点，源点向每一个激光武器连一条容量为$time*b[i]$的边，表示该激光武器在$time$时间下最多能产生的伤害为$time*b[i]$

每一个机器人向汇点连一条容量为$a[i]$的边，表示每一个机器人最多能承受的伤害

中间每一个激光武器向其能攻击的机器人连一条容量为$inf$的边

然后跑最大流，若可以跑出$\sum a[i]$，那么该时间可以将所有与机器人摧毁，更新二分边界

  1 #include <bits/stdc++.h>
  2 #define inf 1e9
  3 using namespace std;
  4 const int N=210,M=6100;
  5 int n,m,s,t,tot,first[N],nxt[M],point[M],nrl[N],d[N];
  6 int vi[N][N];
  7 double data[M],a[N],b[N],sum;
  8 void add_edge(int x,int y,double z)
  9 {
 10     tot++;
 11     nxt[tot]=first[x];
 12     first[x]=tot;
 13     point[tot]=y;
 14     data[tot]=z;
 15     tot++;
 16     nxt[tot]=first[y];
 17     first[y]=tot;
 18     point[tot]=x;
 19     data[tot]=0;
 20 }
 21 bool bfs()
 22 {
 23     queue <int> q;
 24     for (int i=s;i<=t;i++)
 25     {
 26         d[i]=inf;
 27         nrl[i]=first[i];
 28     }
 29     d[s]=0;q.push(s);
 30     while (!q.empty())
 31     {
 32         int x=q.front();
 33         q.pop();
 34         for (int i=first[x];i!=-1;i=nxt[i])
 35         {
 36             int u=point[i];
 37             if (data[i]>0 && d[u]>d[x]+1)
 38             {
 39                 d[u]=d[x]+1;
 40                 if (u==t) return true;
 41                 q.push(u);
 42             }
 43         }
 44     }
 45     return false;
 46 }
 47 double dfs(int x,double flow)
 48 {
 49     if (x==t) return flow;
 50     double sum=0;
 51     for (int i=nrl[x];i!=-1;i=nxt[i])
 52     {
 53         int u=point[i];
 54         if (data[i]>0 && d[u]==d[x]+1)
 55         {
 56             double tmp=dfs(u,min(flow,data[i]));
 57             flow-=tmp;sum+=tmp;
 58             data[i]-=tmp;data[i^1]+=tmp;
 59             if (flow<=0) return sum;
 60         }
 61         nrl[x]=nxt[i];
 62     }
 63     return sum;
 64 }//使用dinic跑最大流
 65 bool check(double mid)
 66 {
 67     tot=-1;
 68     memset(first,-1,sizeof(first));
 69     memset(nxt,-1,sizeof(nxt));
 70     s=0;t=n+m+1;//同上建边
 71     for (int i=1;i<=m;i++) add_edge(s,i,mid*b[i]);
 72     for (int i=1;i<=n;i++) add_edge(m+i,t,a[i]);
 73     for (int i=1;i<=m;i++)
 74     {
 75         for (int j=1;j<=n;j++)
 76         {
 77             if (vi[i][j]) add_edge(i,m+j,inf);
 78         }
 79     }
 80     double ans=0;
 81     while (bfs()) ans+=dfs(s,inf);
 82     return ans==sum;//判断是否跑满
 83 }
 84 int main()
 85 {
 86     scanf("%d%d",&n,&m);
 87     for (int i=1;i<=n;i++) scanf("%lf",&a[i]);
 88     for (int i=1;i<=m;i++) scanf("%lf",&b[i]);
 89     for (int i=1;i<=m;i++)
 90     {
 91         for (int j=1;j<=n;j++)
 92           scanf("%d",&vi[i][j]);
 93     }
 94     sum=0;
 95     for (int i=1;i<=n;i++) sum+=a[i];
 96     double l=0,r=1e8;
 97     while ((r-l)>=1e-8)
 98     {
 99         double mid=(l+r)/2;
100         if (check(mid)) r=mid;
101         else l=mid;
102     }
103     printf("%.6lf\n",r);
104 }