1. 程式人生 > >[BZOJ]5415: [Noi2018]歸程 Kruskal重構樹

[BZOJ]5415: [Noi2018]歸程 Kruskal重構樹

題解:

複習一下這個東西,做道題爽爽。
關於Kruskal重構樹的知識可以看popoqqq的部落格。這道題目就是個裸題,利用重構樹可以求出 v v 可以免費到哪些點,這些點在一個子樹內,然後看這些點與 1 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); } } }