MyException - 我的异常网
当前位置:我的异常网» VSTS » Codevs 1183 泥泞的路经

Codevs 1183 泥泞的路经

www.MyException.Cn  网友分享于:2013-10-22  浏览:0次
Codevs 1183 泥泞的道路

codevs1183 泥泞的道路

时间限制: 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
3
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

spfa+二分。对于每一个ans,代入原式spfa检查是否满足。跑最长路。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <queue>
#include <cstring>
using namespace std;
#define maxn 105
#define ljt(i,a,b) for(int i=(a);i<=(b);i++)

int P[maxn][maxn],T[maxn][maxn],n;
double W[maxn][maxn],dis[maxn];
bool inq[maxn];
int inv[maxn];
queue <int> que;

int read()
{
    int sum=0;
    char ch=getchar();
    while(!(ch>='0' && ch<='9')) ch=getchar();
    while(ch>='0' && ch<='9') sum=sum*10+ch-'0',ch=getchar();
    return sum;
}

bool spfa()
{
    memset(inv,0,sizeof(inv));
    memset(dis,-0x3f,sizeof(dis));
    memset(inq,0,sizeof(inq));
    que.push(1),dis[1]=0,inq[1]=1,inv[1]=1;
    while(!que.empty())
    {
    int now=que.front();
    que.pop();
    inq[now]=0;
    ljt(i,1,n)
    {
        if(P[now][i] && dis[i]<W[now][i]+dis[now])
        {
        dis[i]=W[now][i]+dis[now];
        if(!inq[i])
        {
            que.push(i);
            inq[i]=1;
            inv[i]++;
            if(inv[i]>n) return 1;
        }
        }
    }
    }
    if(dis[n]>0) return 1;
    return 0;
}

bool slove(double x)
{
    memset(W,0,sizeof(W));
    ljt(i,1,n) ljt(j,1,n) W[i][j]=P[i][j]-x*T[i][j];
    return spfa();
}

double binary()
{
    double mid,l=0,r=100000,ans;
    while(r-l>0.0001)
    {
    mid=(r+l)/2;
    if(slove(mid)) l=mid,ans=mid;
    else r=mid;
    }
    return ans;
}

int main()
{
    n=read();
    ljt(i,1,n) ljt(j,1,n) P[i][j]=read();
    ljt(i,1,n) ljt(j,1,n) T[i][j]=read();
    printf("%.3lf",binary());
    return 0;
}

文章评论

鲜为人知的编程真相
鲜为人知的编程真相
程序员最害怕的5件事 你中招了吗?
程序员最害怕的5件事 你中招了吗?
那些性感的让人尖叫的程序员
那些性感的让人尖叫的程序员
每天工作4小时的程序员
每天工作4小时的程序员
为啥Android手机总会越用越慢?
为啥Android手机总会越用越慢?
总结2014中国互联网十大段子
总结2014中国互联网十大段子
程序员的鄙视链
程序员的鄙视链
程序员眼里IE浏览器是什么样的
程序员眼里IE浏览器是什么样的
程序员周末都喜欢做什么?
程序员周末都喜欢做什么?
代码女神横空出世
代码女神横空出世
初级 vs 高级开发者 哪个性价比更高?
初级 vs 高级开发者 哪个性价比更高?
那些争议最大的编程观点
那些争议最大的编程观点
我的丈夫是个程序员
我的丈夫是个程序员
什么才是优秀的用户界面设计
什么才是优秀的用户界面设计
如何区分一个程序员是“老手“还是“新手“?
如何区分一个程序员是“老手“还是“新手“?
一个程序员的时间管理
一个程序员的时间管理
如何成为一名黑客
如何成为一名黑客
5款最佳正则表达式编辑调试器
5款最佳正则表达式编辑调试器
科技史上最臭名昭著的13大罪犯
科技史上最臭名昭著的13大罪犯
程序员都该阅读的书
程序员都该阅读的书
不懂技术不要对懂技术的人说这很容易实现
不懂技术不要对懂技术的人说这很容易实现
Java 与 .NET 的平台发展之争
Java 与 .NET 的平台发展之争
聊聊HTTPS和SSL/TLS协议
聊聊HTTPS和SSL/TLS协议
60个开发者不容错过的免费资源库
60个开发者不容错过的免费资源库
亲爱的项目经理,我恨你
亲爱的项目经理,我恨你
漫画:程序员的工作
漫画:程序员的工作
中美印日四国程序员比较
中美印日四国程序员比较
 程序员的样子
程序员的样子
2013年中国软件开发者薪资调查报告
2013年中国软件开发者薪资调查报告
10个帮程序员减压放松的网站
10个帮程序员减压放松的网站
我是如何打败拖延症的
我是如何打败拖延症的
程序员应该关注的一些事儿
程序员应该关注的一些事儿
十大编程算法助程序员走上高手之路
十大编程算法助程序员走上高手之路
Web开发者需具备的8个好习惯
Web开发者需具备的8个好习惯
编程语言是女人
编程语言是女人
程序猿的崛起——Growth Hacker
程序猿的崛起——Growth Hacker
10个调试和排错的小建议
10个调试和排错的小建议
为什么程序员都是夜猫子
为什么程序员都是夜猫子
老程序员的下场
老程序员的下场
2013年美国开发者薪资调查报告
2013年美国开发者薪资调查报告
团队中“技术大拿”并非越多越好
团队中“技术大拿”并非越多越好
看13位CEO、创始人和高管如何提高工作效率
看13位CEO、创始人和高管如何提高工作效率
做程序猿的老婆应该注意的一些事情
做程序猿的老婆应该注意的一些事情
程序员和编码员之间的区别
程序员和编码员之间的区别
写给自己也写给你 自己到底该何去何从
写给自己也写给你 自己到底该何去何从
Web开发人员为什么越来越懒了?
Web开发人员为什么越来越懒了?
“懒”出效率是程序员的美德
“懒”出效率是程序员的美德
程序员的一天:一寸光阴一寸金
程序员的一天:一寸光阴一寸金
要嫁就嫁程序猿—钱多话少死的早
要嫁就嫁程序猿—钱多话少死的早
软件开发程序错误异常ExceptionCopyright © 2009-2015 MyException 版权所有