bzoj2200 [Usaco2011 Jan]道路和航线 最短路

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2200

2200: [Usaco2011 Jan]道路和航线

Time Limit: 10 Sec   Memory Limit: 259 MB
Submit: 724   Solved: 255
[ 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号城镇。

HINT

这题就是裸最短路，只不过有双向边和单向边之分罢了。

听说裸spfa会T，所以加上SLL优化就AC了。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<map>
#include<set>
#include<stack>
using namespace std;
typedef long long LL;
const int inf=1e9;  
inline int read()  
{  
    int x=0,f=1; char ch=getchar();  
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}  
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}  
    return x*f;  
}
int n,m1,m2,s;
int u,v,w;
bool vis[25005];
int dis[25005];
struct node{
	int to,val;
};
vector<node> edge[25005];
void ins(int x,int y,int c)
{
	node xxx;
	xxx.to=y;xxx.val=c;
	edge[x].push_back(xxx);
}
void spfa(int s)
{
	for(int i=1;i<=n;i++)
	dis[i]=inf;
	deque<int> q;
	q.push_back(s);
	dis[s]=0;
	vis[s]=1;
	while(!q.empty())
	{
		int now=q.front();
		q.pop_front();
		vis[now]=0;
		int lens=edge[now].size();
		for(int i=0;i<lens;i++)
		{
			node nex=edge[now][i];
			if(dis[nex.to]>dis[now]+nex.val)
			{
				dis[nex.to]=dis[now]+nex.val;
				if(vis[nex.to])
				continue;  
	            vis[nex.to]=1;
	            if(!q.empty())  
	            {  
	                if(dis[nex.to]<dis[q.front()])
					q.push_front(nex.to);  
	                else q.push_back(nex.to);  
	            }  
	            else q.push_back(nex.to);
			}
		}
		
	}
}
int main()  
{   
    n=read(),m1=read(),m2=read(),s=read();  
    for(int i=1;i<=m1;i++)  
    {  
        u=read(),v=read(),w=read();  
        ins(u,v,w),ins(v,u,w);  
    }
    for(int i=1;i<=m2;i++)  
    {
    	u=read(),v=read(),w=read();
		ins(u,v,w); 
	}
	spfa(s);
    for(int i=1;i<=n;i++)  
    if(dis[i]==inf)
	printf("NO PATH\n");  
    else
	printf("%d\n",dis[i]);  
    return 0;  
}