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
Source
Gold
Solution
这题是卡普通的SPFA,所以用一个小技巧:
if (dis[v]<dis[x[((t==2*n)?1:t+1)]]) swap(x[w],x[((t==2*n)?1:t+1)]);
竟然过了!好像可以快很多!开心!
(正解好像是dijkstra+拓扑)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 |
#include<iostream> #include<cstdio> #include<algorithm> #include<cstring> #include<cmath> #define ll long long using namespace std; int n,tot,Next[800100],head[200100],tree[800100],val[800100],x[400100],dis[200100]; struct node{int x,y,tot;}a[200100]; bool visit[200100]; bool cmpx(node a,node b) { return a.x<b.x;} bool cmpy(node a,node b) { return a.y<b.y;} void add(int x,int y,int z) { tot++; Next[tot]=head[x]; head[x]=tot; tree[tot]=y; val[tot]=z; } int calc(int x) { x++; if (x>2*n) x=1; return x; } int spfa() { for (int i=1;i<=n;i++) visit[i]=false,dis[i]=1<<30; int t=0,w=1; x[1]=1;visit[1]=true;dis[1]=0; while (t!=w) { t=calc(t); int u=x[t]; visit[u]=false; for (int i=head[u];i;i=Next[i]) { int v=tree[i]; if (dis[v]>dis[u]+val[i]) { dis[v]=dis[u]+val[i]; if (!visit[v]) { visit[v]=true,w=calc(w),x[w]=v; if (dis[v]<dis[x[calc(t)]]) swap(x[w],x[calc(t)]); } } } } return dis[n]; } int main() { scanf("%d",&n); for (int i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y),a[i].tot=i; sort(a+1,a+n+1,cmpx); tot=0; for (int i=1;i<n;i++) { add(a[i].tot,a[i+1].tot,a[i+1].x-a[i].x); add(a[i+1].tot,a[i].tot,a[i+1].x-a[i].x); } sort(a+1,a+n+1,cmpy); for (int i=1;i<n;i++) { add(a[i].tot,a[i+1].tot,a[i+1].y-a[i].y); add(a[i+1].tot,a[i].tot,a[i+1].y-a[i].y); } printf("%d\n",spfa()); return 0; } |