About set HDU - 4680

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[];
     int sz,d,r,n;
     int g;
 }nds[N];
 queue<Node*> q;
 int rand1()
 {
     static int x=;
     return x=(48271LL*x+)%;
 }
 void init()
 {
     for(int i=;i<N;i++)    q.push(nds+i);
 }
 Node *getnode()
 {
     Node *t=q.front();q.pop();
     t->ch[]=t->ch[]=;t->sz=t->d=t->g=t->n=;t->r=rand1();
     return t;
 }
 void delnode(Node *x)    {q.push(x);}
 int gsz(Node *o)    {return o?o->sz:;}
 int gg(Node *o)    {return o?o->g:;}
 void upd(Node *o)
 {
     o->g=gcd(gcd(gg(o->ch[]),o->d),gg(o->ch[]));
     o->sz=gsz(o->ch[])+o->n+gsz(o->ch[]);
 }
 Node *merge(Node *a,Node *b)
 {
     if(!a)    return b;
     if(!b)    return a;
     if(a->r<b->r)
     {
         a->ch[]=merge(a->ch[],b);upd(a);
         return a;
     }
     else
     {
         b->ch[]=merge(a,b->ch[]);upd(b);
         return b;
     }
 }
 typedef pair<Node*,Node*> pnn;
 pnn split_x(Node *a,int x)//<=x,>x
 {
     if(!a)    return pnn(,);
     pnn y;
     if(x<a->d)
     {
         y=split_x(a->ch[],x);
         a->ch[]=y.se;upd(a);y.se=a;
     }
     else
     {
         y=split_x(a->ch[],x);
         a->ch[]=y.fi;upd(a);y.fi=a;
     }
     return y;
 }
 void add(Node *&a,int x,int k)
 {
     pnn t1=split_x(a,x-),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==)    delnode(t2.fi),t2.fi=;
     a=merge(t1.fi,merge(t2.fi,t2.se));
 }
 int qg(Node *&a,int l,int r)
 {
     pnn t1=split_x(a,l-),t2=split_x(t1.se,r);
     int ans=t2.fi?t2.fi->g:-;
     a=merge(t1.fi,merge(t2.fi,t2.se));
     return ans;
 }
 Node *find(Node *a,int x)
 {
     if(!a)    return ;
     if(x<a->d)    return find(a->ch[],x);
     else if(x==a->d)    return a;
     else    return find(a->ch[],x);
 }
 int upb(Node *a,int x)
 {
     if(!a)    return 0x7fffffff;
     if(a->d>x)    return min(a->d,upb(a->ch[],x));
     else return upb(a->ch[],x);
 }
 void mergeto(Node *&a,Node *b)//b合并到a上
 {
     if(!b)    return;
     mergeto(a,b->ch[]);mergeto(a,b->ch[]);
     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[]);deltree(o->ch[]);
     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=;i<=n+m;i++)    fa[i]=i,rt[i]=;
     for(i=;i<=n;i++)
     {
         p[i]=i;
         add(rt[i],a[i],);
     }
 }
 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],-);
     ++num;p[x]=num;unionn(p[x],p[y]);
     add(rt[find(x)],a[x],);
 }
 void change(int x,int y)//x权值改y
 {
     add(rt[find(x)],a[x],-);
     a[x]=y;
     add(rt[find(x)],a[x],);
 }
 }
 int main()
 {
     Node *tt,*r;
     int T,TT=,i,ans,u,v,idx,last,now,t,L,R;
     S::init();
     scanf("%d",&T);
     while(T--)
     {
         scanf("%d%d",&n,&m);
         for(i=;i<=n;i++)    scanf("%d",&a[i]);
         ++TT;
         printf("Case #%d:\n",TT);
         D::init();
         for(i=;i<=m;i++)
         {
             scanf("%d",&idx);
             if(idx==)
             {
                 scanf("%d%d",&u,&v);
                 D::union1(u,v);
             }
             else if(idx==)
             {
                 scanf("%d%d",&u,&v);
                 D::move(u,v);
             }
             else if(idx==)
             {
                 scanf("%d%d",&u,&v);
                 D::change(u,v);
             }
             else if(idx==)
             {
                 scanf("%d",&u);r=rt[D::find(u)];
                 if(r->sz<=)    ans=r->sz;
                 else
                 {
                     last=upb(r,-);tt=S::find(r,last);
                     if(tt->n>)    now=last;
                     else    now=upb(r,last);
                     ans=;
                     while()
                     {
                         t=last;last=now;
                         now=upb(r,t+now-);
                         if(now==0x7fffffff)    break;
                         ans++;
                     }
                 }
                 printf("%d\n",ans);
             }
             else if(idx==)
             {
                 scanf("%d%d%d",&u,&L,&R);
                 printf("%d\n",S::qg(rt[D::find(u)],L,R));
             }
         }
         for(i=;i<=n+m;i++)
             if(D::findd(i)==i)
                 S::deltree(rt[i]);
         //printf("%d\n",S::q.size());
     }
     return ;
 }