About set HDU

https://vjudge.net/problem/HDU-4680

一直想写，终于写完了。。。

要点：

这个set不需要去重

操作4的做法就是暴力枚举取的数（最开始两个取set中最小两个，设这次取的值为now，前一次取的值为last，那么下次要取大于等于now+last的最小数），由于每一次取的值大于等于last的值的两倍，因此枚举次数是logn级别的

错误记录：

1.100行少upd

2.多组数据忘消去内存占用

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
int n,m;
int gcd(int a,int b)
{
    if(!a||!b)    return a+b;
    int t;
    while(b)
    {
        t=a;
        a=b;
        b=t%b;
    }
    return a;
}
namespace S
{
#define N 400100
struct Node
{
    Node *ch[2];
    int sz,d,r,n;
    int g;
}nds[N];
queue<Node*> q;
int rand1()
{
    static int x=233;
    return x=(48271LL*x+1)%2147483647;
}
void init()
{
    for(int i=0;i<N;i++)    q.push(nds+i);
}
Node *getnode()
{
    Node *t=q.front();q.pop();
    t->ch[0]=t->ch[1]=0;t->sz=t->d=t->g=t->n=0;t->r=rand1();
    return t;
}
void delnode(Node *x)    {q.push(x);}
int gsz(Node *o)    {return o?o->sz:0;}
int gg(Node *o)    {return o?o->g:0;}
void upd(Node *o)
{
    o->g=gcd(gcd(gg(o->ch[0]),o->d),gg(o->ch[1]));
    o->sz=gsz(o->ch[0])+o->n+gsz(o->ch[1]);
}
Node *merge(Node *a,Node *b)
{
    if(!a)    return b;
    if(!b)    return a;
    if(a->r<b->r)
    {
        a->ch[1]=merge(a->ch[1],b);upd(a);
        return a;
    }
    else
    {
        b->ch[0]=merge(a,b->ch[0]);upd(b);
        return b;
    }
}
typedef pair<Node*,Node*> pnn;
pnn split_x(Node *a,int x)//<=x,>x
{
    if(!a)    return pnn(0,0);
    pnn y;
    if(x<a->d)
    {
        y=split_x(a->ch[0],x);
        a->ch[0]=y.se;upd(a);y.se=a;
    }
    else
    {
        y=split_x(a->ch[1],x);
        a->ch[1]=y.fi;upd(a);y.fi=a;
    }
    return y;
}
void add(Node *&a,int x,int k)
{
    pnn t1=split_x(a,x-1),t2=split_x(t1.se,x);
    if(!t2.fi)
    {
        Node *t=getnode();t->d=x;upd(t);
        t2.fi=t;
    }
    t2.fi->n+=k;upd(t2.fi);//upd
    if(t2.fi->n==0)    delnode(t2.fi),t2.fi=0;
    a=merge(t1.fi,merge(t2.fi,t2.se));
}
int qg(Node *&a,int l,int r)
{
    pnn t1=split_x(a,l-1),t2=split_x(t1.se,r);
    int ans=t2.fi?t2.fi->g:-1;
    a=merge(t1.fi,merge(t2.fi,t2.se));
    return ans;
}
Node *find(Node *a,int x)
{
    if(!a)    return 0;
    if(x<a->d)    return find(a->ch[0],x);
    else if(x==a->d)    return a;
    else    return find(a->ch[1],x);
}
int upb(Node *a,int x)
{
    if(!a)    return 0x7fffffff;
    if(a->d>x)    return min(a->d,upb(a->ch[0],x));
    else return upb(a->ch[1],x);
}
void mergeto(Node *&a,Node *b)//b合并到a上
{
    if(!b)    return;
    mergeto(a,b->ch[0]);mergeto(a,b->ch[1]);
    add(a,b->d,b->n);
    delnode(b);
}
void mgto(Node *&a,Node *&b)//b合并到a上，启发式
{
    if(gsz(a)<gsz(b))    swap(a,b);
    mergeto(a,b);
}
void deltree(Node *o)
{
    if(!o)    return;
    deltree(o->ch[0]);deltree(o->ch[1]);
    delnode(o);
}
}
using S::Node;
using S::add;
using S::upb;
Node *rt[N];
int a[N];
namespace D
{
int fa[N],p[N],num;
void init()
{
    int i;num=n;
    for(i=1;i<=n+m;i++)    fa[i]=i,rt[i]=0;
    for(i=1;i<=n;i++)
    {
        p[i]=i;
        add(rt[i],a[i],1);
    }
}
int findd(int x){return x==fa[x]?x:fa[x]=findd(fa[x]);}
int find(int x){return findd(p[x]);}
void unionn(int x,int y)
{
    x=findd(x);y=findd(y);
    if(x==y)    return;
    fa[x]=y;S::mgto(rt[y],rt[x]);
}
void union1(int x,int y){unionn(p[x],p[y]);}
void move(int x,int y)//将x移到y所在集合中
{
    if(x==y)    return;
    add(rt[find(x)],a[x],-1);
    ++num;p[x]=num;unionn(p[x],p[y]);
    add(rt[find(x)],a[x],1);
}
void change(int x,int y)//x权值改y
{
    add(rt[find(x)],a[x],-1);
    a[x]=y;
    add(rt[find(x)],a[x],1);
}
}
int main()
{
    Node *tt,*r;
    int T,TT=0,i,ans,u,v,idx,last,now,t,L,R;
    S::init();
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        for(i=1;i<=n;i++)    scanf("%d",&a[i]);
        ++TT;
        printf("Case #%d:
",TT);
        D::init();
        for(i=1;i<=m;i++)
        {
            scanf("%d",&idx);
            if(idx==1)
            {
                scanf("%d%d",&u,&v);
                D::union1(u,v);
            }
            else if(idx==2)
            {
                scanf("%d%d",&u,&v);
                D::move(u,v);
            }
            else if(idx==3)
            {
                scanf("%d%d",&u,&v);
                D::change(u,v);
            }
            else if(idx==4)
            {
                scanf("%d",&u);r=rt[D::find(u)];
                if(r->sz<=2)    ans=r->sz;
                else
                {
                    last=upb(r,-1);tt=S::find(r,last);
                    if(tt->n>1)    now=last;
                    else    now=upb(r,last);
                    ans=2;
                    while(1)
                    {
                        t=last;last=now;
                        now=upb(r,t+now-1);
                        if(now==0x7fffffff)    break;
                        ans++;
                    }
                }
                printf("%d
",ans);
            }
            else if(idx==5)
            {
                scanf("%d%d%d",&u,&L,&R);
                printf("%d
",S::qg(rt[D::find(u)],L,R));
            }
        }
        for(i=1;i<=n+m;i++)
            if(D::findd(i)==i)
                S::deltree(rt[i]);
        //printf("%d
",S::q.size());
    }
    return 0;
}