[BZOJ]5415: [Noi2018]归程 Kruskal重构树

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/baidu_36797646/article/details/85231453

题解：

复习一下这个东西，做道题爽爽。
关于Kruskal重构树的知识可以看popoqqq的博客。这道题目就是个裸题，利用重构树可以求出 $v$可以免费到哪些点，这些点在一个子树内，然后看这些点与 $1$的最短路就行了，十分简单。

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=200010;
const int Maxm=400010;
const int inf=2147483647;
int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
	return x*f;
}
int n,m,Q,K,S,ans,tot;
int val[Maxn<<1],lc[Maxn<<1],rc[Maxn<<1],dep[Maxn<<1],fa[Maxn<<1][19],minf[Maxn<<1];
struct Edge{int x,y,d,next;}e[Maxm<<1],E[Maxm];
bool cmp(Edge a,Edge b){return a.d>b.d;}
int rt[Maxn];
int findrt(int x){return((x==rt[x])?x:rt[x]=findrt(rt[x]));}
int last[Maxn],len;
void ins(int x,int y,int d)
{
	int t=++len;
	e[t].y=y;e[t].d=d;e[t].next=last[x];last[x]=t;
}
struct Node
{
	int x,t;
	Node(int _x,int _t){x=_x,t=_t;}
};
bool operator < (Node a,Node b){return a.t>b.t;}
priority_queue<Node>q;
int f[Maxn];bool vis[Maxn];
void dijkstra(int st)
{
	memset(f,63,sizeof(f));f[st]=0;
	memset(vis,false,sizeof(vis));
	q.push(Node(st,0));
	while(!q.empty())
	{
		Node tmp=q.top();q.pop();
		int x=tmp.x;
		if(vis[x])continue;
		vis[x]=true;
		for(int i=last[x];i;i=e[i].next)
		{
			int y=e[i].y;
			if(f[x]+e[i].d<f[y])f[y]=f[x]+e[i].d,q.push(Node(y,f[y]));
		}
	}
}
void build()
{
	sort(E+1,E+1+m,cmp);
	for(int i=1;i<=n;i++)rt[i]=i,val[i]=inf;tot=n;
	for(int i=1;i<=m;i++)
	{
		int x=E[i].x,y=E[i].y;
		int fx=findrt(x),fy=findrt(y);
		if(fx==fy)continue;
		rt[fx]=rt[fy]=++tot;rt[tot]=tot;
		val[tot]=E[i].d;lc[tot]=fx,rc[tot]=fy;
	}
}
void dfs(int x,int ff)
{
	if(x<=n)minf[x]=f[x];
	fa[x][0]=ff;dep[x]=dep[ff]+1;
	for(int i=1;(1<<i)<=dep[x];i++)fa[x][i]=fa[fa[x][i-1]][i-1];
	if(x>n)dfs(lc[x],x),dfs(rc[x],x),minf[x]=min(minf[lc[x]],minf[rc[x]]);
}
int get(int x,int lim)
{
	for(int i=18;i>=0;i--)
	if((1<<i)<=dep[x]&&val[fa[x][i]]>lim)x=fa[x][i];
	return x;
}
int main()
{
	int T=read();
	while(T--)
	{
		memset(last,0,sizeof(last));len=0;
		n=read(),m=read();
		for(int i=1;i<=m;i++)
		{
			int x=read(),y=read(),d=read(),h=read();
			E[i].x=x;E[i].y=y;E[i].d=h;
			ins(x,y,d),ins(y,x,d);
		}
		Q=read(),K=read(),S=read(),ans=0;
		dijkstra(1);
		build();
		dep[0]=-1;dfs(tot,0);
		while(Q--)
		{
			int v=read(),p=read();
			v=(v+K*ans%n-1)%n+1,p=(p+K*ans)%(S+1);
			ans=minf[get(v,p)];
			printf("%d\n",ans);
		}
	}
}