【noi2018】归程

最新推荐文章于 2021-12-03 09:32:55 发布

anweiyun4308

最新推荐文章于 2021-12-03 09:32:55 发布

阅读量191

点赞数

原文链接：http://www.cnblogs.com/rlddd/p/9657710.html

版权

题目描述

本题的故事发生在魔力之都，在这里我们将为你介绍一些必要的设定。魔力之都可以抽象成一个

Yazid 是一名来自魔力之都的OIer，刚参加完ION2018 的他将踏上归程，回到他温暖的家。 Yazid 的家恰好在魔力之都的

车会在新的出发点被准备好。
Yazid 不能利用之前在某处停放的车。 Yazid 非常讨厌在雨天步行，因此他希望在完成回家这一目标的同时，最小化他步行经过的边的总长度。请你帮助 Yazid 进行计算。本题的部分测试点将强制在线，具体细节请见【输入格式】和【子任务】。

输入

单个测试点中包含多组数据。输入的第一行为一个非负整数

接下来依次描述每组数据，对于每组数据：

第一行

接下来

接下来一行

接下来

对于输入中的每一行，如果该行包含多个数，则用单个空格将它们隔开。

输出

依次输出各组数据的答案。对于每组数据：

输出

样例输入

样例输出

题解

kruskal重构树终于有用了。先dijstra求一遍最短路，以海拔建最大生成树，然后kruskal重构树。很容易想到，我们求一遍 x 到1的边中海拔最小的h，如果h大于水位线，那么一定是可以开车回家的；如果h小于等于水位线，那么在到这条边之前就得下车走路。如何找到该下车的点？重构树上的每个点都对应在原生成树上的一段区间，且重构树的深度越小，对应的区间越大。那么找到重构树上的大于水位线的深度最小的点，答案就等于它所对应的生成树上的点距离1号点的最短路的最小值。做法：最短路+kruskal重构树+树上倍增。

#include<queue>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long

const int maxn=800010+50;
const int maxm=2*maxn;

int fir[maxn],to[maxm],nex[maxm],wi[maxm],ecnt,vis[maxn];
int val[maxn],fat[maxn],d[maxn],mi[maxn],anc[maxn][23];
int T,cnt,n,m,u,v,l,aa,Q,K,S,lastans,v0,p0,vv,pp;

priority_queue< pair<ll,int> > q;

struct node{int x,y,z;}a[maxm];
int cmp(const node &a,const node &b){return a.z>b.z;}

#define mem(a,b) memset(a,b,sizeof(a))
void init(){
    mem(fir,0);mem(to,0);mem(nex,0);mem(wi,0);ecnt=0;lastans=0;
    mem(val,0);mem(fat,0);mem(mi,0);mem(a,0);mem(anc,0);
}

void init2(){
    mem(fir,0);mem(to,0);mem(nex,0);mem(wi,0);ecnt=0;
}

void add_edge(int u,int v,int w){
    nex[++ecnt]=fir[u];fir[u]=ecnt;to[ecnt]=v;wi[ecnt]=w;
}

int father(int x){
    if(x!=fat[x]) fat[x]=father(fat[x]);
    return fat[x];
}

template<typename T>void read(T& aa){
    char cc; ll ff;aa=0;cc=getchar();ff=1;
    while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
    if(cc=='-') ff=-1,cc=getchar();
    while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
    aa*=ff;
}

void dfs(int u){
    mi[u]=d[u];
    for(int e=fir[u];e;e=nex[e]){
        int v=to[e];
        anc[v][0]=u;
        dfs(v);
        mi[u]=min(mi[u],mi[v]);
    }
}

void dij(){
    memset(vis,0,sizeof(vis));
    memset(d,111,sizeof(d)); d[1]=0;
    q.push(make_pair(0,1));
    while(!q.empty()){
        int u=q.top().second; q.pop();
        if(vis[u]) continue;
        vis[u]=1;
        for(int e=fir[u];e;e=nex[e]){
            int v=to[e];
            if(d[u]+wi[e]<d[v]){
                d[v]=d[u]+wi[e];
                q.push(make_pair(-d[v],v));
            }
        }
    }
}

void kruskal(){
    sort(a+1,a+1+m,cmp);
    for(int i=1;i<=m;i++){
        int fa=father(a[i].x),fb=father(a[i].y);
        if(fa!=fb){
            val[++cnt]=a[i].z;
            fat[cnt]=fat[fa]=fat[fb]=cnt;
            add_edge(cnt,fa,0);
            add_edge(cnt,fb,0);
        }
    }
    dfs(cnt);
}

int main(){
    read(T);
    while(T--){
        init();
        read(n),read(m);cnt=n;
        for(int i=1;i<=n;i++) fat[i]=i;
        for(int i=1;i<=m;i++){
            read(u),read(v),read(l),read(aa);
            a[i].x=u;a[i].y=v;a[i].z=aa;
            add_edge(u,v,l);
            add_edge(v,u,l);
        }
        dij();
        init2();
        kruskal();
        for(int i=1;(1<<i)<=cnt;i++)
        for(int u=1;u<=cnt;u++)
        anc[u][i]=anc[anc[u][i-1]][i-1];
        read(Q),read(K),read(S);
        while(Q--){
            read(v0),read(p0);
            vv=(v0+K*lastans-1)%n+1;
            pp=(p0+K*lastans)%(S+1);
            for(int i=22;i>=0;i--)
            if(anc[vv][i]&&val[anc[vv][i]]>pp) vv=anc[vv][i];
            printf("%d\n",mi[vv]);lastans=mi[vv];
        }
    }
    return 0;
}