BZOJ 2200: [Usaco2011 Jan]道路和航线

Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 900 Solved: 338
[Submit][Status][Discuss]
Description
Farmer John正在一个新的销售区域对他的牛奶销售方案进行调查。他想把牛奶送到T个城镇 (1 <= T <= 25,000)，编号为1T。这些城镇之间通过R条道路 (1 <= R <= 50,000，编号为1到R) 和P条航线 (1 <= P <= 50,000，编号为1到P) 连接。每条道路i或者航线i连接城镇A_i (1 <= A_i <= T)到B_i (1 <= B_i <= T)，花费为C_i。对于道路，0 <= C_i <= 10,000;然而航线的花费很神奇，花费C_i可能是负数(-10,000 <= C_i <= 10,000)。道路是双向的，可以从A_i到B_i，也可以从B_i到A_i，花费都是C_i。然而航线与之不同，只可以从A_i到B_i。事实上，由于最近恐怖主义太嚣张，为了社会和谐，出台 了一些政策保证：如果有一条航线可以从A_i到B_i，那么保证不可能通过一些道路和航线从B_i回到A_i。由于FJ的奶牛世界公认十分给力，他需要运送奶牛到每一个城镇。他想找到从发送中心城镇S(1 <= S <= T) 把奶牛送到每个城镇的最便宜的方案，或者知道这是不可能的。
Input
* 第1行：四个空格隔开的整数: T, R, P, and S * 第2到R+1行：三个空格隔开的整数（表示一条道路）：A_i, B_i 和 C_i * 第R+2到R+P+1行：三个空格隔开的整数（表示一条航线）：A_i, B_i 和 C_i
Output
* 第1到T行：从S到达城镇i的最小花费，如果不存在输出”NO PATH”。
Sample Input

6 3 3 4

1 2 5

3 4 5

5 6 10

3 5 -100

4 6 -100

1 3 -10

样例输入解释：

一共六个城镇。在1-2，3-4，5-6之间有道路，花费分别是5，5，10。同时有三条航线：3->5，

4->6和1->3，花费分别是-100，-100，-10。FJ的中心城镇在城镇4。

Sample Output

NO PATH

NO PATH

5

0

-95

-100

样例输出解释：

FJ的奶牛从4号城镇开始，可以通过道路到达3号城镇。然后他们会通过航线达到5和6号城镇。

但是不可能到达1和2号城镇。

解题思路

可以用SLF优化的spfa强行水过，SLF就是用一个双端队列，如果当前被松弛后的dis小于队头，就把当前的点在队头入队，否则在队尾入队。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<deque>

using namespace std;
const int MAXN = 25005;
const int MAXM = 100005;
const int inf = 0x3f3f3f3f;

inline int rd(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int T,R,P,S;
int head[MAXN],cnt;
int to[MAXM<<1],nxt[MAXM<<1],val[MAXM<<1];
int dis[MAXN];
deque<int> Q;
bool vis[MAXN];

inline void add(int bg,int ed,int v){
    to[++cnt]=ed,nxt[cnt]=head[bg],val[cnt]=v,head[bg]=cnt;
}

inline void spfa(){
    for(register int i=1;i<=T;i++) dis[i]=inf;
    Q.push_back(S);vis[S]=1;dis[S]=0;
    while(!Q.empty()){
        int x=Q.front();Q.pop_front();
        vis[x]=0;
        for(register int i=head[x];i;i=nxt[i]){
            int u=to[i];
            if(dis[u]>dis[x]+val[i]){
                dis[u]=dis[x]+val[i];
                if(!vis[u]){
                    if(dis[u]<dis[Q.front()]) Q.push_front(u);
                    else Q.push_back(u);
                    vis[u]=1;
                }   
            }
        }
    }
}

int main(){
    T=rd(),R=rd(),P=rd(),S=rd();
    for(register int i=1;i<=R;i++){
        int x=rd(),y=rd(),z=rd();
        add(x,y,z);add(y,x,z);
    }
    for(register int i=1;i<=P;i++){
        int x=rd(),y=rd(),z=rd();
        add(x,y,z);
    }spfa();
    for(register int i=1;i<=T;i++){
        if(dis[i]==inf) puts("NO PATH");
        else printf("%d\n",dis[i]);
    }
    return 0;
}