洛谷P4768 [NOI2018]归程 [可持久化并查集，Dijkstra]

　　题目传送门

归程

格式难调，题面就不放了。

---

　　分析：

　　之前同步赛的时候反正就一脸懵逼，然后场场暴力大战，现在呢，还是不会$Kruskal$重构树，于是就拿可持久化并查集做。

　　但是之前做可持久化并查集的时候感觉掌握的并不熟，还是需要参照别人的题解，不过至少现在对可持久化的理解更深了一步，而且终于这题给调对了。

　　Code：

//It is made by HolseLee on 23rd Aug 2018
//Luogu.org 4768
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define Min(a,b) (a)<(b)?(a):(b)
#define Max(a,b) (a)>(b)?(a):(b)
using namespace std;

const int N=2e5+,M=8e5+,S=2e7+;
typedef unsigned int ui;
ui T,n,m,Q,K,s,lastans,L,cnt;
ui head[N],nxt[M],to[M],d[M],a[M],b[M],dis[N];
ui rt[M],lc[S],rc[S],dep[S],fa[S],mn[S];
bool vis[N];
struct Node{
    ui id,val;
    bool operator < (const Node x) const {
        return val>x.val;
    }
    Node(ui x=,ui y=){
        id=x,val=y;
    }
};
struct Edge{
    ui u,v,h;
    bool operator < (const Edge x) const {
        return h<x.h;
    }
    Edge(ui x=,ui y=,ui z=){
        u=x,v=y,h=z;
    }
}e[M];
priority_queue<Node> t; 

inline ui read()
{
    char ch=getchar();ui num=;
    while(ch<''||ch>'')ch=getchar();
    while(ch>=''&&ch<=''){
        num=(num<<)+(num<<)+(ch^);ch=getchar();
    }
    return num;
}

void build(ui &root,ui l,ui r)
{
    root=++cnt;
    if(l==r){mn[root]=dis[fa[root]=l];return;}
    ui mid=(l+r)>>;
    build(lc[root],l,mid);
    build(rc[root],mid+,r);
}

ui ins(ui *root,ui las,ui tag)
{
    ui l=,r=n,mid;
    while(l!=r){
        *root=++cnt,mid=(l+r)>>;
        if(tag<=mid)r=mid,rc[*root]=rc[las],root=lc+*root,las=lc[las];
        else l=mid+,lc[*root]=lc[las],root=rc+*root,las=rc[las];
    }
    return *root=++cnt;
}

ui find(ui root,ui tag)
{
    ui now,l,r,mid;
    while(){
        now=root,l=,r=n;
        while(l!=r){
            mid=(l+r)>>;
            if(tag<=mid)r=mid,now=lc[now];
            else l=mid+,now=rc[now];
        }
        if(tag==fa[now])break;
        tag=fa[now];
    }
    return now;
}

int main()
{
    T=read();
    ui i,j;
    while(T--){
        n=read();m=read();
        ui x,y,w,cnte=;cnt=;
        for(i=;i<=m;++i){
            x=read();y=read();
            to[++cnte]=y,nxt[cnte]=head[x],head[x]=cnte;
            to[++cnte]=x,nxt[cnte]=head[y],head[y]=cnte;
            d[cnte]=d[cnte-]=read();
            e[i]=Edge(x,y,a[cnte]=a[cnte-]=read());
        }
        memset(dis,-,(n+)<<);
        dis[]=;t.push(Node(,));
        while(!t.empty()){
            Node now=t.top();t.pop();
            if(vis[x=now.id])continue;
            vis[x]=true;
            for(i=head[x];i;i=nxt[i]){
                y=to[i];
                if(dis[y]>dis[x]+d[i])
                t.push(Node(y,dis[y]=dis[x]+d[i]));
            }
        }
        Q=read(),K=read(),s=read();lastans=;
        sort(e+,e+m+);
        for(i=;i<=m;++i)b[i]=e[i].h;
        b[m+]=s+;
        L=unique(b+,b+m+)-b-;
        build(rt[L],,n);
        for(i=L-,j=m;i;--i){
            rt[i]=rt[i+];
            for(;j&&e[j].h==b[i];--j){
                if((x=find(rt[i],e[j].u))==(y=find(rt[i],e[j].v)))continue;
                if(dep[x]>dep[y])swap(x,y);
                fa[ins(&rt[i],rt[i],fa[x])]=fa[y];
                w=ins(&rt[i],rt[i],fa[y]);
                fa[w]=fa[y];mn[w]=Min(mn[x],mn[y]);
                dep[w]=dep[y]+(dep[y]==dep[x]);
            }
        }
        for(;Q;--Q){
            x=(read()+K*lastans-)%n+;
            y=(read()+K*lastans)%(s+);
            printf("%u\n",lastans=mn[find(rt[upper_bound(b+,b+L+,y)-b],x)]);
        }
        memset(vis,,n+);
        memset(head,,(n+)<<);
        memset(rt,,(L+)<<);
        memset(lc,,(cnt+)<<);
        memset(rc,,(cnt+)<<);
        memset(fa,,(cnt+)<<);
        memset(mn,,(cnt+)<<);
        memset(dep,,(cnt+)<<);
    }
    return ;
}