bzoj 4819: [Sdoi2017]新生舞会 01分数规划 二分图最优匹配

Description

学校组织了一次新生舞会，Cathy作为经验丰富的老学姐，负责为同学们安排舞伴。有n个男生和n个女生参加舞会

买一个男生和一个女生一起跳舞，互为舞伴。Cathy收集了这些同学之间的关系，比如两个人之前认识没计算得出 

a[i][j] ，表示第i个男生和第j个女生一起跳舞时他们的喜悦程度。Cathy还需要考虑两个人一起跳舞是否方便，

比如身高体重差别会不会太大，计算得出 b[i][j]，表示第i个男生和第j个女生一起跳舞时的不协调程度。当然，

还需要考虑很多其他问题。Cathy想先用一个程序通过a[i][j]和b[i][j]求出一种方案，再手动对方案进行微调。C

athy找到你，希望你帮她写那个程序。一个方案中有n对舞伴，假设没对舞伴的喜悦程度分别是a'1,a'2,...,a'n，

假设每对舞伴的不协调程度分别是b'1,b'2,...,b'n。令

C=(a'1+a'2+...+a'n)/(b'1+b'2+...+b'n),Cathy希望C值最大。

Input

第一行一个整数n。

接下来n行，每行n个整数，第i行第j个数表示a[i][j]。

接下来n行，每行n个整数，第i行第j个数表示b[i][j]。

1<=n<=100,1<=a[i][j],b[i][j]<=10^4

Output

一行一个数，表示C的最大值。四舍五入保留6位小数，选手输出的小数需要与标准输出相等

01分数规划+最优匹配都是除了模板的第一题，感觉很多∑ xi/ ∑ yi的最优值问题都可以转换成这种方法来解决。

首先我们令答案k=∑ ai/ ∑ bi，移项得到∑ ai-k*∑ bi=0，反之如果∑ ai-k*∑ bi>=0，移项可以得到k<=∑ ai/ ∑ bi，那么如果我们构造e[i]=a[i]-k*b[i]，那么如果各个男性与舞伴配对之后最大的∑ei之和都<=0那么就说明了答案是小于等于当前构造的k的，所以我们就找到了二分的性质，对k进行二分，便可以得到最优的k了。

剩下的一个问题是如何求解各个男性与舞伴配对之后最大的∑ei之和，由于男性与女性性别之内没有边，只有性别之外右边，而新的e[i]就相当于男性对女性有一条长度为边，那么就是裸的二分图最优匹配了，跑km算法即可判断。

下附AC代码。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define maxn 205
#define yyxzhj 250
using namespace std;
const double eps=1e-7;
int n,tim;
int match[maxn],vislef[maxn],visrig[maxn];
double exlef[maxn],exrig[maxn],need[maxn];
double a[maxn][maxn],b[maxn][maxn],e[maxn][maxn];
double myabs(double now)
{
	return now>0.0 ? now : -now; 
}
int dfs(int now)
{
	vislef[now]=1;
	for(int nex=1;nex<=n;nex++)
	{
		if(visrig[nex]) continue;
		double gap=exlef[now]+exrig[nex]-e[now][nex];
		if(myabs(gap)<=eps)
		{
			visrig[nex]=1;
			if(!match[nex] || dfs(match[nex]))
			{
				match[nex]=now;
				return true;
			}
		}
		else
		{
			need[nex]=min(need[nex],gap);
		}
	}
	return false;
}
int check(double now)
{
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		e[i][j]=a[i][j]-now*b[i][j];
	}
	
	for(int i=1;i<=n;i++)	
	{
		match[i]=0,exrig[i]=0;
		exlef[i]=e[i][1];
		for(int j=2;j<=n;j++)
		exlef[i]=max(exlef[i],e[i][j]);
	}
	
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)	need[j]=123456789.0;
		
		while(yyxzhj)
		{
			for(int j=1;j<=n;j++)	vislef[j]=visrig[j]=0;
			
			if(dfs(i)) break;
		
			double minn=123456789;
			
			for(int j=1;j<=n;j++)
			if(!visrig[j])
			minn=min(minn,need[j]);
			
			for(int j=1;j<=n;j++)
			{
				if(vislef[j]) exlef[j]-=minn;
				if(visrig[j]) exrig[j]+=minn;
				else need[j]-=minn;
			}
		}
	}
	double ans=0;
	for(int i=1;i<=n;i++)
	ans+=e[match[i]][i];
	return ans>=-eps;
}
int main()
{
	double l=0.0,r=10000.0;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	{
		scanf("%lf",&a[i][j]);
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	scanf("%lf",&b[i][j]);
	while(myabs(r-l)>=eps)
	{
		double mid=(l+r)/2.0;
		if(check(mid)) l=mid;
		else r=mid;
	}
	printf("%.6lf\n",l);
}