【CSP模拟赛】益智游戏（最短路（DJSPFA）&拓扑排序）

题目描述

小P和小R在玩一款益智游戏。游戏在一个正权有向图上进行。小P 控制的角色要从A 点走最短路到B 点，小R 控制的角色要从C 点走最短路到D 点。一个玩家每回合可以有两种选择，移动到一个相邻节点或者休息一回合。假如在某一时刻，小P 和小R 在相同的节点上，那么可以得到一次特殊奖励，但是在每个节点上最多只能得到一次。求最多能获得多少次特殊奖励

输入格式

第一行两个整数n，m 表示有向图的点数和边数。接下来m 行每行三个整数xi，yi，li，表示从xi 到yi 有一条长度为li 的边。最后一行四个整数A，B，C，D，描述小P 的起终点，小R 的起终点。

输出格式

输出一个整数表示最多能获得多少次特殊奖励。若小P不能到达B点或者小R不能到达D点则输出-1。

输入样例

5 5
1 2 1
2 3 2
3 4 4
5 2 3
5 3 5
1 3 5 4

输出样例

2
提示
对于30%的数据，满足n≤50 对于60%的数据，满足n≤1000，m≤5000 对于100%的数据，满足n≤50000，m≤200000，1≤li≤500000000

分析

A到B，C到D的最短路有很多条。

但对于任意一条A到B的最短路和任意一条C到D的最短路如果他们有交点，那交点一定是连续的，否则一定存在交点更多的两条最短路。

如图

如果A到B的路径（红色）与C到D的路径（蓝色）有不连续的交点E与F

那么E通过路径1到达F与E通过路径2到达F的距离应该是一样，否则路径A-B，路径C-D中有一条不是最短路

此时，若A走到达E后走路径2到达F会与路径C-D有更多交点

即走绿色这条路<-也是最短路

所以最优方案中特殊奖励的点一定是连续的。

于是我们可以用同时在A-B最短路上和C-D最短路上的有向边构成一个新图

那么新图一定是一个有向无环图，但不一定联通

最后对新图拓扑排序dp最长链即可

代码（给你们演示一下堆优化的SPFA）

#include<queue>

#include<cstdio>

#include<cstring>

#include<algorithm>

using namespace std;

const int maxn=;

const int maxm=;

struct node{int id;long long d;};

bool operator <(node a,node b){return a.d>b.d;}

int info[maxn],nx[maxm<<],v[maxm<<],w[maxm<<],re[maxm<<];long long d[][maxn];

int n,m,en,ans,ecnt,p[],dp[maxn],vis[maxn],inp[maxn],top[maxn],lef[maxm<<];priority_queue<node>q;

void add(int u1,int v1,int w1,int r){nx[++ecnt]=info[u1];info[u1]=ecnt;v[ecnt]=v1;w[ecnt]=w1;re[ecnt]=r;}

void SPFA(int k)

{

    memset(vis,,sizeof vis);

    d[k][p[k]]=;q.push((node){p[k],});

    while(!q.empty())

    {

        node nw=q.top();q.pop();

        if(vis[nw.id]==)continue;vis[nw.id]=;d[k][nw.id]=nw.d;

        for(int i=info[nw.id];i;i=nx[i])

        if(!((k%)^re[i])&&d[k][v[i]]>nw.d+w[i])q.push((node){v[i],nw.d+w[i]});

    }

}

int main()

{

    scanf("%d%d",&n,&m);memset(d,0x7f,sizeof d);

    for(int i=,u1,v1,w1;i<=m;i++)scanf("%d%d%d",&u1,&v1,&w1),add(u1,v1,w1,),add(v1,u1,w1,);

    for(int k=;k<=;k++)scanf("%d",&p[k]),SPFA(k);

    if(d[][p[]]==d[][]||d[][p[]]==d[][]){puts("-1");return ;}

    for(int x=;x<=n;x++)

    {

        for(int e=info[x];e;e=nx[e])

        if(re[e]&&d[][x]+w[e]+d[][v[e]]==d[][p[]]&&d[][x]+w[e]+d[][v[e]]==d[][p[]])

        lef[e]=,inp[v[e]]++;

    }

    for(int x=;x<=n;x++)if(inp[x]==)top[++en]=x,dp[x]=;

    for(int i=;i<=en;i++)

    {

        int nw=top[i];

        for(int e=info[nw];e;e=nx[e])

            if(lef[e])

            {

                dp[v[e]]=max(dp[v[e]],dp[nw]+);

                ans=max(ans,dp[v[e]]);inp[v[e]]--;

                if(!inp[v[e]])top[++en]=v[e];

            }

    }

    printf("%d\n",ans);

}