洛谷P1768 天路

题目描述

“那是一条神奇的天路诶~，把第一个神犇送上天堂~”，XDM先生唱着这首“亲切”的歌曲，一道猥琐题目的灵感在脑中出现了。

和C_SUNSHINE大神商量后，这道猥琐的题目终于出现在本次试题上了，旨在难到一帮大脑不够灵活的OIer们（JOHNKRAM真的不是说你……）。

言归正传，小X的梦中，他在西藏开了一家大型旅游公司，现在，他要为西藏的各个景点设计一组铁路线。但是，小X发现，来旅游的游客都很挑剔，他们乘火车在各个景点间游览，景点的趣味当然是不用说啦，关键是路上。试想，若是乘火车一圈转悠，却发现回到了游玩过的某个景点，花了一大堆钱却在路上看不到好的风景，那是有多么的恼火啊。

所以，小X为所有的路径定义了两个值，Vi和Pi，分别表示火车线路的风景趣味度和乘坐一次的价格。现在小X想知道，乘客从任意一个景点开始坐火车走过的一条回路上所有的V之和与P之和的比值的最大值。以便为顾客们推荐一条环绕旅游路线（路线不一定包含所有的景点，但是不可以存在重复的火车路线）。

于是，小X梦醒之后找到了你……

输入输出格式

输入格式：

第一行两个正整数N，M，表示有N个景点，M条火车路线，火车路线是单向的。

以下M行，每行4个正整数，分别表示一条路线的起点，终点，V值和P值。

注意，两个顶点间可能有多条轨道，但一次只能走其中的一条。

输出格式：

一个实数，表示一条回路上最大的比值，保留1位小数。

若没有回路，输出-1。

输入输出样例

输入样例#1：  复制

5 6
1 2 1 1
4 1 6 2
5 4 8 1
2 3 2 2
5 2 4 1
3 5 6 4

输出样例#1：  复制

2.3

说明

对于30%的数据，1≤N≤100,1≤M≤20;

对于60%的数据，1≤N≤3,000,1≤M≤2,000;

对于100%的数据，1≤N≤7,000,1≤M≤20,000,1≤Vi,Pi≤1,000.

保证答案在200以内.

这一题直接找环肯定超时，甚至60分都不好拿。我们这样想：设ans为答案，则存在任意环E保证∑V(i)/∑P(i)<=ans，可以得到∑V(i)<=ans*∑P(i)，于是∑{P(i)*ans}-∑V(i)>=0。

最后得到∑{P(i)*ans-V(i)}>=0，即以P(i)*ans-V(i)为边权的图中不存在负权环。（我没想到这点导致我不会做。。。）

一道dpfa应用的好题

代码：

#include<bits/stdc++.h>
using namespace std;
int h[10001],next[130001],v[130001],n,m;
float p[130001],w[130001],d[130001],low,high,mid;
bool in[130001],flag;
void dfs(int now)
{
    in[now]=1;
    for(int i=h[now];i;i=next[i])
      if(d[now]+p[i]*mid-w[i]<=d[v[i]])
      {
        if(in[v[i]]){flag=0;return;}
        d[v[i]]=d[now]+p[i]*mid-w[i];
        dfs(v[i]);
        if(!flag)return;
      }
    in[now]=0;
}
void init()
{
    scanf("%d%d",&n,&m);
    int i,j,k;
    for(j=1;j<=m;j++)
    {
        scanf("%d%d%f%f",&i,v+j,w+j,p+j);
        if(i!=v[m])
        {
          next[j]=h[i];
          h[i]=j;
        }
    }
    for(i=1;i<=n;i++)
    {
        v[i+m]=i;
        next[i+m]=h[0];
        h[0]=i+m;
    }
}
int main()
{
    init();
    low=0;
    high=200;
    while(low+0.01<high)
    {
        mid=(low+high)/2;
        flag=1;
        memset(in,0,sizeof(in));
        memset(d,127,sizeof(d));
        d[0]=0;
        dfs(0);
        if(!flag)low=mid;
        else high=mid;
    }
    if(low==0)printf("-1");
    else printf("%.1f",high);
    return 0;
}