「Poetize4」上帝造题的七分钟2

描述 Description

"第一分钟，X说，要有数列，于是便给定了一个正整数数列。
第二分钟，L说，要能修改，于是便有了对一段数中每个数都开平方(下取整)的操作。
第三分钟，k说，要能查询，于是便有了求一段数的和的操作。
第四分钟，彩虹喵说，要是noip难度，于是便有了数据范围。
第五分钟，诗人说，要有韵律，于是便有了时间限制和内存限制。
第六分钟，和雪说，要省点事，于是便有了保证运算过程中及最终结果均不超过64位有符号整数类型的表示范围的限制。
第七分钟，这道题终于造完了，然而，造题的神牛们再也不想写这道题的程序了。"
——《上帝造题的七分钟·第二部》
所以这个神圣的任务就交给你了。

输入格式 InputFormat

第一行一个整数n，代表数列中数的个数。
第二行n个正整数，表示初始状态下数列中的数。
第三行一个整数m，表示有m次操作。
接下来m行每行三个整数k,l,r，k=0表示给[l,r]中的每个数开平方(下取整)，k=1表示询问[l,r]中各个数的和。
UPD：注意数据中有可能l>r，所以遇到这种情况请交换l和r。

题解：
以前就做过，因为每一个数至多被开log2log2 x 次就变成1了，所以之前我们可以暴力开方，然后成了1之后就不用管他，直接把跳过它去开放没有变成1的
这一步可以用并查集来搞
代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 110000

 #define maxm 500+100

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 1000000007

 using namespace std;

 inline ll read()

 {

     ll x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }
 int n,m,fa[maxn];
 ll s[maxn],a[maxn];
 inline void add(int x,ll y)
 {
     for(;x<=n;x+=x&(-x))s[x]+=y;
 }
 inline ll sum(int x)
 {
     ll t=;
     for(;x;x-=x&(-x))t+=s[x];
     return t;
 }
 inline int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     n=read();
     for1(i,n)add(i,a[i]=read());
     for1(i,n+)fa[i]=i;
     m=read();
     while(m--)
     {
         int ch=read(),l=read(),r=read();if(l>r)swap(l,r);
         if(ch)printf("%lld\n",sum(r)-sum(l-));
         else
         {
             for(int i=find(l);i<=r;i=find(i+))
             {
                 int t=(int)sqrt(a[i]);
                 add(i,t-a[i]);
                 a[i]=t;
                 if(a[i]<=)fa[i]=find(i+);
             }
         }
     }

     return ;

 }