01分数规划问题:
给你n个物品,a[i]表示第i个物品的收益,b[i]表示代价,x[i]表示选或不选,求一个最佳方案使得下式取值最大
通用解:
当然是从式子上考虑,定义函数F[cnt]如下:
其中cnt就是目标式
仔细观察上式,可以发现
①如果F[cnt]>0,那么很显然存在一个cnt',使得cnt'>cnt满足F[cnt']>=0,也就说当前不是最优解
②当cnt递增,F[cnt]一定随着递减
③出现F[cnt]<0没有意义,因为当cnt最优时F[]=0(将cnt的式子带入F[]式子中,你会发现可以抵消为0)
既然这样,很容易想到二分或者随机+迭代两种方法
其实都是没问题的,那就看具体问题吧
1690: [Usaco2007 Dec]奶牛的旅行
Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1032 Solved: 547
[ Submit][ Status][ Discuss]
Description
作为对奶牛们辛勤工作的回报,Farmer John决定带她们去附近的大城市玩一天。旅行的前夜,奶牛们在兴奋地讨论如何最好地享受这难得的闲暇。 很幸运地,奶牛们找到了一张详细的城市地图,上面标注了城市中所有L(2 <= L <= 1000)座标志性建筑物(建筑物按1..L顺次编号),以及连接这些建筑物的P(2 <= P <= 5000)条道路。按照计划,那天早上Farmer John会开车将奶牛们送到某个她们指定的建筑物旁边,等奶牛们完成她们的整个旅行并回到出发点后,将她们接回农场。由于大城市中总是寸土寸金,所有的道路都很窄,政府不得不把它们都设定为通行方向固定的单行道。 尽管参观那些标志性建筑物的确很有意思,但如果你认为奶牛们同样享受穿行于大城市的车流中的话,你就大错特错了。与参观景点相反,奶牛们把走路定义为无趣且令她们厌烦的活动。对于编号为i的标志性建筑物,奶牛们清楚地知道参观它能给自己带来的乐趣值F_i (1 <= F_i <= 1000)。相对于奶牛们在走路上花的时间,她们参观建筑物的耗时可以忽略不计。 奶牛们同样仔细地研究过城市中的道路。她们知道第i条道路两端的建筑物 L1_i和L2_i(道路方向为L1_i -> L2_i),以及她们从道路的一头走到另一头所需要的时间T_i(1 <= T_i <= 1000)。 为了最好地享受她们的休息日,奶牛们希望她们在一整天中平均每单位时间内获得的乐趣值最大。当然咯,奶牛们不会愿意把同一个建筑物参观两遍,也就是说,虽然她们可以两次经过同一个建筑物,但她们的乐趣值只会增加一次。顺便说一句,为了让奶牛们得到一些锻炼,Farmer John要求奶牛们参观至少2个建筑物。 请你写个程序,帮奶牛们计算一下她们能得到的最大平均乐趣值。
Input
* 第1行: 2个用空格隔开的整数:L 和 P
* 第2..L+1行: 第i+1行仅有1个整数:F_i * 第L+2..L+P+1行: 第L+i+1行用3个用空格隔开的整数:L1_i,L2_i以及T_i, 描述了第i条道路。
Output
* 第1行: 输出1个实数,保留到小数点后2位(直接输出,不要做任何特殊的取 整操作),表示如果奶牛按题目中描述的一系列规则来安排她们的旅 行的话,她们能获得的最大平均乐趣值
Sample Input
Sample Output
题意:
一张n个点m条边的有向图,边有花费,点有价值,点可以多次经过但价值不叠加,边花费叠加
求一个环满足收益和/花费和最大
如果知道01分数规划的话,可以一眼发现就是这种题
公式:
其中a[i]就是每个建筑物的观赏价值,b[i]就是每条单向边要花费的时间
二分cnt也就是二分答案,如果存在正权环,则说明当前不是最优解
把d[i]取反,就成简单的判负权环了,用spfa dfs即可
#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
using namespace std;
vector<int> G[1005];
int ok, n, val[1005], vis[1005], road[1005][1005];
double bet[1005];
void SPFA(int u, double cnt)
{
int i, v;
vis[u] = 1;
for(i=0;i<G[u].size();i++)
{
v = G[u][i];
if(cnt*road[u][v]-val[v]+bet[u]<bet[v])
{
if(vis[v]==1)
{
ok = 1;
return;
}
bet[v] = cnt*road[u][v]-val[v]+bet[u];
SPFA(v, cnt);
}
}
vis[u] = 0;
}
int Jud(double cnt)
{
int i;
ok = 0;
memset(bet, 0, sizeof(bet));
memset(vis, 0, sizeof(vis));
for(i=1;i<=n;i++)
{
SPFA(i, cnt);
if(ok)
return 1;
}
return 0;
}
int main(void)
{
int m, i, x, y, len;
double l, r, mid;
while(scanf("%d%d", &n, &m)!=EOF)
{
for(i=1;i<=n;i++)
{
scanf("%d", &val[i]);
G[i].clear();
}
memset(road, 62, sizeof(road));
for(i=1;i<=m;i++)
{
scanf("%d%d%d", &x, &y, &len);
G[x].push_back(y);
road[x][y] = min(road[x][y], len);
}
l = 0, r = 100000;
while(r-l>0.0001)
{
mid = (l+r)/2;
if(Jud(mid))
l = mid;
else
r = mid;
}
printf("%.2f\n", l);
}
return 0;
}