bzoj3040最短路 road dij+heap

最新推荐文章于 2019-07-13 22:50:33 发布

EMber _

最新推荐文章于 2019-07-13 22:50:33 发布

阅读量641

点赞数

CC 4.0 BY-SA版权

分类专栏： bzoj 最短路堆

本文链接：https://blog.youkuaiyun.com/qq_35866453/article/details/70477401

bzoj 同时被 3 个专栏收录

452 篇文章

订阅专栏

最短路

29 篇文章

订阅专栏

堆

11 篇文章

订阅专栏

本文记录了一次使用Dijkstra算法解决图论问题时遇到的空间限制挑战。通过逐步尝试不同类型的优先队列实现，包括标准库配对堆、手写斐波那契堆等，并最终采用HZWER的STL配对堆成功解决问题。文章详细描述了优化过程中的心得和技术细节。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

tmd这题的空间卡的我蛋都碎了。。
1e6的点,1e7的边，我开3个1e7的数组就MLE？~~你tm有本事再开小点？~~
。。经历了从优先队列——>stl配对堆——>斐波那契堆（手写）——>手写配对堆的过程。。tmd一个dij我写了一晚上我也是醉了。。
最后还是得用hzwer的STL配对堆，，~~空间能开的再走心一点？~~

#include<iostream>
#include<cstdio>
#include<cstring>
#include<ext/pb_ds/priority_queue.hpp>
#define ll long long
#define pa pair<ll,int>
#define llinf 9000000000000000000LL
using namespace std;
using namespace __gnu_pbds;
typedef __gnu_pbds::priority_queue<pa,greater<pa>,pairing_heap_tag > heap;
int n,m,cnt,last[1000005];
int T,rxa,rxc,rya,ryc,rp;
heap::point_iterator id[1000005];
int x,y,z;
ll dis[1000005];
struct data{int to,next,v;}e[10000005];
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void insert(int u,int v,int w)
{
    e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;e[cnt].v=w;
}
void dijkstra()
{
    heap q;
    for(int i=1;i<=n;i++)dis[i]=llinf;
    dis[1]=0;id[1]=q.push(make_pair(0,1));
    while(!q.empty())
    {
        int now=q.top().second;q.pop();
        for(int i=last[now];i;i=e[i].next)
            if(e[i].v+dis[now]<dis[e[i].to])
            {
                dis[e[i].to]=e[i].v+dis[now];
                if(id[e[i].to]!=0)
                    q.modify(id[e[i].to],make_pair(dis[e[i].to],e[i].to));
                else id[e[i].to]=q.push(make_pair(dis[e[i].to],e[i].to));
            }
    }
}
int main()
{
    n=read();m=read();
    T=read();rxa=read();rxc=read();rya=read();ryc=read();rp=read();
    int a,b;
    for(int i=1;i<=T;i++)
    {
        x=((ll)x*rxa+rxc)%rp;
        y=((ll)y*rya+ryc)%rp;
        a=min(x%n+1,y%n+1);
        b=max(y%n+1,y%n+1);
        insert(a,b,100000000-100*a);
    }
    for(int i=1;i<=m-T;i++)
    {
        x=read(),y=read(),z=read();
        insert(x,y,z);
    }
    dijkstra();
    printf("%lld",dis[n]);
    return 0;
}