F. 数学上来先打表

题解：

搞这题搞了一天

思维不是很难

就是暴力压位bitset

分块做法速度更快

但是stl里的不能实现这个功能

所以手动实现

64位压一位

到65535跑一下1的个数

然后(x>>16)&65535+...计算出1的个数

眼查错果然还是不靠谱以后还是写对拍吧

另外的难点在于这题很卡空间

我们建立dfs树

这样把总空间变为最大消耗空间

合并bitset的时候我们利用启发式合并

并且bitset只存出现了的位置

就是那些64位都没有出现的就不存了

空间是nlogn的，但是还是挺卡的吧

使用动态开数组的技巧替代vector

*int=calloc(大小，sizeof（））

另外这个东西不会自动释放空间

所以我们要free（）这个东西

代码：

#include <bits/stdc++.h>
using namespace std;
#define ull unsigned long long
#define ll long long
#define rint int
#define IL inline
#define rep(i,h,t) for (rint i=h;i<=t;i++)
#define dep(i,t,h) for (rint i=t;i>=h;i--)
const int N=1.1e5;
struct bi{
  ull *v;
  int *pos,cnt;
  bi() { cnt=;}
}B[N];
int head[N],count2[N],count3[N],pos[N],pos2[N],pos3[N],pos4[N],fa[N],ans[N];
IL void js2(bi &a,bi b,int p)
{
  if (p==)
  {
    int sad;
    sad=;
  }
    bi c;
    int tmp=a.cnt+b.cnt;
    c.v=(ull*)calloc(tmp+,sizeof(ull));
    c.pos=(int*)calloc(tmp+,sizeof(int));
    rint x=a.cnt,y=b.cnt;
    c.cnt=;
    rint i=,j=;
    while (i<=x&&j<=y)
    {
      if (a.pos[i]<b.pos[j]) c.v[++c.cnt]=a.v[i],c.pos[c.cnt]=a.pos[i],i++;
      else if (a.pos[i]>b.pos[j]) c.v[++c.cnt]=b.v[j],c.pos[c.cnt]=b.pos[j],j++;
      else
      {
  //      if (a.v[i]^b.v[j])
        c.v[++c.cnt]=a.v[i]^b.v[j],c.pos[c.cnt]=b.pos[j];
        i++; j++;
      }
    }
    while (j<=y) c.v[++c.cnt]=b.v[j],c.pos[c.cnt]=b.pos[j],j++;
    while (i<=x) c.v[++c.cnt]=a.v[i],c.pos[c.cnt]=a.pos[i],i++;
    free(a.v); free(a.pos);
    a=c;
}
struct re{
  int a,b,c,d;
}e[N*],a[N];
bool cmp(re x,re y)
{
  return(x.a<y.a);
}
int l;
void arr(int x,int y,int x1,int y1)
{
  e[++l].a=head[x];
  e[l].b=y;
  e[l].c=x1;
  e[l].d=y1;
  head[x]=l;
}
#define js(x) (count2[x&65535]+count2[(x>>16)&65535]+count2[(x>>32)&65535]+count2[(x>>48)&65535])
int find(int x)
{
  return fa[x]==x?x:find(fa[x]);
}
vector<re> ve[N];
void dfs(int x)
{
  int l=(int)(ve[x].size())-;
    rep(i,,l)
    {
      int x1=pos4[find(ve[x][i].a)],x2=ve[x][i].b,x3=ve[x][i].c;
      int cnt=;
      rep(i,,B[x1].cnt)
        if (cnt+js(B[x1].v[i])<x2)
          cnt+=js(B[x1].v[i]);
        else
        {
          ull o=B[x1].v[i];
          rep(j,,)
            if ((o>>j)&)
            {
              cnt++;
              if (cnt==x2)
              {
                ans[x3]=*(B[x1].pos[i]-)+j+; break;
              }
            }
          break;
        }
      if (cnt<x2) ans[x3]=-;
    }
  for (rint u=head[x];u;u=e[u].a)
  {
    rint v=e[u].b;
    rint x1=find(e[u].c),x2=find(e[u].d),tmp=;
  //  bi tmp2;
    bool tt=;
    if (x1!=x2)
    {
      tt=;
      if (count3[x1]>count3[x2]) swap(x1,x2);
      fa[x1]=x2; count3[x2]+=count3[x1];
      int k1=B[pos4[x1]].cnt,k2=B[pos4[x2]].cnt;
      if (k1>k2)
      {
        tmp=pos4[x2];
      //  tmp2=B[pos4[x1]];
        js2(B[pos4[x1]],B[pos4[x2]],);
        pos4[x2]=pos4[x1];
      } else
      js2(B[pos4[x2]],B[pos4[x1]],pos4[x2]);
    }
    dfs(v);
    if (tt)
    {
      fa[x1]=x1; count3[x2]-=count3[x1];
      if (tmp)
      {
        pos4[x2]=tmp; js2(B[pos4[x1]],B[tmp],);
      } else
      {
        js2(B[pos4[x2]],B[pos4[x1]],pos4[x2]);
      }
    }
  }
}
int main()
{
  ios::sync_with_stdio(false);
  rep(i,,) count2[i]=count2[i>>]+(i&);
  int n,m;
  cin>>n>>m;
  rep(i,,n) cin>>a[i].a,a[i].b=i;
  rep(i,,n) pos4[i]=i,count3[i]=;
  sort(a+,a+n+,cmp);
  rep(i,,n) pos[a[i].b]=i,pos2[i]=a[i].a;
  rep(i,,n)
  {
    B[i].v=(ull*)calloc(,sizeof(ull)); B[i].pos=(int*)calloc(,sizeof(int));
    int t=pos[i]; t--;
    B[i].pos[]=(t/)+; B[i].v[]=1ull<<(t%);
    fa[i]=i;
    B[i].cnt=;
  }
  int now,cnt,cnt2=;
  now=cnt=;
  pos3[]=;
  rep(i,,m)
  {
    int kk,x,y;
    cin>>kk;
    if (kk==)
    {
      cin>>x>>y;
      arr(now,++cnt,x,y);
      pos3[i]=cnt;
      now=cnt;
    }
    if (kk==)
    {
      cin>>x;
      now=pos3[x];
    }
    if (kk==)
    {
      cin>>x>>y;
      ve[now].push_back((re){x,y,++cnt2});
    }
    pos3[i]=now;
  }
  dfs();
  rep(i,,cnt2)
  if (ans[i]==-) cout<<-<<endl;
  else cout<<pos2[ans[i]]<<endl;
  return ;
}