分数规划类问题选讲

0-1分数规划

模型

给定一些二元组$(x_i,y_i)$，让你选择其中的一些二元组，使$x_i$的和与$y_i$的和比值最大（小），即求
$$\max 或\min (\frac{\sum x_i}{\sum y_i})$$

通用解法

设$t=\frac{\sum x_i}{\sum y_i}$，那么我们有$t\ast \sum y_i=\sum x_i,即\sum x_i-t\ast \sum y_i=0$

当原式的值最大（小）时，$t$的值也最大（小），所以我们可以随便选取一个数$t$，然后去验证是否有一些二元组满足$\sum x_i-t\ast \sum y_i\ge 0$，如果有，那么该$t$则为一个合法值，然后找出所有合法值里的最大值就好了。

容易想到用二分实现该选数过程，那么如何验证我们选取的$t$值合法呢？

由$\sum x_i-t\ast \sum y_i\ge 0$得

$(x_1-t\ast y_1)+(x_2-t\ast y_2)+...+(x_n-t\ast y_n)\ge 0$

所以我们可以$O(n)$处理出每一个$x_i-t\ast y_i$，然后排序选择前$k$大的元素求和，如果和大于$0$，那么该值合法。

基本题型

0-1分数规划

POJ2976

题意简述：你有$n$场考试，每场考试有$b$道题，而你答对了$a$道题，现在你可以取消$k$场考试的成绩，问你的准确率最大是多少。

解法：和上面说的一样，但注意留下来的成绩是$n-k$场考试的成绩，二分是注意实数域的精度。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define eps 1e-6
using namespace std;
int n,k;
double a[100005],b[100005],s[100005];
bool check(double x){
	for(int i=1;i<=n;i++){
		s[i]=a[i]-b[i]*x;
	}
	sort(s+1,s+n+1);
	double ans=0;
	for(int i=n;i>=k+1;i--){
		ans+=s[i];
	}
	return ans>=0;
}
int main(){
	while(~scanf("%d%d",&n,&k)){
		if(!n&&!k)  break;
		for(int i=1;i<=n;i++){
			scanf("%lf",&a[i]);
		}
		for(int i=1;i<=n;i++){
			scanf("%lf",&b[i]);
		}
		double l=0.0,r=100.0;
		while(fabs(r-l)>eps){
			double mid=(l+r)/2.0;
			if(check(mid))  l=mid;
			else r=mid;
		}
		printf("%.0lf\n",l*100.0);
	}
}

最优比率生成树

POJ2728

题意简述：给你$n$个点，坐标为$(x,y)$，高度为$z$，你可以在任意两个点之间连边，边长为两点间直线距离，成本为两点间高度差。你需要用$n-1$条边将所有点和$1$号点连通并且这些边的总成本和总长度的比必须最小化。

模型：给你一个图，每条边有一个花费$p$和一个权值$s$，求出这张图的一棵生成树，使$\frac{\sum p}{\sum s}$最小。

观察到我们要求的这个式子和上文所讲的差不多，那么我们是否可以转化一下呢？

依然先二分一个值$t$，使$\sum p-\sum s\ast t=0$

用与上文相同的方法，将图中每条边的权值变为$p_i-s_i\ast t$，因为我们要求花费最小，所以我们直接求最小生成树，如果边权和大于$0$，那么不满足题意，返回$false$，反之则返回$true$

值得注意的一点是，由题目的要求，这张图的边数是接近$n^2$的，故此处用$Prim$算法求最小生成树最优，时间复杂度$O(n^2\ast log10^7)$

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define inf 1000000007
#define eps 1e-6
#pragma GCC optimize(3)
using namespace std;
int n;
double x[1005],y[1005],z[1005],cur[1005],vis[1005];
double mp[1005][1005],height[1005][1005],dis[1005][1005];
double calc(int i,int j){
	return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}
double prim(){
	double ans=0;
	for(int i=1;i<=n;i++){
		vis[i]=0;cur[i]=mp[1][i];
	}
	vis[1]=1;
	for(int i=1;i<n;i++){
		double minn=inf;
		int pos=0;
		for(int j=1;j<=n;j++){
			if(!vis[j]&&minn>cur[j])  pos=j,minn=cur[j];
		}
		if(minn==inf)  break;
		vis[pos]=1;
		ans+=minn;
		for(int j=1;j<=n;j++){
			if(!vis[j])  cur[j]=min(cur[j],mp[pos][j]);
		}
	}
	return ans;
}
bool check(double x){
	for(int i=1;i<=n;i++){
		for(int j=i+1;j<=n;j++){
			mp[i][j]=mp[j][i]=height[i][j]-x*dis[i][j];
		}
	}
	return  prim()<0;
}
int main(){
	while(1){
		scanf("%d",&n);
		if(!n)  break;
		for(int i=1;i<=n;i++){
			scanf("%lf%lf%lf",&x[i],&y[i],&z[i]);
		}
		for(int i=1;i<=n;i++){
			for(int j=i+1;j<=n;j++){
				if(i==j)  continue;
				dis[i][j]=dis[j][i]=calc(i,j);
				height[i][j]=height[j][i]=fabs(z[i]-z[j]);
			}
		}
		double l=0.0,r=1e7;
		while(fabs(r-l)>eps){
			double mid=(l+r)/2;
			if(check(mid))  r=mid;
			else  l=mid;
		}
		printf("%.3lf\n",l);
	}
}

最优比率生成环

天路

模型大体同最优比率生成树，只是我们需要找到一个环。

但我们有一些不同的地方：

求最大值时，我们设最大值为$t$，那么我们有$\sum s\ast t-\sum p=0$

求最小值时，我们设最小值为$t$，那么我们有$\sum p-\sum s\ast t=0$

求最大值按照最优比率生成树的方法处理，然后用$spfa$找负环，如果我们二分的值过大，那么就会出现负环，最小值时同理。

#include<bits/stdc++.h>
using namespace std;
int first[20005],nxt[20005],to[20005],tot=0;
double w[200005],val[200005];
void Add(int x,int y,double z,double d){
	nxt[++tot]=first[x];
	first[x]=tot;
	to[tot]=y;
	w[tot]=z;
	val[tot]=d;
}
double d[100005];
int vis[100005];
bool check(int u,double x){
	vis[u]=1;
	for(int e=first[u];e;e=nxt[e]){
		int v=to[e];
		if(d[v]>d[u]+val[e]*x-w[e]){
			if(vis[v])  return 1;
			d[v]=d[u]+val[e]*x-w[e];
			if(check(v,x))  return 1;
		}
	}
	vis[u]=0;
	return 0;
}
int main(){
	int n,m;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++){
		int x,y;
		double z,d;
		scanf("%d%d%lf%lf",&x,&y,&z,&d);
		Add(x,y,z,d);
	}
	double l=0,r=1e7;
	for(int i=1;i<=n;i++)  Add(0,i,0,0);
	while(fabs(r-l)>1e-6){
		double mid=(l+r)/2;
		memset(vis,0,sizeof(vis));
		memset(d,0x7f,sizeof(d));
		d[0]=0;
		if(check(0,mid))  l=mid;
		else  r=mid;
	}
	if(l==0)  cout<<"-1\n";  
	else printf("%.1lf\n",l);
}