[LOJ6281]数列分块入门 5

题目大意：
　　给你一个长度为$n(nleq50000)$的序列$A(0leq A_i<2^{31})$，支持进行以下两种操作：
　　　　1.将区间$[l,r]$中所有数开方；
　　　　2.询问区间$[l,r]$的和。
思路：
　　分块。
　　因为当一个区间内的所有$A_i$都$leq1$时，操作1对答案没有影响，因此我们可以记录每个区间是否还有操作的必要。如果还有大于$1$的数，就暴力进行操作，否则跳过。易证每个区间最多开方4次。

#include<cmath>
#include<cstdio>
#include<cctype>
#include<algorithm>
typedef long long int64;
inline int getint() {
    register char ch;
    register bool neg=false;
    while(!isdigit(ch=getchar())) if(ch=='-') neg=true;
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return neg?-x:x;
}
const int N=50001;
int val[N],bel[N],begin[N],end[N],max[N];
int64 sum[N];
inline void modify(const int &l,const int &r,const int &c) {
    if(bel[l]==bel[r]) {
        for(register int i=l;i<=r;i++) val[i]=sqrt(val[i]);
        sum[bel[l]]=max[bel[l]]=0;
        for(register int i=begin[bel[l]];i<=end[bel[l]];i++) {
            sum[bel[l]]+=val[i];
            max[bel[l]]=std::max(max[bel[l]],val[i]);
        }
        return;
    }
    for(register int i=l;bel[i]==bel[l];i++) val[i]=sqrt(val[i]);
    for(register int i=r;bel[i]==bel[r];i--) val[i]=sqrt(val[i]);
    for(register int i=bel[l]+1;i<bel[r];i++) {
        if(max[i]==1) continue;
        for(register int j=begin[i];j<=end[i];j++) val[j]=sqrt(val[j]);
    }
    for(register int i=bel[l];i<=bel[r];i++) {
        if(max[i]==1) continue;
        sum[i]=max[i]=0;
        for(register int j=begin[i];j<=end[i];j++) {
            sum[i]+=val[j];
            max[i]=std::max(max[i],val[j]);
        }
    }
}
inline int64 query(const int &l,const int &r) {
    int64 ret=0;
    if(bel[l]==bel[r]) {
        for(register int i=l;i<=r;i++) ret+=val[i];
        return ret;
    }
    for(register int i=l;bel[i]==bel[l];i++) ret+=val[i];
    for(register int i=r;bel[i]==bel[r];i--) ret+=val[i];
    for(register int i=bel[l]+1;i<bel[r];i++) ret+=sum[i];
    return ret;
}
int main() {
    const int n=getint(),block=sqrt(n);
    for(register int i=1;i<=n;i++) {
        val[i]=getint();
        bel[i]=i/block;
        if(!begin[bel[i]]) begin[bel[i]]=i;
        end[bel[i]]=i;
        sum[bel[i]]+=val[i];
        max[bel[i]]=std::max(max[bel[i]],val[i]);
    }
    for(register int i=0;i<n;i++) {
        const int opt=getint(),l=getint(),r=getint(),c=getint();
        if(opt) {
            printf("%lld
",query(l,r));
        } else {
            modify(l,r,c);
        }
    }
    return 0;
}