[bsoj3135] 慢跑问题

---

题目描述

　　这一天小K到附近的山上去慢跑，这座山上有N个凉亭，简单地编号为1到N（山顶的凉亭编号为N，山底的凉亭编号为1）。有趣的是，如果凉亭i和凉亭j满足i 这些凉亭中存在一些道路，连接两个海拔不同的凉亭，穿过一条道路需要花一定的时间。但这里还有一个十分诡异的事件，这里存在一些特殊道路。从这些道路经过 小K可以回到过去，也就是说小K需要花的时间为负值。通过这些道路，小K可以花较短的时间跑完他需要跑的路程。这样他就可以节约下许多时间，回到机房刷水题。
　　小K想从凉亭1跑到凉亭N去，但他不想走下坡路，也就是说他想要走一条海拔递增的路径，因为 小K始终认为递增的东西是具有美感的。
　　小K他想知道他最少需要花多少时间就可以从山底的凉亭到达山顶的凉亭？有一点需要注意，小K到达山顶的凉亭的时间可能会小于他出发的时间， 这就意味着他穿越了时空。这是一件多么美妙的事，于是 小K也想知道他是否能穿越时空。
他把这个任务交给了你，他希望你能告诉他最少需要花多少时间才能到达山顶，同时你也需要告诉他，他是否能穿越时空。

---

输入格式

　　输入数据第一行包含三个整数N，ML，MD（分别用一个空格隔开），分别表示凉亭的数目，普通道路的数目和特殊道路的数目；
　　第2行到第ML+1行，每行三个整数Ai，Bi，Ci（分别用一个空格隔开，1≤Ai,Bi≤N，Ai≠Bi），表示普通道路连接凉亭Ai和Bi，穿过这条道路需要花时间Ci；
　　第ML+2行到第ML+MD+1行，每行三个整数Ai，Bi，Ci（分别用一个空格隔开，1≤Ai,Bi≤N，Ai≠Bi），表示特殊道路连接凉亭Ai和Bi，穿过这条道路会回到相对时间Ci之前。

　　

输出格式

　　输出数据包含两行。
　　第一行包含一个整数T，表示小K最少花多少时间能到达山顶的凉亭；
　　第二行包含一个字符串“YES”或“NO”（不包含引号），如果小K可以回到过去那么输出“YES”，否则输出“NO”。

---

样例数据

样例输入

样例数据#1
4 3 2
1 2 1
1 3 4
3 4 3
2 4 1
1 4 2
样例数据#2
5 4 2
2 4 2
2 3 3
3 5 4
4 5 1
1 2 2
1 3 1

样例输出

样例数据#1
-2
YES
样例数据#2
1
NO

---

数据范围

　　对于10%的数据，N≤100，0≤ML,MD≤1000，Ci≤100；
　　对于30%的数据，N≤10000，0≤ML,MD≤30000，Ci≤500；
　　对于100%的数据，N≤100000，0≤ML,MD≤1000000，Ci≤2000。
　　输入数据保证从1到N会有一条路径。

---

题目分析

SPFA，插入正边负边，注意从海拔低的建单向边到海拔高的，否则边多了会TLE，最好加SLF优化和LLL优化。
当然也可以动规做。

---

源代码

#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
inline const int Get_Int() {
    int num=0,bj=1;
    char x=getchar();
    while(x<'0'||x>'9') {
        if(x=='-')bj=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9') {
        num=num*10+x-'0';
        x=getchar();
    }
    return num*bj;
}
const int maxn=200005;
struct Edge {
    int from,to;
    long long dist;
};
struct Spfa {
    int n,m;
    vector<Edge>edges;
    vector<int>G[maxn];
    bool inque[maxn];
    long long dist[maxn],used[maxn],path[maxn],pathe[maxn];
    void init(int n) {
        this->n=n;
        edges.clear();
        for(int i=1; i<=n; i++)G[i].clear();
    }
    void AddEdge(int from,int to,int dist) {
        edges.push_back((Edge) {
            from,to,dist
        });
        m=edges.size();
        G[from].push_back(m-1);
    }
    bool main(int s) {
        for(int i=1; i<=n; i++)dist[i]=0x7fffffff/2;
        memset(inque,0,sizeof(inque));
        deque<int>Q;
        Q.push_back(s);
        dist[s]=0;
        path[s]=s;
        inque[s]=1;
        used[s]++;
        while(!Q.empty()) {
            int Now=Q.front();
            Q.pop_front();
            inque[Now]=0;
            for(int i=0; i<G[Now].size(); i++) {
                Edge& e=edges[G[Now][i]];
                int Next=e.to;
                if(Next<=Now)continue;
                if(dist[Next]>dist[Now]+e.dist) {
                    dist[Next]=dist[Now]+e.dist;
                    path[Next]=Now;
                    pathe[Next]=G[Now][i];
                    if(!inque[Next]) {
                        used[Next]++;
                        if(used[Next]==edges.size())return false; //负权回环
                        if(!Q.empty()&&dist[Next]<dist[Q.front()])Q.push_front(Next); //SLF优化
                        else Q.push_back(Next);
                        inque[Next]=1;
                    }
                }
            }
        }
        return true;
    }
};
Spfa spfa;
int n,m1,m2;
int main() {
    n=Get_Int();
    m1=Get_Int();
    m2=Get_Int();
    spfa.init(n);
    for(int i=1; i<=m1; i++) {
        int x=Get_Int(),y=Get_Int(),v=Get_Int();
        if(y>x)spfa.AddEdge(x,y,v);
        else spfa.AddEdge(y,x,v);
    }
    for(int i=1; i<=m2; i++) {
        int x=Get_Int(),y=Get_Int(),v=Get_Int();
        if(y>x)spfa.AddEdge(x,y,-v);
        else spfa.AddEdge(y,x,-v);
    }
    spfa.main(1);
    printf("%lld\n",spfa.dist[n]);
    if(spfa.dist[n]<0)puts("YES");
    else puts("NO");
    return 0;
}

---