[bzoj3040][dijkstra][配对堆]最短路(road)

Description

N个点，M条边的有向图，求点1到点N的最短路（保证存在）。 1<=N<=1000000，1<=M<=10000000

Input

第一行两个整数N、M，表示点数和边数。 第二行六个整数T、rxa、rxc、rya、ryc、rp。

前T条边采用如下方式生成：
1.初始化x=y=z=0。
2.重复以下过程T次： x=(x*rxa+rxc)%rp; y=(y*rya+ryc)%rp; a=min(x%n+1,y%n+1); b=max(y%n+1,y%n+1); 则有一条从a到b的，长度为1e8-100*a的有向边。

后M-T条边采用读入方式： 接下来M-T行每行三个整数x,y,z，表示一条从x到y长度为z的有向边。

1<=x,y<=N，0＜z,rxa,rxc,rya,ryc,rp<2^31

Output

一个整数，表示1~N的最短路。

Sample Input

3 3
0 1 2 3 5 7
1 2 1
1 3 3
2 3 1

Sample Output

2

HINT

【注释】

请采用高效的堆来优化Dijkstra算法。

题解

不死心的我一开始试了一发spfa结果被卡到惨绝人寰
不会Dijkstra啊快学啊。。
其实也很简单。。每次当前扩展距离最少的那个优先扩展，然后扔进堆里。堆顶遇见终点就跳出。
其实是给A*做准备的这题
既然说的这么清楚了，那为什么不能直接A呢？
因为这题会卡优先队列。。
不会啊那就%题解。其实就是背代码。用一个叫啥配对堆的优化这个东西。。
然后。还是卡空间的这题。。我开大了一点点就MLE了
好题

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;
using namespace __gnu_pbds;
typedef long long LL;
typedef __gnu_pbds::priority_queue<pair<LL,int> > heap;
heap::point_iterator id[1000005];
struct node
{
    int x,y,c,next;
}a[10000005];int len,last[1000005];
void ins(int x,int y,int c)
{
    len++;
    a[len].x=x;a[len].y=y;a[len].c=c;
    a[len].next=last[x];last[x]=len;
}
LL d[1000005];
int n,m;
int T,rxa,rxc,rya,ryc,rp;
int tx,ty;
void dijkstra()
{
    heap q;
    memset(d,127,sizeof(d));
    id[1]=q.push(make_pair(0,1));
    d[1]=0;
    while(!q.empty())
    {
        int x=q.top().second;
        q.pop();
        if(x==n)break;
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(d[y]>d[x]+a[k].c)
            {
                d[y]=d[x]+a[k].c;
                if(id[y]!=0)q.modify(id[y],make_pair(-d[y],y));
                else id[y]=q.push(make_pair(-d[y],y));
            }
        }
    }
    printf("%lld\n",d[n]);
}
int main()
{
    scanf("%d%d",&n,&m);
    scanf("%d%d%d%d%d%d",&T,&rxa,&rxc,&rya,&ryc,&rp);
    tx=ty=0;int aa,bb;
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=T;i++)
    {
        tx=((LL)tx*rxa+rxc)%rp;
        ty=((LL)ty*rya+ryc)%rp;
        aa=min(tx%n+1,ty%n+1);
        bb=max(ty%n+1,ty%n+1);
        ins(aa,bb,100000000-100*aa);
    }
    for(int i=1;i<=m-T;i++)
    {
        int x,y,c;
        scanf("%d%d%d",&x,&y,&c);
        ins(x,y,c);
    }

    dijkstra();
    return 0;
}