【BZOJ】3993: [SDOI2015]星际战争

题意

\(m\)个人\(n\)个物品，第\(i\)个物品生命值为\(A_i\)，第\(i\)个人每秒可以减少一个物品\(B_i\)的生命值，给出一个\(m \times n\)的矩阵，如果\(i\)行\(j\)列为\(1\)，则表示第\(i\)个人可以攻击第\(j\)个物品，否则不能攻击，问至少需要多少秒，能干掉所有物品。一个物品被干掉当且仅当生命值小于等于\(0\)。（\(n, m \le 50, 1 \le A_i \le 10^5, 1 \le B_i \le 10^3\) ）

分析

我们可以二分时间，然后判定。首先我们能列出一些约束，然后发现可以用上下界可行流来做。

题解

二分时间\(t\)，令\(C(i, j)\)表示第\(i\)个人减少了\(j\)物品的血量，容易列出：

\[\sum_{j} C(i, j) \le t \times B_i
\]

\[\sum_{i} C(i, j) \ge A_j
\]

然后裸上上下界可行流即可。

（即使有\(2500\)个点....

#include <bits/stdc++.h>
using namespace std;
typedef long double lf;
const lf eps=1e-5, oo=1e150;
const int N=55, M=55, nN=3005, nM=nN*8;
inline int dcmp(lf a) {
	return a>-eps?a>=eps:-1;
}
int ihead[nN], cnt, A[N], B[M], ed[M][N], n, m;
struct E {
	int next, from, to;
	lf cap;
}e[nM];
void add(int x, int y, lf cap) {
	e[++cnt]=(E){ihead[x], x, y, cap}; ihead[x]=cnt;
	e[++cnt]=(E){ihead[y], y, x, 0  }; ihead[y]=cnt;
}
lf isap(int s, int t, int n) {
	static int gap[nN], cur[nN], d[nN], p[nN];
	memset(gap, 0, sizeof(int)*(n+1));
	memset(d,	0, sizeof(int)*(n+1));
	memcpy(cur, ihead, sizeof(int)*(n+1));
	gap[0]=n;
	lf ret=0, f;
	for(int i, x=s; d[s]<n; ) {
		for(i=cur[x]; i; i=e[i].next) if(dcmp(e[i].cap)>0 && d[x]==d[e[i].to]+1) break;
		if(i) {
			p[e[i].to]=cur[x]=i; x=e[i].to;
			if(x==t) {
				for(f=oo; x!=s; x=e[p[x]].from) f=min(f, e[p[x]].cap);
				for(x=t;  x!=s; x=e[p[x]].from) e[p[x]].cap-=f, e[p[x]^1].cap+=f;
				ret+=f;
			}
		}
		else {
			if(!--gap[d[x]]) break;
			d[x]=n;
			cur[x]=ihead[x];
			for(i=cur[x]; i; i=e[i].next) if(dcmp(e[i].cap)>0 && d[x]>d[e[i].to]+1) d[x]=d[e[i].to]+1;
			++gap[d[x]];
			if(x!=s) x=e[p[x]].from;
		}
	}
	return ret;
}
bool check(lf R) {
	int in=0, all=n*m+n+m+2+2, s=n*m+n+m+1, t=s+1, S=t+1, T=S+1;
	memset(ihead, 0, sizeof(int)*(all+1));
	cnt=1;
	for(int i=1; i<=m; ++i) {
		int idi=n*m+i;
		add(s, idi, R*B[i]);
		for(int j=1; j<=n; ++j) {
			int idj=n*m+m+j;
			int id=(i-1)*n+j;
			add(idi, id, oo);
			if(ed[i][j]) {
				add(id, idj, oo);
			}
		}
	}
	for(int i=1; i<=n; ++i) {
		int id=n*m+m+i;
		add(id, T, A[i]);
		in+=A[i];
	}
	add(S, t, in);
	add(t, s, oo);
	lf sum=isap(S, T, all);
	if(dcmp(sum-in)!=0) {
		return 0;
	}
	return 1;
}
int main() {
	scanf("%d%d", &n, &m);
	for(int i=1; i<=n; ++i) {
		scanf("%d", &A[i]);
	}
	for(int i=1; i<=m; ++i) {
		scanf("%d", &B[i]);
	}
	for(int i=1; i<=m; ++i) {
		for(int j=1; j<=n; ++j) {
			scanf("%d", &ed[i][j]);
		}
	}
	lf l=0, r=1e6;
	while(dcmp(r-l)>0) {
		lf mid=(l+r)/2;
		if(check(mid)) {
			r=mid;
		}
		else {
			l=mid;
		}
	}
	printf("%.9Lf\n", r);
	return 0;
}