codevs1183 泥泞的道路(01分数规划)

1183 泥泞的道路

时间限制: 1 s

空间限制: 128000 KB

题目等级 : 钻石 Diamond

题目描述 Description

CS有n个小区，并且任意小区之间都有两条单向道路（a到b，b到a）相连。因为最近下了很多暴雨，很多道路都被淹了，不同的道路泥泞程度不同。小A经过对近期天气和地形的科学分析，绘出了每条道路能顺利通过的时间以及这条路的长度。

现在小A在小区1，他希望能够很顺利地到达目的地小区n，请帮助小明找出一条从小区1出发到达小区n的所有路线中（总路程/总时间）最大的路线。请你告诉他这个值。

输入描述 Input Description

第一行包含一个整数n，为小区数。

接下来n*n的矩阵P，其中第i行第j个数表示从小区i到小区j的道路长度为Pi,j。第i行第i个数的元素为0，其余保证为正整数。

接下来n*n的矩阵T，第i行第j个数表示从小区i到小区j需要的时间Ti,j。第i行第i个数的元素为0，其余保证为正整数。

输出描述 Output Description

写入一个实数S，为小区1到达n的最大答案，S精确到小数点后3位。

样例输入 Sample Input

0 8 7

9 0 10

5 7 0

0 7 6

6 0 6

6 2 0

样例输出 Sample Output

2.125

数据范围及提示 Data Size & Hint

【数据说明】

30%的数据，n<=20

100%的数据，n<=100，p,t<=10000

/*
01分数规划
总路程/总时间=平均速度
可以二分v,建图spfa最长路 判断S总-v*T总=0为最优解,>0当前v满足。
若有正权回路则当前v一定满足(求最长路可以判正环)。
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>#define N 1001
#define inf 0x3f3f3f3fusing namespace std;
int n,m,ans,cnt,flag;
int s[N][N],t[N][N],in[N];
double dis[N][N],d[N];
bool inq[N];
queue<int>q;void clear()
{memset(inq,0,sizeof inq);while(!q.empty()) q.pop();memset(in,0,sizeof in);    for(int i=1;i<=n;i++) d[i]=-inf;
}bool spfa(int u)
{clear();q.push(u);inq[u]=1;d[u]=0;in[u]=0;while(!q.empty()){int now=q.front();q.pop();inq[now]=0; in[now]++;if(in[now]>=n) return true;for(int i=1;i<=n;i++){if(i==now) continue;if(d[i]<d[now]+dis[now][i] && s[now][i]){d[i]=d[now]+dis[now][i];if(!inq[i]) q.push(i),inq[i]=1;} }}return d[n]>0;
}bool check(double x)
{for(int i=1;i<=n;i++) for(int j=1;j<=n;j++)dis[i][j]=s[i][j]-t[i][j]*x;return spfa(1);
}int main()
{scanf("%d",&n);for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) scanf("%d",&s[i][j]);for(int i=1;i<=n;i++) for(int j=1;j<=n;j++)scanf("%d",&t[i][j]);double eps=0.00001,l=0,r=10000;while(r-l>eps){double mid=(l+r)/2;if(check(mid)) l=mid;else r=mid;}printf("%.3lf\n",l);
}

转载于:https://www.cnblogs.com/L-Memory/p/9748383.html

codevs1183 泥泞的道路(01分数规划)相关推荐

bzoj 1690: [Usaco2007 Dec]奶牛的旅行（01分数规划--最优比率环）
01分数规划问题: 给你n个物品,a[i]表示第i个物品的收益,b[i]表示代价,x[i]表示选或不选,求一个最佳方案使得下式取值最大通用解: 当然是从式子上考虑,定义函数F[cnt]如下: 其中c ...
BZOJ1690 Usaco2007 Dec 奶牛的旅行【01分数规划】
BZOJ1690 Usaco2007 Dec 奶牛的旅行题目描述作为对奶牛们辛勤工作的回报,Farmer John决定带她们去附近的大城市玩一天.旅行的前夜,奶牛们在兴奋地讨论如何最好地享受这难得 ...
poj2976（0-1分数规划）
0-1分数规划设x[i]等于1或0. 则我们所求的比率 rate = ∑(cost[i] * x[i]) / ∑(cost[i] * x[i]), 0≤i<m . z( rate ) = ∑( ...
点分治问题 ----------- luoguP2942 [WC2010]重建计划 [点分治 + bfs + 单调队列 + 预处理建树 + 二分 + 01分数规划]
题目链接解题思路: 1.对于这个Avgvalue=∑e∈sv(e)∣s∣Avgvalue = \frac{\sum_{e\in s}v(e)}{|s|}Avgvalue=∣s∣∑e∈sv(e) ...
解题报告：AcWing 1165. 单词环（01分数规划、hash、经验优化）
本题的关键在于: 建图 01分数规划本题的数据过大,如果直接spfa判断会TLE,因此我们使用经验优化,就是如果所有的点入队的次数过多,比如大于100000,那么我们直接认为它是存在正环的.(免去T ...
模板 - 判断负环（超时高效优化技巧）、01分数规划
整理的算法模板合集: ACM模板判断负环判正环求最长路,判负环求最短路 int n; // 总点数 int h[N], w[N], e[N], ne[N], idx; // 邻接表存储所有边 in ...
【图论专题】负环与01分数规划
整理的算法模板合集: ACM模板题目列表: 题目算法 AcWing 904. 虫洞 spfa判负环 AcWing 361. 观光奶牛最优比率环.01分数规划 AcWing 1165. 单词环 0 ...
POJ2728 Desert King ——01分数规划Dinkelbach迭代法+最小生成树prim算法
首先,纪念我用Linux系统AC的第一题- 安装这个万恶的NOI Linux系统费了6小时的时间,不过好在最后终于装上了,但是因为我安装的Linux系统比较烂,还遭到了小花儿和js的鄙视,唉,本人 ...
POJ 2728 01分数规划
题意: 最优比率生成树,要求生成树中的所有边的花费与所有边的长度的比值最小题解: 01分数规划,详见http://www.cnblogs.com/proverbs/archive/2013/01/0 ...

codevs1183 泥泞的道路(01分数规划)

1183 泥泞的道路

codevs1183 泥泞的道路(01分数规划)相关推荐

最新文章

热门文章