【BZOJ 1880】 Elaxia的路线

来源：互联网发布：帆布背包知乎编辑：程序博客网时间：2024/05/16 08:30

【分析】

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

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

这样就构成了一个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 300010struct 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;}

0 0