[Luogu1342] 请柬 - 最短路模板

Description

在电视时代,没有多少人观看戏剧表演。Malidinesia古董喜剧演员意识到这一事实，他们想宣传剧院,尤其是古色古香的喜剧片。他们已经打印请帖和所有必要的信息和计划。许多学生被雇来分发这些请柬。每个学生志愿者被指定一个确切的公共汽车站，他或她将留在那里一整天,邀请人们参与。

这里的公交系统是非常特殊的:所有的线路都是单向的，连接两个站点。公共汽车离开起始点，到达目的地之后又空车返回起始点。学生每天早上从总部出发，乘公交车到一个预定的站点邀请乘客。每个站点都被安排了一名学生。在一天结束的时候,所有的学生都回到总部。现在需要知道的是，学生所需的公交费用的总和最小是多少。

Input&Output

Input

第1行有两个整数n、m(1<=n,m<=1000000)，n是站点的个数，m是线路的个数。

然后有m行，每行描述一个线路，包括3个整数，起始点，目的地和价格。

总部在第1个站点，价钱都是整数，且小于1000000000。

Output

输出一行，表示最小费用。

Sample

Input

Output

Solution

本题是一道裸的最短路。需要特殊处理的问题在于，由于路线都是单向的，我们在求返程最短路时需要建一张反图。不得不说数据范围比较毒瘤，堆优化dijkstra可以稳定在O(nlogn)，相较于O(玄学)的SPFA还是稳一些，另外建议读入优化。

代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cctype>
#include<cstring>
#define maxn 1000005
#define INF_ 214748000000
using namespace std;
typedef unsigned long long ull;
struct edge{
    int to,nxt,v;
}e[maxn],re[maxn];
int n,m,edgenum,redgenum,lnk[maxn],rlnk[maxn],u,v,w;
void add(int bgn,int end,int val)
{
    e[++edgenum].to=end;
    e[edgenum].v=val;
    e[edgenum].nxt=lnk[bgn];
    lnk[bgn]=edgenum;
}
void radd(int bgn,int end,int val)
{
    re[++redgenum].to=end;
    re[redgenum].v=val;
    re[redgenum].nxt=rlnk[bgn];
    rlnk[bgn]=redgenum;
}
inline int rd()
{
    int x=0;char c=getchar();
    while(!isdigit(c))c=getchar();
    while(isdigit(c))x=(x<<1)+(x<<3)+(c^48),c=getchar();
    return x;
}
struct node{
    ull dis;
    int id;
    node(ull d,int i):dis(d),id(i){}
    bool operator < (const node &b)const{
        return dis>b.dis;
    }
};
ull dis[maxn];
bool vis[maxn];
void dij()
{
    for(int i=1;i<=n;++i)
        dis[i]=INF_;
    priority_queue<node> q;
    dis[1]=0;
    q.push(node(0,1));
    while(!q.empty()){
        node u=q.top();q.pop();
        if(vis[u.id])continue;
        vis[u.id]=true;
        for(int p=lnk[u.id];p;p=e[p].nxt){
            int y=e[p].to;
            if(dis[y]>dis[u.id]+e[p].v){
                dis[y]=dis[u.id]+e[p].v;
                q.push(node(dis[y],y));
            }
        }
    }
}
void rdij()
{
    for(int i=1;i<=n;++i)
        dis[i]=INF_;
    memset(vis,0,n+1);
    priority_queue<node> q;
    dis[1]=0;
    q.push(node(0,1));
    while(!q.empty()){
        node u=q.top();q.pop();
        if(vis[u.id])continue;
        vis[u.id]=true;
        for(int p=rlnk[u.id];p;p=re[p].nxt){
            int y=re[p].to;
            if(dis[y]>dis[u.id]+re[p].v){
                dis[y]=dis[u.id]+re[p].v;
                q.push(node(dis[y],y));
            }
        }
    }
}
int main()
{
    n=rd();m=rd();
    for(int i=1;i<=m;++i)
        u=rd(),v=rd(),w=rd(),add(u,v,w),radd(v,u,w);
    dij();
    ull sum=0;
    for(int i=2;i<=n;++i)
        sum+=dis[i];
    rdij();
    for(int i=2;i<=n;++i)
        sum+=dis[i];
    printf("%lld\n",sum);
    return 0;
}