【NOI2018DAY1T1】【洛谷P4768】—归程[Return](Kruscal重构树)
重构树简单题
考虑到低于水位线的边都会被淹没
那显然我们做最大生成树后树上路径就是最优的情况
多次询问考虑重构树
显然对于当前一个要求的的水位和点
我们倍增找到第一个刚好大于水位线的点
那显然答案就是该点的子树中所有叶子节点到的最短的一条
直接最短路后树形预处理一下就可以了
复杂度
#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))<<'\n';
}
}
}