20140710 loop-I SPFA 二分答案

考试的时候写的dfs找环然后求平均值

当时感觉复杂度不是太高

结果Wa和T了一些后只有40分。。。

正解是二分答案再将所有边减去答案后用SPFA找有无负环即可

#include <cstdio>
#include <queue>
#include <vector>
using namespace std;
#define N 660
#define INF 0x3f3f3f+0.0

struct edge {
    int u,v;
};
struct edge1 {
    int u;
    double len;
};

int n,m;
int map[N][N];
vector <edge1> G[N];
bool vis[N];
double dis[N];
int outqueue[N];
edge e[1010];

bool spfa(double k) {
    for (int i=1;i<=n;i++) {
        dis[i]=INF;
        vis[i]=false;
        outqueue[i]=0;
    }
    vis[1]=true;
    dis[1]=0;
    queue <int> q;
    q.push(1);
    while (!q.empty()) {
        int v=q.front();
        q.pop(); vis[v]=false;
        outqueue[v]++;
        if (outqueue[v]>n) return false;
        for (int i=0;i<G[v].size();i++) {
            edge1 e=G[v][i];
            if (dis[v]+e.len-k<dis[e.u]) {
                dis[e.u]=dis[v]+e.len-k;
                if (!vis[e.u]) {
                    vis[e.u]=true;
                    q.push(e.u);
                }
            }
        }
    }
    return true;
}

double min(double a,double b) {
    return a<b?a:b;
}

int main() {
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) map[i][j]=INF;
    for (int i=1;i<=m;i++) {
        int u,v;
        double cost;
        scanf("%d%d%lf",&u,&v,&cost);
        map[u][v]=min(map[u][v],cost);
        e[i].u=u; e[i].v=v;
    }
    for (int i=1;i<=m;i++) {
        edge1 e1;
        e1.u=e[i].v; e1.len=map[e[i].u][e[i].v];
        G[e[i].u].push_back(e1);
    }
    
    double l=0.0,r=INF;
    for (int k=0;k<100;k++) {
        double mid=(l+r)/2;
        if (spfa(mid)) l=mid;
        else r=mid;
    }
    printf("%.2lf",l);
}