[SDOI2015]星际战争

一、题目

点此看题

二、解法

可以把激光看成$X$部，把机器人看成$Y$部，考虑二分答案后建图。

• 建一个起点，把所有的激光连向起点，容量为$t\times b[i]$，也就是把伤害都给激光，然后让它去伤害机器人。
• 建一个终点，把所有机器人连向终点，容量为$a[i]$，表示一个机器人能最多接受$b[i]$的伤害，检查受到的伤害是否等于$a[i]$之和，也就是最大流。
• 每一个机器人和激光之间连一条边，不需要限流，直接为$inf$。

所以时间复杂度$O(n^4\log )$，实际远远跑不满。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
#define int long long
#define inf (1ll<<40)
const int MAXN = 105;
const int C = 1000;
using namespace std;
int read()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,S,T,tot,cur[MAXN],a[MAXN],b[MAXN],f[MAXN],dis[MAXN],c[MAXN][MAXN];
int ans,sum;
queue<int> q;
struct edge
{
	int v,c,next;
}e[MAXN*MAXN];
void add_edge(int u,int v,int c)
{
	e[++tot]=edge{v,c,f[u]},f[u]=tot;
	e[++tot]=edge{u,0,f[v]},f[v]=tot;
}
int bfs()
{
	memset(dis,0,sizeof dis);
	dis[S]=1;
	q.push(S);
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=f[u];i;i=e[i].next)
		{
			int v=e[i].v;
			if(e[i].c>0 && !dis[v])
			{
				dis[v]=dis[u]+1;
				q.push(v);
			}
		}
	}
	if(!dis[T]) return 0;
	return 1;
}
int dfs(int u,int ept)
{
	if(u==T) return ept;
	int flow=0,tmp=0;
	for(int &i=cur[u];i;i=e[i].next)
	{
		int v=e[i].v;
		if(dis[u]+1==dis[v] && e[i].c>0)
		{
			tmp=dfs(v,min(e[i].c,ept));
			if(!tmp) continue;
			ept-=tmp;
			e[i].c-=tmp;
			e[i^1].c+=tmp;
			flow+=tmp;
			if(!ept) break;
		}
	}
	return flow;
}
int check(int x)
{
	tot=1;int ans=0;
	memset(f,0,sizeof f);
	for(int i=1;i<=m;i++)
		for(int j=1;j<=n;j++)
			if(c[i][j])
				add_edge(1+i,1+m+j,inf);
	for(int i=1;i<=m;i++)
		add_edge(1,1+i,x*b[i]);
	for(int i=1;i<=n;i++)
		add_edge(1+m+i,n+m+2,a[i]);
	S=1;T=n+m+2;
	while(bfs())
	{
		for(int i=1;i<=n+m+2;i++)
			cur[i]=f[i];
		ans+=dfs(S,inf);
	}
	return ans;
}
void solve(int l,int r)
{
	if(l>r) return ;
	int mid=(l+r)/2;
	if(check(mid)==sum)
	{
		ans=mid;
		solve(l,mid-1);
	}
	else
		solve(mid+1,r);
}
signed main()
{
	n=read();m=read();
	for(int i=1;i<=n;i++)
	{
		a[i]=read()*C;
		sum+=a[i];
	}
	for(int i=1;i<=m;i++)
		b[i]=read();
	for(int i=1;i<=m;i++)
		for(int j=1;j<=n;j++)
			c[i][j]=read();
	solve(1,inf);
	printf("%lf\n",ans*1.0/C);
}