题目大意
(n)((nleq2*10^5))个点,(m)((mleq4*10^5))条边的图,每条边有海拔(a_i(a_ileq10^9))、长度(l_i(l_ileq10^4)),定义两点(a,b)距离为从(a)走到(b)至少要走的长度之和
(q)组询问,强制在线,每次给出(v,p),表示询问不走(a_ileq p)的边,从(v)出发能走到的与(1)号点距离最近的点到(1)号点的距离
题解
预处理每个点到(1)号点的距离
每次询问相当于在问删掉所有(a_ileq p)的边,点(v)所在连通块中到(1)号点最小的距离
发现只考虑原图的最大生成树上的边,不会改变连通性
点(v)所在连通块之所以到不了别的点集,是因为它们之间在最大生成树上的路径中有一条边(a_ileq p)
考虑kruskal的过程,相当于有一次是用一条(a_ileq p)的边合并了点(v)所在连通块与其他点集
这样就可以以这种方法建一新棵树:一开始有(n)个点,没有边,kruskal中每合并两个点集,就新建一个表示当前边的点,并且将两个点集新树中的根变成新建点的儿子
新树中一个点的子树表示这个点对应的原图的最大生成树中一条边kruskal时合并的两个点集,也就是说,这个子树中任意两点的路径中不会出现边权小于该边的边
所以每次询问在新树中找(v)的深度最小的(a_i> p)的祖先
这个新树也叫kruskal重构树
代码
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<iomanip>
#include<iostream>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>
#define rep(i,x,y) for(register int i=(x);i<=(y);++i)
#define dwn(i,x,y) for(register int i=(x);i>=(y);--i)
#define view(u,k) for(int k=fir[u];~k;k=nxt[k])
#define maxn 200010
#define maxm 800010
#define pii pair<int,int>
#define fi first
#define se second
#define mp make_pair
#define LL long long
using namespace std;
int read()
{
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)&&ch!='-')ch=getchar();
if(ch=='-')f=-1,ch=getchar();
while(isdigit(ch))x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
return x*f;
}
void write(int x)
{
if(x==0){putchar('0'),putchar('
');return;}
int f=0;char ch[20];
if(x<0)putchar('-'),x=-x;
while(x)ch[++f]=x%10+'0',x/=10;
while(f)putchar(ch[f--]);
putchar('
');
return;
}
int t,n,m,fir[maxn],nxt[maxm],v[maxm],w[maxm],fa[maxn],anc[maxn<<1][20],val[maxn<<1],mind[maxn<<1],ans;
int cnte,cntnd,dis[maxn],vis[maxn],q,k,s;
struct edge{int u,v,w;}e[maxm>>1];
void ade(int u1,int v1,int w1){v[cnte]=v1,w[cnte]=w1,nxt[cnte]=fir[u1],fir[u1]=cnte++;}
priority_queue<pii >Q;
bool cmp(edge x,edge y){return x.w>y.w;}
void reset()
{
rep(i,1,n)fa[i]=-i,fir[i]=-1,dis[i]=2147483647,vis[i]=0;
rep(i,1,(n<<1)){val[i]=0,mind[i]=2147483647;rep(j,0,19)anc[i][j]=0;}ans=0;
cnte=0,cntnd=n;
}
int f(int x){return fa[x]<0?x:fa[x]=f(fa[x]);}
int main()
{
t=read();
while(t--)
{
n=read(),m=read();
reset();
rep(i,1,m){int x=read(),y=read(),l=read(),a=read();ade(x,y,l),ade(y,x,l),e[i].u=x,e[i].v=y,e[i].w=a;}
sort(e+1,e+m+1,cmp);
dis[1]=0;Q.push(mp(0,1));
while(!Q.empty())
{
int u=Q.top().se;Q.pop();
if(vis[u])continue;vis[u]=1;
view(u,k)if(dis[v[k]]>dis[u]+w[k])
{
dis[v[k]]=dis[u]+w[k];
if(!vis[v[k]])Q.push(mp(-dis[v[k]],v[k]));
}
}
rep(i,1,n)mind[i]=dis[i];
rep(i,1,m)
{
int x=f(e[i].u),y=f(e[i].v);
if(x!=y)
{
cntnd++,val[cntnd]=e[i].w,anc[-fa[x]][0]=anc[-fa[y]][0]=cntnd,mind[cntnd]=min(mind[-fa[x]],mind[-fa[y]]);
fa[x]=y,fa[y]=-cntnd;
}
}
dwn(i,cntnd,1){rep(j,1,19)anc[i][j]=anc[anc[i][j-1]][j-1];}
q=read(),k=read(),s=read();
while(q--)
{
int u=read(),p=read(),tu;
u=(u+k*ans-1)%n+1,p=((LL)p+(LL)k*ans)%(LL)(s+1),tu=u;
dwn(i,19,0)if(anc[tu][i]&&val[anc[tu][i]]>p)tu=anc[tu][i];
write(ans=mind[tu]);
}
}
return 0;
}
/*
2
4 3
1 2 50 1
2 3 100 2
3 4 50 1
5 0 2
3 0
2 1
4 1
3 1
3 2
4 3
1 2 50 1
2 3 100 2
3 4 50 1
5 0 2
3 0
2 1
4 1
3 1
3 2
*/