洛谷P3008「Usaco2011 Jan」道路和航线

乍一看，这题非常简单
不就是一个单源最短路径吗
但是定睛一看，怎么有负边，好像没法写Dijkstra
spfa呢？
靠，usaco日常卡spfa
但实际上，spfa+SLF可以过…因为数据比较水…

但是，既然出题人都刻意的卡我们了，我们自然要按他的意思来写

我们仔细思考这个图的构成，对于道路，全部是正权双向边
航线有负权，但是一定不会构成环

缩完点好像可以跑拓扑排序+DP啊

对于一个环内，因为题目保证，不可能出现负权环，那么我们可以跑Dijkstra
对于不是一个环内的，我们按照拓扑序合并就好了

说是这么说，实际上实现起来细节很多
看代码吧

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<queue>
#include<vector>
#include<cstring>
using namespace std;
const int maxn = 2000007;
const int INF = 0x3f3f3f3f;
int n,m1,m2,s;
struct node{
    int to,next,w;
}edge[maxn];
int cnt,head[maxn];
void add(int from,int to,int w){
    edge[++cnt].to=to;
    edge[cnt].next=head[from];
    head[from]=cnt;
    edge[cnt].w=w;
}
int C[maxn];
vector<int>ma[maxn];
int vis[maxn],tot;
void get_c(int u){
    C[u]=tot;ma[tot].push_back(u);
    vis[u]=1;
    for(int i=head[u];i;i=edge[i].next){
        int to=edge[i].to;
        if(vis[to])continue;
        get_c(to);
    }
}
int rd[maxn];
int dis[maxn];
 
 
int main(){
    scanf("%d%d%d%d",&n,&m1,&m2,&s);
    for(int i=1;i<=m1;i++){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z),add(y,x,z);
    }
    for(int i=1;i<=n;i++){
        if(!vis[i]){
            tot++;
            get_c(i);
        }
    }
    for(int i=1;i<=m2;i++){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z);
        rd[C[y]]++;
    }
     
     
    memset(vis,0,sizeof(vis));
    memset(dis,0x7f,sizeof(dis));
    dis[s]=0;
    queue<int>Q;
    Q.push(C[s]);
    for(int i=1;i<=tot;i++){if(!rd[i])Q.push(i);}
    priority_queue< pair<int,int> >q;
    while(!Q.empty()){
        int now=Q.front();Q.pop();
        for(int i=0;i<ma[now].size();i++){
            q.push(make_pair(-dis[ma[now][i]],ma[now][i]));
        }
        while(!q.empty()){
            int f1=q.top().first;
            int f2=q.top().second;
            q.pop();
            //if(f1!=-dis[f2])continue;
            if(vis[f2])continue;
            vis[f2]=1;
            for(int j=head[f2];j;j=edge[j].next){
                int to=edge[j].to;
                int w=edge[j].w;
                if(dis[to]>dis[f2]+w){
                    dis[to]=dis[f2]+w;
                    if(C[to]==C[f2])q.push(make_pair(-dis[to],to));
                }
                if(C[to]!=C[f2]&&(--rd[C[to]])==0)Q.push(C[to]);
            }
        }
    }
    for(int i=1;i<=n;i++)
    if(dis[i]>INF) cout<<"NO PATH"<<endl;
    else cout<<dis[i]<<endl;
}