題意：給你n個數，q個操作，操作k，l，r，k=1時 區間[l,r]每個數加x，k=2時，區間[l,r]每個數開平方，k=3時，求區間[l,r]的和。

分析：我們知道一個數多次開平方會變成1，但是這裡的1操作會使這個數的值增大，所以直接判斷一個區間是否為1肯定超時。

官方題解加了個優化，就是判斷一個區間內是否為一個值，因為有些區間在操作1和2下會變成一個值。

但是有一組資料把標程hack了：10萬個2，3，2，3，2，3.......,10萬個操作 加6，開根。

這樣每個數相等的區間大小就為1了，標程T了。

我們來考慮這種資料，如果一個區間的極差>1,那麼經過多次1，2操作後，最終這個區間的極差要麼是0，要麼是1，那麼我們只要維護極差<=1的區間就好了。

後來也T了很久。。。要用讀入優化和輸入加到build裡，G++提交。

程式碼：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<vector>
#include<queue>
#include<cmath>
#include<stack>
#include<set>
#include<map>
#define INF 0x3f3f3f3f
#define Mn 100010
#define Mm 2000005
#define mod 1000000007
#define CLR(a,b) memset((a),(b),sizeof((a)))
#define CLRS(a,b,Size) memset((a),(b),sizeof((a[0]))*(Size+1))
#define CPY(a,b) memcpy ((a), (b), sizeof((a)))
#pragma comment(linker, "/STACK:102400000,102400000")
#define ul u<<1
#define ur (u<<1)|1
using namespace std;
typedef long long ll;
int read() {
    char c=getchar();
    int re=0,f=1;
    while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
    while(c>='0'&&c<='9') {re=re*10+c-'0';c=getchar();}
    return re*f;
}
int a[Mn];
ll sum[Mn*4];
ll maxx[Mn*4],minn[Mn*4];
ll maxNum[Mn*4],minNum[Mn*4];
ll all[Mn*4];
ll lazy[Mn*4];
void pushUp(int u) {
    sum[u]=sum[ul]+sum[ur];
    maxx[u]=max(maxx[ul],maxx[ur]);
    minn[u]=min(minn[ul],minn[ur]);
    maxNum[u]=minNum[u]=0;
    if(maxx[u]==maxx[ul]) maxNum[u]+=maxNum[ul];
    if(maxx[u]==maxx[ur]) maxNum[u]+=maxNum[ur];
    if(minn[u]==minn[ul]) minNum[u]+=minNum[ul];
    if(minn[u]==minn[ur]) minNum[u]+=minNum[ur];
}
void pushDown(int u,int l,int r) {
    if(all[u]) {
        int mid=(l+r)>>1;
        all[ul]=all[ur]=all[u];
        sum[ul]=all[u]*(mid-l+1);
        sum[ur]=all[u]*(r-mid);
        maxx[ul]=maxx[ur]=minn[ul]=minn[ur]=all[u];
        maxNum[ul]=minNum[ul]=mid-l+1;
        maxNum[ur]=minNum[ur]=r-mid;
        lazy[ul]=lazy[ur]=0;
        all[u]=0;
    }
    if(lazy[u]) {
        int mid=(l+r)>>1;
        lazy[ul]+=lazy[u];lazy[ur]+=lazy[u];
        sum[ul]+=lazy[u]*(mid-l+1);
        sum[ur]+=lazy[u]*(r-mid);
        maxx[ul]+=lazy[u];minn[ul]+=lazy[u];
        maxx[ur]+=lazy[u];minn[ur]+=lazy[u];
        lazy[u]=0;
    }
}
void build(int l,int r,int u) {
    all[u]=lazy[u]=maxx[u]=minn[u]=sum[u]=0;
    if(l==r) {
        maxx[u]=minn[u]=sum[u]=read();
        maxNum[u]=minNum[u]=1;
        return ;
    }
    int mid=(l+r)>>1;
    build(l,mid,ul);
    build(mid+1,r,ur);
    pushUp(u);
}
int s,t;
void add(int l,int r,int u,int x) {
    if(s<=l&&t>=r) {
        sum[u]+=(ll)x*(r-l+1);
        minn[u]+=x;
        maxx[u]+=x;
        lazy[u]+=x;
        return ;
    }
    pushDown(u,l,r);
    int mid=(l+r)>>1;
    if(s<=mid) add(l,mid,ul,x);
    if(t>mid) add(mid+1,r,ur,x);
    pushUp(u);
}
void Sqrt(int l,int r,int u) {
    if(s<=l&&t>=r&&(maxx[u]-minn[u]<=1)) {
        if(maxx[u]==1) return ;
        int x=maxx[u];
        if(maxx[u]==minn[u]) {
            minn[u]=maxx[u]=floor(sqrt(x));
            lazy[u]+=maxx[u]-x;
            sum[u]=maxx[u]*(r-l+1);
            return ;
        } else if(maxx[u]-minn[u]==1){
            maxx[u]=floor(sqrt(maxx[u]));
            minn[u]=floor(sqrt(minn[u]));
            if(maxx[u]==minn[u]) {
                maxNum[u]=minNum[u]=r-l+1;
                lazy[u]=0;
                all[u]=maxx[u];
                sum[u]=maxx[u]*maxNum[u];
            } else if(maxx[u]-minn[u]==1) {
                lazy[u]+=maxx[u]-x;
                sum[u]=maxx[u]*maxNum[u]+minn[u]*minNum[u];
            }
            return ;
        }
        return ;
    }
    pushDown(u,l,r);
    int mid=(l+r)>>1;
    if(s<=mid) Sqrt(l,mid,ul);
    if(t>mid) Sqrt(mid+1,r,ur);
    pushUp(u);
}
ll query(int l,int r,int u) {
    if(s<=l&&t>=r) {
        return sum[u];
    }
    pushDown(u,l,r);
    int mid=(l+r)>>1;
    ll re=0;
    if(s<=mid) re+=query(l,mid,ul);
    if(t>mid) re+=query(mid+1,r,ur);
    return re;
}
int main() {
    int T=read(),k;
    while(T--) {
        int n=read(),m=read();
        build(1,n,1);
        while(m--) {
            k=read(),s=read(),t=read();
            if(k==1) {
                int x=read();
                add(1,n,1,x);
            } else if(k==2){
                Sqrt(1,n,1);
            } else {
                printf("%I64d\n",query(1,n,1));
            }
        }
    }
    return 0;
}