bzoj 1690: [Usaco2007 Dec]奶牛的旅行

传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1690

思路:0-1分数规划

其实0-1分数规划重要的是它的思想,就像最小乘积XXX(生成树,匹配等)一样,知道了思想,再套一个对应的的算法即可

一个不错的博客:http://www.cnblogs.com/perseawe/archive/2012/05/03/01fsgh.html

说下我的理解:

首先0-1分数规划是用来解决这一类问题:

有N个东西,可以选或不选(用0,1来表示),选了有A[i]的收益,但要付出B[i]的代价

最后求一种方案使得(ΣA[i])/(ΣB[i])最大(或最小)

形式上和最小乘积XXX很像,就是乘号改成除号

但是解决方法却很不同


首先我们用x[i]取1或0表示i是选还是不选(如果x[i]的取值不限制在0,1中就是更一般的分数规划了)

先讨论最大的情况,最小类似

那么假设最优解为R=Σ(A[i]*x[i])/Σ(B[i]*x[i])

移项Σ(A[i]*x[i])-R*Σ(B[i]*x[i])=0

构造一个新函数f(L)=Σ(A[i]*x[i])-L*Σ(B[i]*x[i])

合并同类项令d[i]=Σ((A[i]-L*B[i])*x[i])

f(L)=Σ(d[i]*x[i])

如果有一个方案X使得f(L)>0,

那么Σ(A[i]*x[i])/Σ(B[i]*x[i])>L

这说明这个方案的值比L更大,那我们就选它

d数组是随L单调减的,所以任何一种方案,L减小,f(L)都会下降,所以这个函数就有了单调性

于是我们就可以二分了


这题就是求最优比例环

首先不会有环套环,这个随便证一下就可以了

每次二分L,把边权赋成d[i],把收益也转到边上,判断是否存在正权环即可


#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
const int maxn=1010,maxm=5010;
const double eps=1e-5;
using namespace std;
int n,m,v[maxn],pre[maxm],now[maxn],son[maxm],val[maxm],tot;double ans=0.0,mid,dis[maxn];bool vis[maxn];
void add(int a,int b,int c){pre[++tot]=now[a],now[a]=tot,son[tot]=b,val[tot]=c;}

int spfa(int x){
	vis[x]=1;
	for (int y=now[x];y;y=pre[y]){
		double nval=v[son[y]]-mid*val[y];
		if (dis[son[y]]<dis[x]+nval){
			if (!vis[son[y]]){
				dis[son[y]]=dis[x]+nval;
				if (spfa(son[y])) return 1;
			}
			else return 1;
		}
	}
	return vis[x]=0;
}

int main(){
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++) scanf("%d",&v[i]);
	for (int i=1,x,y,z;i<=m;i++) scanf("%d%d%d",&x,&y,&z),add(x,y,z);
	double l=0.0,r=1000.0;mid=(l+r)/2.0;
	while (r-l>eps){
		memset(dis,254,sizeof(dis));
		memset(vis,0,sizeof(vis));
		dis[1]=0.0;
		if (spfa(1)) ans=mid,l=mid;
		else r=mid;
		//printf("%.2f %.2f\n",l,r);
		mid=(l+r)/2.0;
	}
	printf("%.2f\n",ans);
	return 0;
}





posted @ 2016-01-19 08:33  orzpps  阅读(97)  评论(0编辑  收藏  举报