题目:
题解:
你二分一个答案
然后画一波柿子:
(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);
}