[codevs1183]泥泞的道路（二分+正环）

题目：

我是超链接

题解：

你二分一个答案
然后画一波柿子：
(p[1]+p[i1]+p[i2]+…+p[n])/(t[1]+t[i1]+t[i2]+…+t[n]) < mid 【哇mid好大的所有的路都够不到
(p[1]+p[i1]+p[i2]+…+p[n]) < mid*(t[1]+t[i1]+t[i2]+…+t[n])
p[1]-mid*t[1]+p[i1]-mid*t[i1]+…+p[n]-mid*t[n] < 0———r=mid 【变小变小
p[1]-mid*t[1]+p[i1]-mid*t[i1]+…+p[n]-mid*t[n] >= 0———l=mid 【反过来
这是。。。找了一波最长路，如果有任意一条路>=0就可以return 1，但是最长路如果出现正数怎么停得下来呢？判断正环？
SPFA的话就是判断进入队列（修改长度）的次数，修改超过n次就有正环

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
using namespace std;
const double eps=1e-5;
queue<int>q;
int p[105][105],t[105][105],n,tim[105];
double dis[105];bool vis[105];
bool check(double mid)
{
    memset(tim,0,sizeof(tim));  
    for (int i=1;i<=n;i++) dis[i]=-1e16;
    q.push(1);dis[1]=0;tim[1]++;
    while (!q.empty())
    {
        int now=q.front(); q.pop(); vis[now]=0;
        for (int i=1;i<=n;i++)
          if (i!=now && dis[i] < p[now][i]-mid*t[now][i]+dis[now])
          {
            dis[i]=(double)(p[now][i]-(double)mid*t[now][i])+dis[now];
            if (!vis[i]) q.push(i),vis[i]=1,tim[i]++;
            if (tim[i]>n) return 1;//如果有正环的话，最多循环n个点又到了i，所以如果ta被修改超过n次，相当于又被自己修改了，退出
          }
    }
    return dis[n]>=0;
}
int main()
{
    int i,j;
    scanf("%d",&n);
    for (i=1;i<=n;i++)
      for (j=1;j<=n;j++) scanf("%d",&p[i][j]);
    for (i=1;i<=n;i++)
      for (j=1;j<=n;j++) scanf("%d",&t[i][j]);
    double l=0,r=10000;
    while (r-l>=eps)
    {
        double mid=(l+r)/2.0;
        if (check(mid)) l=mid;
        else r=mid;
    }
    printf("%.3lf",l);
}