HDU 3191 次短路长度和条数

http://www.cnblogs.com/wally/archive/2013/04/16/3024490.html

http://blog.csdn.net/me4546/article/details/6584448

维护最短路长度d[i][0]和次短路d[i][1]，最短路条数dp[i][0]和次短路dp[i][1]

#include <iostream>
#include <string>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <stack>
#include <queue>
#include <cctype>
#include <vector>
#include <iterator>
#include <set>
#include <map>
#include <sstream>
using namespace std;
 
#define mem(a,b) memset(a,b,sizeof(a))
#define pf printf
#define sf scanf
#define spf sprintf
#define pb push_back
#define debug printf("!\n")
#define MAXN 5010
#define MAX(a,b) a>b?a:b
#define blank pf("\n")
#define LL long long
#define ALL(x) x.begin(),x.end()
#define INS(x) inserter(x,x.begin())
#define pqueue priority_queue
#define INF 0x3f3f3f3f
 
struct Edge
{
    int v,w;
};
vector<Edge> g[];
struct node
{
    int v,dist;
    int mark;//标记，0为最短路，1为次短路；
    bool operator < (const node &p) const
    {
        if(p.dist != dist)
            return p.dist<dist;
 
        return p.v<v;//这儿如果不按顶点的大小排序，就wa了。
    }
};
 
int n,m,s,e;
int dist[][];
int dp[][];
bool vis[][];
//dp[i][0]表示到达点i最短的路有多少条，dp[i][1]表示次短的条数
//dist[i][0]表示到达点i最短路的长度,dist[i][1]表示次短路的长度
/*
用v去松驰u时有四种情况 (设当前dist[v][cas])
情况1：dist[v][cas]+w(v,u)<dist[u][0]，找到一个更短的距离，则把原来最短的距离作为次短的距离，同时更新最短的.即
dist[u][1]=dist[u][0]
dist[u][0]=dist[v][cas]+w(v,u);
dp[u][1]=dp[u][0]
dp[u][0]=dp[v][cas]，
把node(dist[u][0],u,0)和node(dist[u][1],u,1)放入队列
情况2：dist[v][cas]+w(v,u)==dist[u][0]，找到一条新的相同距离的最短路，则dp[u][0]+=dp[v][cas],其他不用更新,也不入队
情况3：dist[v][cas]+w(v,u)<dist[u][1]，不可以更新最短距离，但可以更新次短的，则更新dist[u][1]和dp[u][1]
dist[u][1]=dist[v][cas]+w(v,u);
dp[u][1]=dp[v][cas];
把node(dist[u][1],u,1)放入队列
情况4：dist[v][cas]+w(v,u)==dist[u][1] 找到一条新的相同距离的次短路,则dp[u][1]+=dp[v][cas],其他不更新。
*/
 
void dijkstra(int start,int end)
{
    for(int i=;i<n;i++)
    {
        dist[i][]=dist[i][]=INF;
    }
    memset(dp,,sizeof(dp));
    memset(vis,false,sizeof(vis));
    priority_queue<node> Q;
    node p,q;
    dist[start][]=;
    dp[start][]=;
    p.dist=,p.mark=,p.v=start;
    Q.push(p);
    while(!Q.empty())
    {
        p=Q.top();
        Q.pop();
        if(vis[p.v][p.mark]) continue;
        //if(dist[p.v][p.mark]!=p.dist)continue;
        vis[p.v][p.mark]= true;
        for(int i=;i<g[p.v].size();i++)
        {
            int v=g[p.v][i].v;
            int w=g[p.v][i].w;
            if(!vis[v][] && p.dist+w <dist[v][])
            {
                //可能为次短路
                if(dist[v][]!=INF)
                {
                    q.v=v,q.dist=dist[v][],q.mark=;
                    dist[v][]=dist[v][];
                    dp[v][]=dp[v][];
                    Q.push(q);
                }
                dist[v][]=p.dist+w;
                dp[v][]=dp[p.v][p.mark];
                q.v=v,q.dist=dist[v][],q.mark=;
                Q.push(q);
            }
            else if(!vis[v][] && p.dist+w==dist[v][])
            {
                dp[v][]+=dp[p.v][p.mark];
            }
            else if(!vis[v][] && p.dist+w<dist[v][])
            {
                dist[v][]=p.dist+w;
                dp[v][]=dp[p.v][p.mark];
                q.dist=dist[v][],q.v=v,q.mark=;
                Q.push(q);
            }
            else if(!vis[v][]&&p.dist+w==dist[v][])
            {
                dp[v][]+=dp[p.v][p.mark];
            }
        }
    }
}
 
int main(){
    while(~scanf("%d%d%d%d",&n,&m,&s,&e))
    {
        for(int i=;i<n;i++)g[i].clear();
        for(int i=;i<=m;i++)
        {
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            Edge p;
            p.v=v,p.w=w;
            g[u].push_back(p);
        }
        dijkstra(s,e);
        printf("%d %d\n",dist[e][],dp[e][]);
    }
    return ;
}