SPFA找负环----------观光奶牛

给定一张L个点、P条边的有向图,每个点都有一个权值f[i],每条边都有一个权值t[i]。
求图中的一个环,使“环上各点的权值之和”除以“环上各边的权值之和”最大。
输出这个最大值。
注意:数据保证至少存在一个环。
输入格式
第一行包含两个整数L和P。
接下来L行每行一个整数,表示f[i]。
再接下来P行,每行三个整数a,b,t[i],表示点a和b之间存在一条边,边的权值为t[i]。
输出格式
输出一个数表示结果,保留两位小数。
数据范围
2≤L≤10002≤L≤1000,

2≤P≤50002≤P≤5000,

1≤f[i],t[i]≤10001≤f[i],t[i]≤1000
输入样例:
5 7
30
10
10
5
10
1 2 3
2 3 2
3 4 5
3 5 2
4 5 5
5 1 3
5 2 2

输出样例:
6.00

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 1010, M = 5010;
int n, m;
int wf[N];
int h[N], e[M], wl[M], ne[M], idx;
double dist[N];
int q[N], cnt[N];
bool st[N];
void add(int a, int b, int c){
 e[idx] = b, wl[idx] = c, ne[idx] = h[a], h[a] = idx ++;
}
bool check(double mid){
 memset(dist, 0, sizeof dist);
 memset(cnt, 0, sizeof cnt);
 memset(st, 0, sizeof st);
  int hh = 0, tt = 0;
 for (int i = 1; i <= n; i ++){
  q[tt ++] = i;
  st[i] = true;
 }
  while(hh != tt){
  int t = q[hh ++];
  if (hh == N)   hh = 0;
  st[t] = false;
   for (int i = h[t]; ~i; i = ne[i]){
   int j = e[i];
   if (dist[j] < dist[t] + wf[t] - mid * wl[i]){
       dist[j] = dist[t] + wf[t] - mid * wl[i];
   cnt[j] = cnt[t] + 1;
   if (cnt[j] >= n)   return true;
   if (!st[j]){
    q[tt ++] = j;
    if (tt == N)   tt = 0;
    st[j] = true;
   }
  }
 }
 }
 return false;
}
int main(){
 cin >> n >> m;
 for (int i = 1; i <= n; i ++)   scanf("%d", &wf[i]);
  memset(h, -1, sizeof h);
 for (int j = 0; j < m; j ++){
  int a, b, c;
  cin >> a >> b >> c;
  add(a, b, c);
 }
  double l = 0, r = 1010;
 while(r - l > 1e-4){
  double mid = (l + r) / 2;
  if (check(mid))   l = mid;
  else              r = mid;
 }
  printf("%.2lf\n", l);
  return 0;
}

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值