SPFA判负环板题

Link with Game Glitch

分析：

• SPFA判负环板题

• 转换一下题意，就是要找最大的系数w，使得无法成环

  无法成环：环上所有边权的乘积小于1

  二分答案，然后跑spfa判负环即可（注意会卡double的精度）

• 一点小优化，原本是系数乘，可能会很大，可以对所有系数取对数，即转换成求和，但本题并没有卡

#include <bits/stdc++.h>
#define int long long 
#define Pa pair<int,int> 
#define yes cout<<"YES\n"
#define no cout<<"NO\n"
#define double long double
using namespace std;    

const int N=1e3+5;
struct node
{
    int v; double w;
};
vector <node> g[N];
int vis[N], dep[N];
int n,m; 
double d[N];
bool ck(double k)
{
    queue <int> q;
    for(int i=1;i<=n;i++)
    {
        q.push(i); vis[i]=d[i]=dep[i]=1;
    }
    while(!q.empty())
    {
        int u=q.front(); q.pop();
        vis[u]=0;
        for(auto it : g[u])
        {
            int v=it.v ; double w=it.w;
            if(d[v]<d[u]*w*k)
            {
                d[v]=d[u]*w*k;
                dep[v]=dep[u]+1;
                if(dep[v]>n) return 1;
                if(!vis[v])
                {
                    q.push(v);
                    vis[v]=1;
                }
            }
        }
    }
    return 0;
}

void solve()
{
    cin>>n>>m;
    for(int i=1;i<=m;i++)
    {
        int a,b,c,d; 
        cin>>a>>b>>c>>d;
        g[b].push_back({d,1.0*c/a});
    }
    double l=0,r=1, eps=1e-8, ans;
    while(r-l>eps)
    {
        double mid=(l+r)/2;
        if(ck(mid))
        {
            r=mid;
        }
        else 
        {
            l=mid;
            ans=mid;
        }
    }
    cout<<ans<<"\n";
}
signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    cout<<fixed<<setprecision(12);
    int T=1;
    //cin>>T;
    while(T--) solve();
}