【BZOJ1880】[Sdoi2009]Elaxia的路线(最短路)

【BZOJ1880】[Sdoi2009]Elaxia的路线(最短路)

题面

BZOJ
洛谷

题解

假装我们知道了任意两点间的最短路,那么我们怎么求解答案呢?
不难发现公共路径一定是一段连续的路径(如果不连续那么显然可以把中间分开的那段变成一样路径)。
这样子我们只需要\(O(n^2)\)枚举这个路径的起点和终点\(check\)一下就可以知道答案了。
然而并没有办法求任意两点之间的最短路,因为这样子是\(O(n^3)\)的。并不要认为对于每个点跑一遍最短路就行了,\(SPFA\)不说了,它死了。\(Dijkstra\)的复杂度是\(O((n+m)logm)\)的,如果\(m=n^2\)还不如\(Floyd\)
大概的代码如下,我是对于每一个点跑了一遍\(Dij\)(雾

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
#define MAX 1550
#define pi pair<int,int>
#define mp make_pair
#define fr first
#define sd second
inline int read()
{
    int x=0;bool t=false;char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=true,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return t?-x:x;
}
struct Line{int v,next,w;}e[MAX*MAX*2];
int h[MAX],cnt=1;
inline void Add(int u,int v,int w){e[cnt]=(Line){v,h[u],w};h[u]=cnt++;}
int dis[MAX][MAX];
int n,m,S1,S2,T1,T2,ans;
bool vis[MAX];
void SPFA(int S,int *dis)
{
    for(int i=1;i<=n;++i)dis[i]=1e9,vis[i]=false;
    priority_queue<pi,vector<pi>,greater<pi> > Q;
    Q.push(mp(0,S));dis[S]=0;
    while(!Q.empty())
    {
        pi u=Q.top();Q.pop();
        if(vis[u.sd])continue;vis[u.sd]=true;
        for(int i=h[u.sd];i;i=e[i].next)
            if(dis[e[i].v]>u.fr+e[i].w)
            {
                dis[e[i].v]=u.fr+e[i].w;
                Q.push(mp(dis[e[i].v],e[i].v));
            }
    }
}
int Dis1(int i,int j){return dis[S1][i]+dis[i][j]+dis[j][T1];}
int Dis2(int i,int j){return dis[S2][i]+dis[i][j]+dis[j][T2];}
int main()
{
    n=read();m=read();
    S1=read(),T1=read(),S2=read(),T2=read();
    for(int i=1;i<=m;++i)
    {
        int u=read(),v=read(),w=read();
        Add(u,v,w);Add(v,u,w);
    }
    for(int i=1;i<=n;++i)SPFA(i,dis[i]);
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j)
            if(i!=j)
            {
                if(Dis1(i,j)>dis[S1][T1]&&Dis1(j,i)>dis[S1][T1])continue;
                if(Dis2(i,j)>dis[S2][T2]&&Dis2(j,i)>dis[S2][T2])continue;
                ans=max(ans,dis[i][j]);
            }
    printf("%d\n",ans);
    return 0;
}

首先,对于钦定的起点,我们到达终点一定是沿着最短路径\(DAG\)移动,那么我们考虑对于起点构建最短路径\(DAG\),将几个\(DAG\)重叠之后显然是要求的是\(S1\rightarrow T1\)的一条路径上的最长的、在两个最短路径\(DAG\)中都出现过的链。
那么我们把同时出现在两个\(DAG\)中的边全部拿出来建图,跑一个\(DAG\)上最长链就好了。
我代码写的有点冗长。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
#define MAX 1550
#define pi pair<int,int>
#define mp make_pair
#define fr first
#define sd second
inline int read()
{
    int x=0;bool t=false;char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=true,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return t?-x:x;
}
struct Line{int v,next,w;}e[1000000];
int h[MAX],cnt=2,dg[MAX];
inline void Add(int u,int v,int w){e[cnt]=(Line){v,h[u],w};h[u]=cnt++;}
vector<int> E[MAX],W[MAX];
int dis[MAX][MAX];
int n,m,S1,S2,T1,T2,ans;
bool vis[MAX];
void SPFA_heap(int S,int *dis)
{
    for(int i=1;i<=n;++i)dis[i]=1e9,vis[i]=false;
    priority_queue<pi,vector<pi>,greater<pi> > Q;
    Q.push(mp(0,S));dis[S]=0;
    while(!Q.empty())
    {
        pi u=Q.top();Q.pop();
        if(vis[u.sd])continue;vis[u.sd]=true;
        for(int i=h[u.sd];i;i=e[i].next)
            if(dis[e[i].v]>u.fr+e[i].w)
            {
                dis[e[i].v]=u.fr+e[i].w;
                Q.push(mp(dis[e[i].v],e[i].v));
            }
    }
}
int S[MAX],f[MAX],top;
void Topsort()
{
    queue<int> Q;
    for(int i=1;i<=n;++i)if(!dg[i])Q.push(i);
    while(!Q.empty())
    {
        int u=Q.front();Q.pop();S[++top]=u;
        for(int i=0,l=E[u].size();i<l;++i)
            if(!--dg[E[u][i]])Q.push(E[u][i]);
    }
    for(int i=top;i;--i)
        for(int j=0,l=E[S[i]].size();j<l;++j)
            f[S[i]]=max(f[S[i]],f[E[S[i]][j]]+W[S[i]][j]);
    for(int i=1;i<=top;++i)ans=max(ans,f[S[i]]);
}
int main()
{
    n=read();m=read();
    S1=read(),T1=read(),S2=read(),T2=read();
    for(int i=1;i<=m;++i)
    {
        int u=read(),v=read(),w=read();
        Add(u,v,w);Add(v,u,w);
    }
    SPFA_heap(S1,dis[S1]);SPFA_heap(T1,dis[T1]);SPFA_heap(S2,dis[S2]);SPFA_heap(T2,dis[T2]);
    for(int u=1;u<=n;++u)
        for(int i=h[u];i;i=e[i].next)
            if(i&1)
            {
                int u=e[i].v,v=e[i^1].v;
                if(min(dis[S1][u]+dis[T1][v]+e[i].w,dis[S1][v]+dis[T1][u]+e[i].w)>dis[S1][T1])continue;
                if(min(dis[S2][u]+dis[T2][v]+e[i].w,dis[S2][v]+dis[T2][u]+e[i].w)>dis[S2][T2])continue;
                if(dis[S1][u]+dis[T1][v]+e[i].w==dis[S1][T1])E[u].push_back(v),W[u].push_back(e[i].w),dg[v]+=1;
                else E[v].push_back(u),W[v].push_back(e[i].w),dg[u]+=1;
            }
    Topsort();
    printf("%d\n",ans);
    return 0;
}

转载于:https://www.cnblogs.com/cjyyb/p/9800218.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值