Spfa+DP【p2149】[SDOI2009]Elaxia的路线

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 ≤ y2 ≤ N),分别表示Elaxia的宿舍和实验室及w**的宿舍和实验室的标号(两对点分别 x1,y1和x2,y2)。

接下来M行:每行三个整数,u、v、l(1 ≤ u ≤ N,1 ≤ v ≤ N,1 ≤ l ≤ 10000),表 u和v之间有一条路,经过这条路所需要的时间为l。

Output

一行,一个整数,表示每天两人在一起的时间(即最长公共路径的长度)

一句话题意,求两对点间最短路的最长公共路径.

首先需要明确的是,我们需要跑最短路.

先跑出第一组点对之间的最短路,标记最短路上的边,然后再跑第二组点对的最短路,最后需要再小小的记录一下答案即可.

注意\(to[u]\)代表到达\(u\)这个节点的最长公共路径.

每次取\(max\)即可,应该不是很难理解,就不过多解释了.

代码

#include<cstdio>
#include<cctype>
#include<cstring>
#include<queue>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#define int long long
#define clear(a,b) memset(a,b,sizeof a)
#define N 2005
#define R register
using namespace std;
inline void in(int &x)
{
    int f=1;x=0;char s=getchar();
    while(!isdigit(s)){if(s=='-')f=-1;s=getchar();}
    while(isdigit(s)){x=x*10+s-'0';s=getchar();}
    x*=f;
}
int n,m,a,b,c,d,dis[N],head[N],tot=-1;
struct cod{int u,v,w,flg;}edge[N*N];
bool vis[N];
inline void add(int x,int y,int z)
{
    edge[++tot].u=head[x];
    edge[tot].v=y;
    edge[tot].w=z;
    head[x]=tot;
}
inline void spfa(int s)
{
    clear(dis,0x7f);clear(vis,0);
    dis[s]=0;vis[s]=true;
    queue<int>q;q.push(s);
    while(!q.empty())
    {
        int u=q.front();q.pop();vis[u]=false;
        for(R int i=head[u];i!=-1;i=edge[i].u)
        {
            if(dis[edge[i].v]>dis[u]+edge[i].w)
            {
                dis[edge[i].v]=dis[u]+edge[i].w;
                if(!vis[edge[i].v])
                {
                    vis[edge[i].v]=true;
                    q.push(edge[i].v);
                }
            }
        }
    }
}
inline void calc()
{
    clear(vis,0);
    queue<int>q;q.push(b);vis[b]=true;
    while(!q.empty())
    {
        int u=q.front();q.pop();
        for(R int i=head[u];i!=-1;i=edge[i].u)
        {
            if(dis[edge[i].v]+edge[i].w==dis[u])
            {
                edge[i].flg=edge[i^1].flg=1;
                if(!vis[edge[i].v])
                {
                    vis[edge[i].v]=true;
                    q.push(edge[i].v);
                }
            }
        }
        
    }
}
bool inq[N];
int to[N];
inline int ans()
{
    R int ans=0;clear(vis,0);
    queue<int>q;q.push(d);vis[d]=inq[d]=true;
    while(!q.empty())
    {
        int u=q.front();q.pop();inq[u]=false;
        for(R int i=head[u];i!=-1;i=edge[i].u)
        {
            if(dis[edge[i].v]+edge[i].w==dis[u])
            {
                if(edge[i].flg and to[edge[i].v]<to[u]+edge[i].w)
                {
                    to[edge[i].v]=to[u]+edge[i].w;
                    ans=max(ans,to[edge[i].v]);
                    if(!inq[edge[i].v])
                    {
                        inq[edge[i].v]=true;
                        q.push(edge[i].v);
                    }
                }
                if(!vis[edge[i].v])
                {
                    vis[edge[i].v]=inq[edge[i].v]=true;
                    q.push(edge[i].v);
                }
            }
        }
    }
    return ans;
}
signed main()
{
    in(n),in(m);clear(head,-1);
    in(a),in(b),in(c),in(d);
    for(R int i=1,x,y,z;i<=m;i++)
    {
        in(x),in(y),in(z);
        add(x,y,z);add(y,x,z);
    }
    spfa(a);;calc();spfa(c);
    printf("%lld\n",ans());
}

猜你喜欢

转载自www.cnblogs.com/-guz/p/9822639.html
今日推荐