【NOI2018DAY1T1】【洛谷P4768】—归程[Return]（Kruscal重构树）

传送门

重构树简单题

和$Peaks$相似，还要简单一些

考虑到低于水位线的边都会被淹没
那显然我们做最大生成树后树上路径就是最优的情况

多次询问考虑$Kruscal$重构树

显然对于当前一个要求的的水位和点

我们倍增找到第一个刚好大于水位线的点

那显然答案就是该点的子树中所有叶子节点到$1$的最短的一条

直接最短路后树形$dp$预处理一下就可以了

复杂度$O(T*(mlogm+mlogn+nlogn+qlogn))$

#include<bits/stdc++.h>
using namespace std;
inline int read(){
	char ch=getchar();
	int res=0,f=1;
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch))res=(res<<3)+(res<<1)+(ch^48),ch=getchar();
	return res*f;
}
const int N=200005;
const int M=400005;
const int Log=22;
int n,m,q,s,k,a[N<<1],ans;
struct edge{
	int u,v,w;
}e[M];
namespace Kruscal{
	int tot,fa[N<<1],val[N<<1],f[N<<1][Log],lc[N<<1],rc[N<<1],cnt;
	int dis[N],g[N<<1];
	inline bool comp(const edge &a,const edge &b){
		return a.w>b.w;
	}
	int find(int x){
		return fa[x]==x?x:fa[x]=find(fa[x]);
	}
	inline void kruscal(){
		sort(e+1,e+m+1,comp);
		for(int i=1;i<=m;i++){
			int f1=find(e[i].u),f2=find(e[i].v);
			if(f1!=f2){
				fa[f1]=fa[f2]=fa[++tot]=tot;
				val[tot]=e[i].w;
				f[f1][0]=f[f2][0]=tot;
				lc[tot]=f1,rc[tot]=f2;
			}
		}
	}
	void dfs(int u){
		for(int i=1;i<=20;i++)
			f[u][i]=f[f[u][i-1]][i-1];
		if(u<=n){
			g[u]=dis[u];
			return;
		}
		dfs(lc[u]),dfs(rc[u]);
		g[u]=min(g[lc[u]],g[rc[u]]);
	}
	inline int get(int u,int k){
		for(int i=20;i>=0;i--){
			if(f[u][i]&&val[f[u][i]]>k)
				u=f[u][i];
		}
		return g[u];
	}
}
using namespace Kruscal;
namespace Pre_Operator{
	priority_queue<pair<int,int>,vector<pair<int,int> >,greater<pair<int,int> > > que;
	#define mp make_pair
	int vis[N],adj[N<<1],nxt[N<<2],to[N<<2],cost[N<<2];
	inline void addedge(int u,int v,int w){
		nxt[++cnt]=adj[u],adj[u]=cnt,to[cnt]=v,cost[cnt]=w;
	}
	inline void dijkstra(){
		dis[1]=0;que.push(mp(0,1));
		while(!que.empty()){
			int u=que.top().second;que.pop();
			if(vis[u])continue;vis[u]=1;
			for(int e=adj[u];e;e=nxt[e]){
				int v=to[e];
				if(dis[u]+cost[e]<dis[v]){
					dis[v]=dis[u]+cost[e];
					que.push(mp(dis[v],v));
				}
			}
		}
	}
	inline void init(){
		memset(dis,127,sizeof(dis));
		memset(adj,0,sizeof(adj));
		memset(vis,0,sizeof(vis));
		tot=cnt=ans=0;
	}
}
using namespace Pre_Operator;
int main(){
	int T=read();
	while(T--){
		init();
		n=tot=read(),m=read();
		for(int i=1;i<=n;i++)fa[i]=i;
		for(int i=1;i<=m;i++){
			int u=read(),v=read(),w=read(),h=read();
			e[i].u=u,e[i].v=v,e[i].w=h;
			addedge(u,v,w),addedge(v,u,w);
		}
		dijkstra();
		kruscal();
		dfs(tot);
		q=read(),k=read(),s=read();
		for(int i=1;i<=q;i++){
			int v=read(),p=read();
			v=(v+k*ans-1)%n+1,p=(p+k*ans)%(s+1);
			cout<<(ans=get(v,p))<<'
';
		}
	}
}