1. 程式人生 > >洛谷.4768.[NOI2018]歸程(Kruskal重構樹 倍增)

洛谷.4768.[NOI2018]歸程(Kruskal重構樹 倍增)

queue 就是 break 不能 gis fin algo define 容易

題目鏈接

容易想到按高度Kruskal重構樹+預處理到點1的距離dis。
建一棵最大生成樹,如果隨便建的話,如果非樹邊能走,整棵樹都能走答案當然是0...;如果有些樹邊不能走,那麽可走範圍被限制在了某個連通塊。
然而被限制在某個連通塊和圖(還要暴力,難道樹分塊?)沒什麽區別,所以我們可以讓生成樹邊的高度由葉子向上遞減,這樣每次詢問 找到深度最小的可行點後,答案就是其子樹dis最小值(樹形態顯然不會影響什麽)。
就是在Kruskal合並兩個集合時,新建一個節點作為兩集合的代表節點,最低高度mn為這條邊權(當然不會比兩集合中的大),dis為兩集合dis的min。新樹葉子節點即為原所有節點。

昨天一時zz覺得Kruskal不對。。

//4860ms    50.49MB + 4392ms    50.33MB(後4組數據)
#include <queue>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 400000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define mp std::make_pair
#define pr std::pair<int,int>
const int N=2e5+5,M=8e5+5,INF=0x7fffffff;

int n,m,tot,Enum,H[N],nxt[M],to[M],len[M],dis[N],fa[N<<1][19],mn[N<<1],anc[N<<1],Ans[N<<1];
std::priority_queue<pr> q;
char IN[MAXIN],*SS=IN,*TT=IN;
struct Edge
{
    int fr,to,h;
    Edge() {}
    Edge(int fr,int to,int h):fr(fr),to(to),h(h) {}
    bool operator <(const Edge &x)const{
        return h>x.h;
    }
}e[M>>1];

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
inline void AddEdge(int _h,int w,int u,int v)
{
    to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum, len[Enum]=w;
    to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum, len[Enum]=w;
    e[Enum>>1]=Edge(u,v,_h);
}
void Dijkstra()
{
    static bool vis[N];
    memset(dis,0x3f,sizeof dis), memset(vis,0,sizeof vis);
    dis[1]=0, q.push(mp(0,1));
    while(!q.empty())
    {
        int x=q.top().second; q.pop();
        if(vis[x]) continue;
        vis[x]=1;
        for(int v,i=H[x]; i; i=nxt[i])
            if(dis[v=to[i]]>dis[x]+len[i])
                dis[v]=dis[x]+len[i], q.push(mp(-dis[v],v));
    }
}
int Get_fa(int x){
    return x==anc[x]?x:anc[x]=Get_fa(anc[x]);
}
void Kruskal()
{
    for(int i=1; i<=n; ++i) anc[i]=i, Ans[i]=dis[i];
    int m=Enum>>1; std::sort(e+1,e+1+m);
    for(int r1,r2,k=1,i=1; i<=m; ++i)
    {
        if((r1=Get_fa(e[i].fr))==(r2=Get_fa(e[i].to))) continue;
        anc[r1]=anc[r2]=fa[r1][0]=fa[r2][0]=++tot, anc[tot]=fa[tot][0]=tot/*!*/;//清空新建的fa[tot]!(可能作為根節點)
        mn[tot]=e[i].h, Ans[tot]=std::min(Ans[r1],Ans[r2]);
        if(++k==n) break;
    }
}
void Init_ST()
{
    for(int i=1; i<=18; ++i)
        for(int x=1; x<=tot; ++x) fa[x][i]=fa[fa[x][i-1]][i-1];
}
inline int Solve(int p,int ht)
{
    for(int i=18; ~i; --i)
        if(mn[fa[p][i]]>ht) p=fa[p][i];
    return Ans[p];
}

int main()
{
//  freopen("return.in","r",stdin);
//  freopen("return.out","w",stdout);

    int Case=read();
    while(Case--)
    {
        Enum=0, memset(H,0,sizeof H);
        tot=n=read(), m=read();
        while(m--) AddEdge(read(),read(),read(),read());
        Dijkstra(), Kruskal(), Init_ST();

        int Q=read(),K=read(),S=read(),ans=0,pos,ht;
        if(K) while(Q--)
            pos=(read()+ans-1)%n+1,ht=(read()+ans)%(S+1),printf("%d\n",ans=Solve(pos,ht));
        else while(Q--)
            pos=read(),ht=read(),printf("%d\n",Solve(pos,ht));
    }
    return 0;
}

洛谷.4768.[NOI2018]歸程(Kruskal重構樹 倍增)