【BZOJ 1880】 Elaxia的路线

【分析】

我们分别从两个起点，两个终点，做四次最短路。

然后我们可以知道哪些点和边是两人都经过的，并且按照经过的次序使边变为有向的。

这样就构成了一个DAG。毫无疑问，我们要求的是这个DAG的最长路。

用拓扑排序+DP即可完成。

(PS：不知道标解是怎样，感觉我的做法很奇葩，但是在BZOJ上排rank4)

【代码】

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>
#include <map>
#include <utility>
#include <functional>
#include <vector>

using namespace std;

#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define rrep(i,b,a) for(int i=(b);i>=(a);--i)
#define INF (~0U>>3)
#define pii pair<int,int> 
#define m_p make_pair
#define p_b push_back
#define sf scanf
#define pf printf
#define MAXN 1510
#define MAXM 300010

struct edge{
	int v,w;
	edge* next;
}*head[2][MAXN],Edge[2][MAXM<<1];
int Idx[2];
int N,M;
int S1,T1;
int S2,T2;
int Dis1[MAXN];
int Dit1[MAXN];
int Dis2[MAXN];
int Dit2[MAXN];
int degree[MAXN];
int f[MAXN];
bool vis[MAXN];

void Read(int& x,bool mark=0){
	char tt=getchar();
	for(;tt<'0'||'9'<tt;tt=getchar()) if(tt=='-') mark=1;
	for(x=0;'0'<=tt&&tt<='9';x=(x<<1)+(x<<3)+tt-'0',tt=getchar());
	x=mark?-x:x;
}

void Addedge(int u,int v,int w,int k){
	int idx=++Idx[k];
	Edge[k][idx].v=v;
	Edge[k][idx].w=w;
	Edge[k][idx].next=head[k][u];
	head[k][u]=Edge[k]+idx;
}

void Init(){
	Read(N);
	Read(M);
	Read(S1);Read(T1);
	Read(S2);Read(T2);
	int u,v,w;
	while(M--){
		Read(u);
		Read(v);
		Read(w);
		Addedge(u,v,w,0);
		Addedge(v,u,w,0);
	}
}

void Dijkstra(int* dis,int S){
	priority_queue<pii,vector<pii>,greater<pii> >Pq;
	rep(i,1,N) dis[i]=INF;
	dis[S]=0;
	Pq.push(m_p(dis[S],S));
	while(!Pq.empty()){
		int u=Pq.top().second;
		Pq.pop();
		if(vis[u]) continue;
		vis[u]=true;
		for(edge* p=head[0][u];p;p=p->next){
			int v=p->v;
			if((!vis[v])&&dis[v]>dis[u]+p->w){
				dis[v]=dis[u]+p->w;
				Pq.push(m_p(dis[v],v));
			}
		}
	}
	memset(vis,0,sizeof vis);
}

int Dp(){
	int rt=0;
	queue <int> Q;
	rep(i,1,N)
		if(vis[i]&&(!degree[i]))
			Q.push(i);
	while(!Q.empty()){
		int u=Q.front();
		Q.pop();
		for(edge* p=head[1][u];p;p=p->next){
			int v=p->v;
			f[v]=max(f[v],f[u]+p->w);
			rt=max(rt,f[v]);
			if(!(--degree[v])) Q.push(v);
		}
	}
	return rt;
}

void Solve(){
	if(S1==T1||S2==T2){
		puts("0");
		return;
	}
	Dijkstra(Dis1,S1);
	Dijkstra(Dit1,T1);
	Dijkstra(Dis2,S2);
	Dijkstra(Dit2,T2);
	
	rep(i,1,N)
		vis[i]=(Dis1[i]+Dit1[i]==Dis1[T1]&&Dis2[i]+Dit2[i]==Dis2[T2]);
	rep(i,1,N)
		if(vis[i])
			for(edge* p=head[0][i];p;p=p->next){
				int v=p->v;
				if(vis[v]&&Dis1[v]==Dis1[i]+p->w){
					Addedge(i,v,p->w,1);
					degree[v]++;
				}
			}
	
	pf("%d\n",Dp());
}

int main(){
	Init();
	Solve();
	return 0;
}