道路和航路

【题目描述】
农夫约翰正在针对一个新区域的牛奶配送合同进行研究。他打算分发牛奶到T个城镇（标号为1..T），这些城镇通过R条标号为（1..R）的道路和P条标号为（1..P）的航路相连。
每一条公路i或者航路i表示成连接城镇Ai（1<=A_i<=T）和Bi（1<=Bi<=T）代价为Ci。每一条公路，Ci的范围为0<=Ci<=10,000；由于奇怪的运营策略，每一条航路的Ci可能为负的，也就是-10,000<=Ci<=10,000。
每一条公路都是双向的，正向和反向的花费是一样的，都是非负的。
每一条航路都根据输入的Ai和Bi进行从Ai->Bi的单向通行。实际上，如果现在有一条航路是从Ai到Bi的话，那么意味着肯定没有通行方案从Bi回到Ai。
农夫约翰想把他那优良的牛奶从配送中心送到各个城镇，当然希望代价越小越好，你可以帮助他嘛？配送中心位于城镇S中（1<=S<=T）。
【输入格式】
输入的第一行包含四个用空格隔开的整数T，R，P，S。
接下来R行，描述公路信息，每行包含三个整数，分别表示Ai，Bi和Ci。
接下来P行，描述航路信息，每行包含三个整数，分别表示Ai，Bi和Ci。
【输出格式】
输出T行，分别表示从城镇S到每个城市的最小花费，如果到不了的话输出NO PATH。
【样例输入】

6 3 3 4
1 2 5
3 4 5
5 6 10
3 5 -100
4 6 -100
1 3 -10

【样例输出】

NO PATH
NO PATH
5
0
-95
-100

【数据范围】
对于20%的数据，T<=100，R<=500，P<=500；
对于30%的数据，R<=1000，R<=10000，P<=3000；
对于100%的数据，1<=T<=25000，1<=R<=50000，1<=P<=50000。
【分析】
先上SPFA的代码：

#include<bits/stdc++.h>
using namespace std;
#define MAXN 200000
#define INF  1<<30
struct Edge{
  int to,next,dis;
}a[MAXN];
int m=1,d[25005],head[25005],vis[25005],T,R,P,S;

inline int read(){
    int x=0,p=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-')p=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*p;
}

inline void write(int x){
    char ch[15]={};
    int i=0;
    if (x<0) putchar('-'),x=-x;
    do {ch[++i]='0'+x%10;} while (x/=10);
    while (i) putchar(ch[i--]);
    putchar('\n');
}

int main(){
    T=read(),R=read(),P=read(),S=read();
    for(int i=0;i<=T;i++) head[i]=-1,vis[i]=0,d[i]=INF;
    int u,v,dis;
    for(int i=1;i<=R;i++){
        u=read(),v=read(),dis=read();
      a[m].to=v,a[m].dis=dis,a[m].next=head[u],head[u]=m++;
        a[m].to=u,a[m].dis=dis,a[m].next=head[v],head[v]=m++;
    }
    for(int i=1;i<=P;i++){
        u=read(),v=read(),dis=read();
      a[m].to=v,a[m].dis=dis,a[m].next=head[u],head[u]=m++;
    }

    vis[S]=1;
    d[S]=0;
    deque<int> q;
    q.push_front(S);
    int sum=0,cnt=!;
    while(!q.empty()){
        int u=q.front();
        q.pop_front();
        if (dis[u]*cnt<=sum){
          q.push_back(u);
          continue;
        }
        vis[u]=0;
        sum-=dis[u],cnt--;
        for(int i=head[u];i!=-1;i=a[i].next){
            int v=a[i].to;
            if(d[v]>d[u]+a[i].dis){
                d[v]=d[u]+a[i].dis;
                if (q.empty()) q.push_front(v),vis[v]=1;
                else if(!vis[v]){
                    if (d[v]<=d[q.front()]) q.push_front(v); else q.push_back(v);
                    vis[v]=1;
                    sum+=dis[v],cnt++;
                }
            }
        }
    }

    for(int i=1;i<=T;i++)
        if (d[i]==INF) printf("NO PATH\n"); else write(d[i]);
}

然后我们想到了Dijkstra。
众：不是说Dijkstra不能处理含有负权的图吗？
我：这题的图比较特殊啊~
不知道各位有没有思考为什么要读入R条公路和P条航路而不是直接读入2R+P条有向的路呢？其实这就是和SPFA模板题唯一的区别啊突破口。
暂时抛开单向边不谈，我们将无向边连成的各点分成很多块状结构，然后在块与块之间用单向边相连。
然后可以用缩点法将块缩成点，于是现在的图变成了一个DAG。
然后就可以无脑DP咯~~
至于具体做法，可以用并查集维护一下连成的块，然后用Dijkstra求出每个块内部的最短路（从起点开始）。对于第i个块，当所有存在有向边指向该块的块内部的最短路已经求好的时候，再求第i个块的最短路。
可以使用deque（STL里的双向队列）优化程序。