1. 程式人生 > >BZOJ 1880: [Sdoi2009]Elaxia的路線

BZOJ 1880: [Sdoi2009]Elaxia的路線

sdoi esp tar 學校 problems logs tput 兩個人 cst

Time Limit: 4 Sec Memory Limit: 64 MB
Submit: 1461 Solved: 574
[Submit][Status][Discuss]

Description

最近,Elaxia和w**的關系特別好,他們很想整天在一起,但是大學的學習太緊張了,他們 必須合理地安排兩個人在一起的時間。Elaxia和w**每天都要奔波於宿舍和實驗室之間,他們 希望在節約時間的前提下,一起走的時間盡可能的長。 現在已知的是Elaxia和w**所在的宿舍和實驗室的編號以及學校的地圖:地圖上有N個路 口,M條路,經過每條路都需要一定的時間。 具體地說,就是要求無向圖中,兩對點間最短路的最長公共路徑。

Input

第一行:兩個整數N和M(含義如題目描述)。 第二行:四個整數x1、y1、x2、y2(1 ≤ x1 ≤ N,1 ≤ y1 ≤ N,1 ≤ x2 ≤ N,1 ≤ ≤ N),分別表示Elaxia的宿舍和實驗室及w**的宿舍和實驗室的標號(兩對點分別 x1,y1和x2,y2)。 接下來M行:每行三個整數,u、v、l(1 ≤ u ≤ N,1 ≤ v ≤ N,1 ≤ l ≤ 10000),表 u和v之間有一條路,經過這條路所需要的時間為l。 出出出格格格式式式::: 一行,一個整數,表示每天兩人在一起的時間(即最長公共路徑的長度)。

Output

一行,一個整數,表示每天兩人在一起的時間(即最長公共路徑的長度)

Sample Input

9 10
1 6 7 8
1 2 1
2 5 2
2 3 3
3 4 2
3 9 5
4 5 3
4 6 4
4 7 2
5 8 1
7 9 1

Sample Output

3

HINT

對於30%的數據,N ≤ 100;
對於60%的數據,N ≤ 1000;
對於100%的數據,N ≤ 1500,輸入數據保證沒有重邊和自環。

Source

Day2

拓撲排序+最短路 屠龍寶刀點擊就送
#include <cstring>
#include 
<cstdio> #include <cctype> #include <queue> #define M 1000005 #define N 1505 using namespace std; struct Edge { int x,y,z; }e[M]; int dist[2][N],D[N],nextt[M<<1],to[M<<1],head[M<<1],val[M<<1],cnt,n,m,s[3],t[3],f[N],use[M]; inline void Read(int &x) { register char ch=getchar(); for(x=0;!isdigit(ch);ch=getchar()); for(;isdigit(ch);x=x*10+ch-0,ch=getchar()); } inline void ins(int u,int v,int l) { nextt[++cnt]=head[u]; to[cnt]=v; val[cnt]=l; head[u]=cnt; } bool vis[N]; struct node { int x,y; bool operator<(node a)const { return y>a.y; } }; void dijkstra(int s,int *dis) { memset(dis,0x3f,sizeof(dist[0])); memset(vis,0,sizeof(vis)); dis[s]=0; priority_queue<node>q; q.push((node){s,dis[s]}); for(node now;!q.empty();) { now=q.top(); q.pop(); if(vis[now.x]) continue; vis[now.x]=1; for(int i=head[now.x];i;i=nextt[i]) { int v=to[i]; if(dis[v]>dis[now.x]+val[i]) { dis[v]=dis[now.x]+val[i]; q.push((node){v,dis[v]}); } } } } inline void swap(int &x,int &y) {int tmp=y;y=x;x=tmp;} inline int min(int x,int y) {return x>y?y:x;} inline int max(int x,int y) {return x>y?x:y;} void calc(int s,int t) { for(int i=1;i<=m;++i) { int x=e[i].x,y=e[i].y,z=e[i].z; if(min(dist[0][x],dist[0][y])+min(dist[1][x],dist[1][y])+z==dist[0][t]) use[i]++; } } int main(int argc,char *argv[]) { Read(n);Read(m); for(int i=1;i<=2;++i) Read(s[i]),Read(t[i]); for(int a,b,c,i=1;i<=m;++i) { Read(a);Read(b);Read(c); e[i]=(Edge){a,b,c}; ins(a,b,c); ins(b,a,c); } for(int i=1;i<=2;++i) { dijkstra(s[i],dist[0]); dijkstra(t[i],dist[1]); calc(s[i],t[i]); } memset(head,0,sizeof(head)); cnt=0; for(int i=1;i<=m;++i) { if(use[i]==2) { int x=e[i].x,y=e[i].y,z=e[i].z; if(dist[0][x]>dist[0][y]) swap(x,y); ins(x,y,z);D[y]++; } } int h=0,t=0,Que[N]; for(int i=1;i<=n;++i) if(!D[i]) Que[++t]=i; for(int now;h<t;) { now=Que[++h]; for(int i=head[now];i;i=nextt[i]) { int v=to[i]; f[v]=max(f[now]+val[i],f[v]); D[v]--; if(D[v]==0) Que[++t]=v; } } int ans=0; for(int i=1;i<=n;++i) ans=max(ans,f[i]); printf("%d\n",ans); return 0; }

BZOJ 1880: [Sdoi2009]Elaxia的路線