bzoj4819 [Sdoi2017]新生舞会 分数规划（实数二分）+网络流检验

这个题其实就是二分检验。。

稍微移个项就可以得到（a1-c*b1）+（a2-c*b2）+……

这样贡献就化成了具体到每次决策的影响，c显然满足连续性，所以可以二分c

跑最大流以求得到最大的值使当前答案显得小，就可以得出答案

主要就是化动为静，，如果答案反过来影响过程就可以考虑枚举答案分析，然后来判断当前的影响值 和 作为限度的影响值

注：

1、l和r的eps

2、spfa的返回值

3、double的建边

4、前继一定是-1；

码：

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
queue<int>q;
#define N 50000
int tot=-1,hou[N],xia[N],zhong[N],v[N],s,t,i,j,n,m,qj[N];
bool vis[N];
double C[N],d[N],l,r,a[105][105],b[105][105],daan,ans;
void jian(int a,int b,double c,int d)
{
	++tot;hou[tot]=xia[a],xia[a]=tot,zhong[tot]=b,C[tot]=c,v[tot]=d;	
}
void jia(int a,int b,double v)
{
jian(a,b,v,1);
jian(b,a,-v,0);	
}
bool spfa()
{
	int i;
	for(i=1;i<=t;i++)
	{
	d[i]=-99999999;	
	}
	memset(qj,-1,sizeof(qj));
q.push(s);	
d[s]=0;
	while(!q.empty())
	{
	int st=q.front();
	q.pop();
	vis[st]=0;
	for(i=xia[st];i!=-1;i=hou[i])
	{
		int nd=zhong[i];
		if(!v[i]||d[nd]>=d[st]+C[i])continue;
		qj[nd]=i;
		d[nd]=d[st]+C[i];
		if(vis[nd]==0)
		{
			vis[nd]=1;
			q.push(nd);
		}						
	}		
	}
	if(d[t]==-99999999)return 0;
	int o=qj[t],minf=999;
	while(o!=-1)
	{
		minf=min(minf,v[o]);
		o=qj[zhong[o^1]];
	}
		o=qj[t];
		while(o!=-1)
	{
		ans+=C[o]*(1.00000000*v[o]);
		v[o]-=minf;
		v[o^1]+=minf;
		o=qj[zhong[o^1]];
	}
//	cout<<ans<<" ";
	return 1;	
}
void mcmf()
{
while(spfa());
}
int main()
{
	scanf("%d",&n);
	s=n*n+1;
	t=n*n+2;
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=n;j++)
		{
	scanf("%lf",&a[i][j]);
	    }	
	}
		for(i=1;i<=n;i++)
	{
		for(j=1;j<=n;j++)
		{
	scanf("%lf",&b[i][j]);
	    }	
	}
	l=0.0000001;r=1000000.0000001;
	while(l+0.00000001<r)
	{
	double mid=(l+r)/2;
memset(xia,-1,sizeof(xia));
			tot=-1;
	for(i=1;i<=n;i++)
	{jia(s,i,0);
	for(j=1;j<=n;j++)
	{
		jia(i,n+j,a[i][j]-mid*b[i][j]);
	}jia(i+n,t,0);
    }ans=0;
	mcmf();
	//printf("%.8lf  %.8lf\n",mid,ans);
	if(ans<=0){r=mid;daan=mid;}
	else l=mid+0.0000001;	
	}
	printf("%.6lf",daan);
}