SDOI2015 星际战争

最新推荐文章于 2020-08-17 09:48:06 发布

原创最新推荐文章于 2020-08-17 09:48:06 发布 · 139 阅读

CC 4.0 BY-SA版权

本文深入探讨了最小费用流算法在网络流问题中的应用，详细解释了如何通过调整网络中的费用来寻找最优解，同时提供了具体的代码实现，帮助读者理解算法的工作原理。

传送门
?_?不要到 $d a r k b z o j$ 上交。。
几个 $s b$ 错误：
①分层的时候没有 $Q . p u s h (S)$ 。
② $c h e c k$ 的时候 $i n i t$ 放在了 $d i n i c$ 后面。。
③数组没开够。

#include<bits/stdc++.h>
#define re register
#define cs const
using namespace std;
cs int N=55,M=5300;
cs double oo=1e5,eps=1e-8;
double a[N],b[N],sumA=0;
int n,m,stB,edB,S,T,attack;
inline int sgn(double x){
	if(x>eps) return 1;
	if(x<-eps) return -1;
	return 0;
}
namespace NETWORK{
	int Head[N<<1],Next[M],V[M],cnt=1;
	int dep[N<<1],cur[N<<1];double W[M],copy_W[M];
	inline void add(int u,int v,double w){
		Next[++cnt]=Head[u],V[cnt]=v,W[cnt]=w,copy_W[cnt]=w,Head[u]=cnt;
	}
	inline bool bfs(){
		memset(dep,0,sizeof dep);
		memcpy(cur,Head,sizeof Head);
		queue<int> Q;dep[S]=1,Q.push(S);
		while(!Q.empty()){
			int u=Q.front();Q.pop();
			for(int re i=Head[u],v=V[i];i;v=V[i=Next[i]])
				if(W[i]>0&&(!dep[v])) dep[v]=dep[u]+1,Q.push(v);
		}return dep[T];
	}
	inline double dfs(int u,double dist){
		if(u==T) return dist;
		for(int re &i=cur[u],v=V[i];i;v=V[i=Next[i]]){
			if(W[i]>0&& dep[v]==dep[u]+1){
				double d=dfs(v,min(dist,W[i]));
				if(d>0){W[i]-=d,W[i^1]+=d;return d;}
			}
		}return 0.0;
	}
	inline double dinic(double ret=0){
		while(bfs()){
			double d;
			while((d=dfs(S,oo))>0) ret+=d;
		}return ret;
	}
	inline void init(double x){
		for(int re i=stB;i<=edB;i+=2)
			W[i]=x*b[i/2],W[i^1]=0;
		for(int re i=edB+2;i<=cnt;++i) W[i]=copy_W[i];
	}
	inline bool check(double x){
		return init(x),sgn(dinic()-sumA)==0;
	}
}
using NETWORK::add;
using NETWORK::check;
int main(){
//	freopen("2706.in","r",stdin);
	scanf("%d%d",&n,&m),S=m+n+1,T=S+1;
	for(int i=1;i<=n;++i) scanf("%lf",&a[i]),sumA+=a[i];
	for(int i=1;i<=m;++i) scanf("%lf",&b[i]);
	stB=2,edB=m<<1;
	for(int re i=1;i<=m;++i) add(S,i,0),add(i,S,0);
	for(int re i=1;i<=m;++i){
		for(int re j=1;j<=n;++j){
			scanf("%d",&attack);
			if(attack) add(i,m+j,oo),add(m+j,i,0);
		}
	}
	for(int re i=1;i<=n;++i)
		add(m+i,T,a[i]),add(T,m+i,0);
	
	double l=0,r=oo;
	while(r-l>eps){
		double mid=(l+r)/2.0;
		if(check(mid)) r=mid;
		else l=mid;
	}printf("%.6lf\n",l);
}