[BZOJ]5415: [Noi2018]歸程 Kruskal重構樹
阿新 • • 發佈:2018-12-24
題解:
複習一下這個東西,做道題爽爽。
關於Kruskal重構樹的知識可以看popoqqq的部落格。這道題目就是個裸題,利用重構樹可以求出
可以免費到哪些點,這些點在一個子樹內,然後看這些點與
的最短路就行了,十分簡單。
程式碼:
#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);
}
}
}