题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6280
#include<bits/stdc++.h>
#define N 5005
using namespace std;
typedef long long ll;
const ll mod = (ll)1e9+7;
const int maxn=5005;
struct node{
int to,next;
}e[maxn<<1];
int n,m,a,b,x0,y0,f[maxn],ans,sz[maxn];
int head[maxn],now,dep[maxn],fa[maxn][25];
void dfs(int no,int ff,int deps){
dep[no]=deps;
fa[no][0]=ff;
for(int i=head[no];~i;i=e[i].next){
int u=e[i].to;
if(u==ff) continue;
if(!dep[u]){
dfs(u,no,deps+1);
}
}
}
void add(int u,int v){
e[now].to=v;
e[now].next=head[u],head[u]=now++;
}
void deal(){
for(int j=1;(1<<j)<=n;j++)
for(int i=1;i<=n;i++)
fa[i][j]=fa[fa[i][j-1]][j-1];
}
int LCA(int x,int y){
if(dep[x]<dep[y]) swap(x,y);
for(int i=24;i>=0;i--)
if(dep[x]-(1<<i)>=dep[y])
x=fa[x][i];
if(x==y) return x;
for(int i=24;i>=0;i--)
if(fa[x][i]!=fa[y][i]&&fa[x][i]!=-1)
x=fa[x][i],y=fa[y][i];
return fa[x][0];
}
int fin(int x){
return f[x]==x?x:f[x]=fin(f[x]);
}
void Merge(int x,int y){
x=fin(x);
if(dep[fa[x][0]]<=dep[y]||fa[x][0]==-1){
return ;
}
ans=ans^sz[fa[x][0]]^(sz[fa[x][0]]-1);
sz[fa[x][0]]--;
f[x]=fa[x][0];
Merge(fa[x][0],y);
}
int main(){
while(~scanf("%d%d%d%d%d%d",&n,&m,&a,&b,&x0,&y0)){
now=0;
for(int i=0;i<=n;i++) f[i]=i,head[i]=-1,sz[i]=0,dep[i]=0;
memset(fa,-1,sizeof(-1));
for(int i=1;i<n;i++){
int uu,vv;
scanf("%d%d",&uu,&vv);
uu++,vv++;
sz[uu]++,sz[vv]++;
add(uu,vv); add(vv,uu);
}
dfs(1,-1,1);
deal();
ans=0;
for(int i=1;i<=n;i++){
ans^=sz[i];
}
for(int i=1;i<=m;i++){
int nx=(a*x0+b*y0+ans)%n;
int ny=(a*y0+b*x0+ans)%n;
x0=nx,y0=ny;
Merge(x0+1,LCA(x0+1,y0+1));
}
printf("%d %d\n",x0,y0);
}
return 0;
}