bzoj 3993: [SDOI2015]星际战争 二分答案&网络流

       为了学ISAP做了一下这道题。感觉ISAP并没有比DINIC难打多少。。。

       二分答案后源点向激光连边，权为伤害总量；激光向可以达到的机器人连无穷大；机器人向汇点连边，权为装甲值。然后汇点的边满流就满足答案。

       但是奇怪的是把二分的上界设成5e6就能过，改成sum（所有机器人装甲值的和）就过不了。。然后加上一句 if (l>sum)while(1); 之后还是过了。。说明答案不会超过sum啊。。。希望老司机能带带我。

AC代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 105
#define M 5505
#define inf 1000000000
using namespace std;

int n,m,tot=1,gol,fst[N],cur[N],pnt[M],nxt[M],d[N],num[N],h[N],fa[N];
double sum,a[N],b[N],len[M];
void add(int x,int y,double z){
	pnt[++tot]=y; len[tot]=z; nxt[tot]=fst[x]; fst[x]=tot;
}
void bfs(){
	memset(num,0,sizeof(num));
	memset(d,0x3f,sizeof(d)); d[gol]=0;
	int head=0,tail=1; h[1]=gol;
	while (head<tail){
		int x=h[++head],p; num[d[x]]++;
		for (p=fst[x]; p; p=nxt[p]) if (len[p^1]){
			int y=pnt[p];
			if (d[y]>inf){
				d[y]=d[x]+1; h[++tail]=y;
			}
		}
	}
}
double up(){
	int i; double tmp=inf;
	for (i=gol; i; i=pnt[fa[i]^1]) tmp=min(tmp,len[fa[i]]);
	for (i=gol; i; i=pnt[fa[i]^1]){
		len[fa[i]]-=tmp; len[fa[i]^1]+=tmp;
	}
	return tmp;
}
bool ok(double t){
	int p,i,x=0;
	for (i=2; i<=tot; i+=2){
		len[i]+=len[i^1]; len[i^1]=0;
	}
	for (p=fst[0]; p; p=nxt[p]) len[p]=b[pnt[p]]*t;
	memcpy(cur,fst,sizeof(fst));
	double flow=0; bfs(); 
	while (d[0]<=gol){
		if (x==gol){ flow+=up(); x=0; }
		bool flag=1;
		for (p=cur[x]; p; p=nxt[p]) if (len[p]){
			int y=pnt[p];
			if (d[x]==d[y]+1){
				cur[x]=fa[y]=p; x=y;
				flag=0; break;
			}
		}
		if (flag){
			int mn=gol; cur[x]=fst[x];
			for (p=fst[x]; p; p=nxt[p])
				if (len[p]) mn=min(mn,d[pnt[p]]);
			num[d[x]]--; if (!num[d[x]]) break;
			d[x]=mn+1; num[d[x]]++;
			if (x) x=pnt[fa[x]^1];
		}
	}
	return flow>=sum;
}
int main(){
	scanf("%d%d",&n,&m); gol=n+m+1; int i,j,x;
	for (i=1; i<=n; i++){
		scanf("%lf",&a[i]); sum+=a[i];
		add(i+m,gol,a[i]); add(gol,i+m,0);
	}
	for (i=1; i<=m; i++){
		scanf("%lf",&b[i]); add(0,i,b[i]); add(i,0,0);
	}
	for (i=1; i<=m; i++)
		for (j=1; j<=n; j++){
			scanf("%d",&x);
			if (x){ add(i,j+m,inf); add(j+m,i,0); }
		}
	double l=0,r=5e6,mid;
	while (l+1e-4<r){
		mid=(l+r)/2;
		if (ok(mid)) r=mid; else l=mid;
	}
	printf("%.6f\n",l);
	return 0;
}

by lych

2016.4.9