判断负环例题

acwing 904虫洞

#include<bits/stdc++.h>
using namespace std;
#define N 510
#define M 6010
int ne[M],to[M],h[N],w[M];
int m,n,t,cnt;
int rem[N];//the number of the point passed by
void add( int a,int b,int c){
    w[cnt]=c;to[cnt]=b;ne[cnt]=h[a];h[a]=cnt++;
}
queue<int>q;
int dist[N];
int vis[N];
int spfa(){
    while(!q.empty()) q.pop();
    memset(dist,0,sizeof(dist));
    memset(rem,0,sizeof(rem));
    for( int i=1;i<=n;i++){
        q.push(i);
        vis[i]=1;
    }
    while(!q.empty()){
        int now=q.front();
        q.pop();
        vis[now]=0;
        for( int i=h[now];i!=-1;i=ne[i]){
            int j=to[i];
            if(dist[j]>dist[now]+w[i]){
                dist[j]=dist[now]+w[i];
                rem[j]=rem[now]+1;//注意步数统计要加在这里
                if(rem[j]>=n) return 1;
                if(!vis[j]){
                    vis[j]=1;
                    q.push(j);
                }
            }
        }
    }
    return 0;
    
}
int main(){
    int t;
    cin>>t;
    while(t--){
        memset(h,-1,sizeof(h));
        int w;
        cnt=0;
        cin>>n>>m>>w;
        while(m--){
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            add(a,b,c);
            add(b,a,c);
        }
        while(w--){
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            add(a,b,-c);
        }
        if(spfa()) cout<<"YES"<<endl;
        else cout<<"NO"<<endl;
    }
    
}

acwing 361 观光奶牛
01分数规划+求负环
01分数规划主要解决两个和之比的问题
01分数规划一般使用二分来计算。
本题二分的设计思路是：
方法1说明spfa也可以求最大路径


方法1代码

#include<bits/stdc++.h>
using namespace std;
#define N 1010
#define M 10010
#define eps 1e-4
int n,m;
int w[M],to[M],ne[M],h[N];
int cnt;
int f[N];
void add(int a,int b,int c){
    w[cnt]=c;to[cnt]=b;ne[cnt]=h[a];h[a]=cnt++;
}
queue<int>q;
double dist[N];
int vis[N];
int rem[N];
int spfa(double bd){
    while(!q.empty())q.pop();
    memset(dist,0,sizeof(dist));
    memset(rem,0,sizeof(rem));
    for( int i=1;i<=n;i++){
        q.push(i);
        vis[i]=1;
    }
    while(!q.empty()){
        int now=q.front();
        q.pop();
        vis[now]=0;
        for( int i=h[now];i!=-1;i=ne[i]){
            int j=to[i];
            if(dist[j]<dist[now]-bd*w[i]+f[now]){
                dist[j]=dist[now]-bd*w[i]+f[now];
                rem[j]=rem[now]+1;
                if(rem[j]>=n+1) return 1;
                if(!vis[j]){
                    vis[j]=1;
                    q.push(j);
                }
            }
        }
    }
    return 0;
}
int main(){
    cin>>n>>m;
    memset(h,-1,sizeof (h));
    for( int i=1;i<=n;i++){
        cin>>f[i];
    }
    for( int i=1;i<=m;i++){
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
    }
    double l=0,r=1001;
    while(r-l>1e-4){
        double mid=(l+r)/2;
        if(spfa(mid)<=0){
            r=mid;
        }
        else{
            l=mid;
        }
    }
	printf("%.2lf",l);
}