题目描述:
有N个防御机器人,每个机器人有一个防御值K
有M个炮塔,每个炮塔单位时间内可以造成X的伤害
这M个炮塔可以攻击一些机器人,也有一些不能攻击。
求把这些机器人全部打败的时间最少是多少?
题目分析:
看数据范围和题目要求,很明显就是网络流的模型了.
首先二分一个答案 t
那么每个炮塔最多可以造成 t*X 的伤害
由 S -> 炮塔 连容量为 t*X 的 边
对于炮塔可以攻击的机器人,
炮塔 -> 机器人 连容量为 inf 的边
对于每个机器人
机器人 -> T 连一条容量为 K 的 边
很明显,当所有机器人 -> T的边全部满流的时候,当前答案符合要求
答案要求我们不超过
10−3
10
−
3
的精度误差 我们即可以二分实数,也可以把机器人的血量放大10000倍 相应的得出的时间也放大了10000倍,除一下即可
注意二分上界大一点TAT
题目链接:
Ac 代码:
#include <cstdio>
#include <iostream>
#include <queue>
#include <cstring>
#define int long long
#define il inline
using namespace std;
const int inf=1e17+100;
const int maxm=70*70;
int head[1100],to[maxm*2],cap[maxm*2],net[maxm*2],deep[1100],cnt=1;
il void add(int x,int y,int c){cnt++,to[cnt]=y,cap[cnt]=c,net[cnt]=head[x],head[x]=cnt;}
queue <int> dl;
int n,m,a[60],b[60],sum;
bool at[60][60];
il bool BFS(int s,int t)
{
while(!dl.empty()) dl.pop();
memset(deep,-1,sizeof(deep));
dl.push(s),deep[s]=0;
while(!dl.empty())
{
int x=dl.front();dl.pop();
for(int i=head[x];i;i=net[i])
if(cap[i]>0&&deep[to[i]]==-1)
dl.push(to[i]),deep[to[i]]=deep[x]+1;
}
return deep[t]==-1?0:1;
}
int dfs(int now,int flow,int t)
{
if(now==t) return flow;
int w,used=0;
for(int i=head[now];i;i=net[i])
{
int v=to[i];
if(deep[v]==deep[now]+1&&cap[i])
{
w=dfs(v,min(flow-used,cap[i]),t);
cap[i]-=w;
cap[i^1]+=w;
used+=w;
if(used==flow) return flow;
}
}
if(!used) deep[now]=-1;
return used;
}
il int dinic(int s,int t)
{
int maxflow=0;
while(BFS(s,t)) maxflow+=dfs(s,inf,t);
return maxflow;
}
inline void adx(int x,int y,int cax)
{
add(x,y,cax),add(y,x,0);
}
il bool check(int mid)
{
int s=0,t=n+m+1;
memset(head,0,sizeof(head));
cnt=1;
for(int i=1;i<=m;i++)
adx(s,i,b[i]*mid);
for(int i=1;i<=n;i++)
adx(i+m,t,a[i]*10000);
for(int i=1;i<=m;i++)
for(int j=1;j<=n;j++)
if(at[i][j])
adx(i,j+m,inf);
return dinic(s,t)==sum;
}
main()
{
scanf("%lld%lld",&n,&m);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]),sum+=a[i]*10000;
for(int i=1;i<=m;i++)
scanf("%lld",&b[i]);
for(int i=1;i<=m;i++)
for(int j=1;j<=n;j++)
scanf("%lld",&at[i][j]);
int l=0,r=100000000;
int ans=0;
while(l<=r)
{
int mid=(l+r)>>1;
if(check(mid)) ans=mid,r=mid-1;
else l=mid+1;
}
printf("%.4lf\n",ans/10000.0);
return 0;
}