Luogu4768 NOI2018归程（最短路径+kruskal重构树）

　　按海拔从大到小合并建出kruskal重构树，这样就能知道开车能到达哪些点，对这些点到1的最短路取min即可。最难的部分在于多组数据的初始化和数组大小的设置。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cassert>
using namespace std;
#define ll long long
#define N 200010
#define M 400010
#define inf 2000000010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<''||c>'')) c=getchar();return c;}
int gcd(int n,int m){return m==?n:gcd(m,n%m);}
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
int T,n,m,p[N],d[N],fa[N<<],t,cnt;
bool flag[N];
struct data{int to,nxt,len,h;
}edge[M<<];
void addedge(int x,int y,int z,int h){t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].len=z,edge[t].h=h,p[x]=t;}
struct data2
{
    int x,d;
    bool operator <(const data2&a) const
    {
        return d>a.d;
    }
};
struct data3
{
    int x,y,z;
    bool operator <(const data3&a) const
    {
        return z>a.z;
    }
}e[M];
priority_queue<data2> q;
int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
void dijkstra()
{
    while (!q.empty()) q.pop();q.push((data2){,});
    for (int i=;i<=n;i++) d[i]=inf;d[]=;
    memset(flag,,sizeof(flag));
    for (;;)
    {
        while (!q.empty()&&flag[q.top().x]) q.pop();
        if (q.empty()) break;
        data2 x=q.top();q.pop();
        flag[x.x]=;
        for (int i=p[x.x];i;i=edge[i].nxt)
        if (x.d+edge[i].len<d[edge[i].to])
        {
            d[edge[i].to]=x.d+edge[i].len;
            q.push((data2){edge[i].to,d[edge[i].to]});
        }
    }
}
namespace kruskal_tree
{
    int p[N<<],t,fa[N<<][],val[N<<],h[N<<];
    struct data{int to,nxt,len;}edge[N<<];
    void addedge(int x,int y){t++;edge[t].to=y,edge[t].nxt=p[x],p[x]=t;}
    void dfs(int k)
    {
        val[k]=k<=n?d[k]:inf;
        for (int i=p[k];i;i=edge[i].nxt)
        {
            fa[edge[i].to][]=k;
            dfs(edge[i].to);
            val[k]=min(val[k],val[edge[i].to]);
        }
    }
    void build()
    {
        fa[cnt][]=cnt;dfs(cnt);
        for (int j=;j<;j++)
            for (int i=;i<=cnt;i++)
            fa[i][j]=fa[fa[i][j-]][j-];
    }
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("return.in","r",stdin);
    freopen("return.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    T=read();
    while (T--)
    {
        n=read(),m=read();
        for (int i=;i<=n;i++) p[i]=;t=;
        for (int i=;i<=m;i++)
        {
            int x=read(),y=read(),z=read(),h=read();
            e[i].x=x,e[i].y=y,e[i].z=h;
            addedge(x,y,z,h),addedge(y,x,z,h);
        }
        dijkstra();
        sort(e+,e+m+);
        for (int i=;i<=n;i++) fa[i]=i,kruskal_tree::h[i]=inf,kruskal_tree::p[i]=;cnt=n;kruskal_tree::t=;
        for (int i=;i<=m;i++)
        {
            int p=find(e[i].x),q=find(e[i].y);
            if (p!=q)
            {
                cnt++;fa[cnt]=fa[p]=fa[q]=cnt;kruskal_tree::p[cnt]=;
                kruskal_tree::addedge(cnt,p),kruskal_tree::addedge(cnt,q);
                kruskal_tree::h[cnt]=e[i].z;
            }
        }
        kruskal_tree::build();
        int Q=read(),K=read(),S=read(),ans=;
        while (Q--)
        {
            int x=(read()+K*ans-)%n+,y=(read()+K*ans)%(S+);
            for (int j=;~j;j--) if (kruskal_tree::h[kruskal_tree::fa[x][j]]>y) x=kruskal_tree::fa[x][j];
            printf("%d\n",ans=kruskal_tree::val[x]);
        }
    }
    return ;
}