[BZOJ3211]花神游历各国

题目链接:

传送门

题目分析:

题意简化:给你一个元素非负的序列,要求支持区间开方,区间求和。

做法很多,洛谷的题解里线段树/树状数组/分块/...都有(基本就是数据结构的群魔乱舞),不过分块能跑过的话分块就好了。
记录区间\(i\)的元素总和为\(sum_i\),如果\(sum_i=R[i]-L[i]+1\)的话说明这个区间内只有1了,就不用再操作了,否则直接暴力开方即可。
由于数据范围很小,最多开个10次方上下就全都变成1了,所以时间上是没有问题的。
再有一个细节要处理是数据可能包含0(这一点洛谷保证没有0的数据范围顿时变得一万个良心),所以在读入的时候如果读到0就把这个元素加成1,并另外维护一个数组记录此元素的位置,再维护一个数组记录每个块内0的个数,最后求和的时候减去即可。
数据范围很大,不开\(long long\)跑不过。

代码:

#include<bits/stdc++.h>
#define N (200000+5)
using namespace std;
inline long long read(){
    long long cnt=0,f=1;char c;
    c=getchar();
    while(!isdigit(c)){
        if(c=='-') f=-f;
        c=getchar();
    }
    while(isdigit(c)){
        cnt=cnt*10+c-'0';
        c=getchar();
    }
    return cnt*f;
}
long long n,m,k,l,r;
long long a[N],sum[N];
long long L[N],R[N],pos[N];
long long pos_of_zero[N],sum_of_zero[N];
void Sqrt(int l,int r) {
    int p=pos[l];int q=pos[r];
    if(p==q) {
        if(sum[p]==R[p]-L[p]+1) return;
        else {
            for(register int i=l;i<=r;i++) {
                sum[p]-=a[i];
                a[i]=sqrt(a[i]);
                sum[p]+=a[i];
            }
        }
    }
    else {
        if(sum[p]>R[p]-L[p]+1) {
            for(register int i=l;i<=R[p];i++) {
                sum[p]-=a[i];
                a[i]=sqrt(a[i]);
                sum[p]+=a[i];
            }
        }
        if(sum[q]>R[q]-L[q]+1) {
            for(register int i=L[q];i<=r;i++) {
                sum[q]-=a[i];
                a[i]=sqrt(a[i]);
                sum[q]+=a[i];
            }
        }
        for(register int i=p+1;i<=q-1;i++) {
            if(sum[i]==R[i]-L[i]+1) continue;
            else {
                for(register int j=L[i];j<=R[i];j++) {
                    sum[i]-=a[j];
                    a[j]=sqrt(a[j]);
                    sum[i]+=a[j];
                }
            }
        }
    }
}
long long query(int l,int r) {
    int p=pos[l];int q=pos[r];
    long long ans = 0;
    if(p==q) {
        for(register int i=l;i<=r;i++) {
            ans+=a[i];ans-=pos_of_zero[i];
        }
    }
    else {
        for(register int i=l;i<=R[p];i++) {
            ans+=a[i];ans-=pos_of_zero[i];
        }
        for(register int i=L[q];i<=r;i++) {
            ans+=a[i];ans-=pos_of_zero[i];
        }
        for(register int i=p+1;i<=q-1;i++) {
            ans+=sum[i];ans-=sum_of_zero[i];
        }
    }
    return ans;
}
int main(){
    n=read();
    for(register int i=1;i<=n;i++) {
        a[i]=read();
        if(a[i]==0) {
            pos_of_zero[i]=1;
            a[i]++;
        }
    }
    int t=sqrt(n);
    for(register int i=1;i<=t;i++){
        L[i]=(i-1)*t+1;
        R[i]=i*t;
    }
    if(R[t]<n) {
        ++t;
        L[t]=R[t-1]+1;
        R[t]=n;
    }
    for(register int i=1;i<=t;i++)
        for(register int j=L[i];j<=R[i];j++){
            pos[j]=i;
            sum[i]+=a[j];
            if(pos_of_zero[j]) sum_of_zero[i]++;
        }
    m=read();
    for(register int i=1;i<=m;i++) {
        k=read();l=read();r=read();
        if(l>r) swap(l,r);
        if(k==2) Sqrt(l,r);
        else printf("%lld\n",query(l,r));
    }
    return 0;
}

转载于:https://www.cnblogs.com/kma093/p/10336340.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值