TYVJ 1557 MST+LCA

分析：

可以证明，任意两点间的符合题目要求的路径必然是在MST上的。

有了这个，这题就水了。。。我没有想到，是venique神犇告诉我的

然后在MST上做LCA倍增，顺便维护d[i][j]表示i到他的2^j倍祖先的最大权，就好了~

我是蒟蒻。。LCA写萎了。。看了venique的代码，才发现错。。。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <iostream>

#define N 110000
#define M 1000000

using namespace std;

struct KRU
{
    int a,b,e,d;
}kru[M];

int head[N],next[M],to[M],len[M],cnt,n,m,qu,f[N][22],d[N][22],fa[N],dpt[N],q[N];

inline void add(int u,int v,int w)
{
    to[cnt]=v; len[cnt]=w; next[cnt]=head[u]; head[u]=cnt++;
}

inline bool cmp(const KRU &a,const KRU &b)
{
    return a.d<b.d;
}

void read()
{
    memset(head,-1,sizeof head); cnt=0;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++) scanf("%d%d%d",&kru[i].a,&kru[i].b,&kru[i].d);
}

int findfa(int x)
{
    if(fa[x]!=x) fa[x]=findfa(fa[x]);
    return fa[x];
}

void kruskal()
{
    sort(kru+1,kru+1+m,cmp);
    for(int i=1;i<=n;i++) fa[i]=i;
    int num=1;
    for(int i=1;i<=m;i++)
        if(findfa(kru[i].a)!=findfa(kru[i].b))
        {
            fa[findfa(kru[i].a)]=findfa(kru[i].b);
            num++;
            add(kru[i].a,kru[i].b,kru[i].d); add(kru[i].b,kru[i].a,kru[i].d);
            if(num==n) break;
        }
}

void bfs()
{
    int h=1,t=2,sta;
    q[1]=1; dpt[1]=1;
    while(h<t)
    {
        sta=q[h++];
        for(int i=head[sta];~i;i=next[i])
            if(to[i]!=f[sta][0])
            {
                dpt[to[i]]=dpt[sta]+1;
                f[to[i]][0]=sta;
                d[to[i]][0]=len[i];
                q[t++]=to[i];
            }
    }
}

void initlca()
{
    for(int j=1;j<=20;j++)
        for(int i=1;i<=n;i++)
        {
            d[i][j]=max(d[i][j-1],d[f[i][j-1]][j-1]);
            f[i][j]=f[f[i][j-1]][j-1];
        }
}

int lca(int x,int y)
{
    int rt=0; 
    if(dpt[x]<dpt[y]) swap(x,y);
    for(int i=20;i>=0;i--)
        if(dpt[f[x][i]]>=dpt[y])
        {
            rt=max(rt,d[x][i]);
            x=f[x][i];
        }
    if(x==y) return rt;
    for(int i=20;i>=0;i--)
        if(f[x][i]!=f[y][i])
        {
            rt=max(rt,max(d[x][i],d[y][i]));
            x=f[x][i]; y=f[y][i];
        }
    return max(rt,max(d[x][0],d[y][0]));
}

void go()
{
    kruskal();
    bfs();
    initlca();
    scanf("%d",&qu);
    int si=1,x,y,p1,p2,q1,q2;
    for(int i=1;i<=qu;i++)
    {
        scanf("%d%d%d%d",&p1,&p2,&q1,&q2);
        x=(long long)(si+p1)*p2%n+1;
        y=(long long)(si+q1)*q2%n+1;
        si=lca(x,y);
        printf("%d\n",si);
    }
}

int main()
{
    read();
    go();
    return 0;
}