loj #6281. 数列分块入门 5

#6281. 数列分块入门 5

题目描述

给出一个长为 nnn 的数列 a1…ana_1\ldots a_na1​​an​​,以及 nnn 个操作,操作涉及区间开方,区间求和。

输入格式

第一行输入一个数字 nnn。

第二行输入 nnn 个数字,第 i 个数字为 aia_iai​​,以空格隔开。

接下来输入 nnn 行询问,每行输入四个数字 opt,l,r,c\mathrm{opt}, l, r, copt,l,r,c,以空格隔开。

若 opt=0\mathrm{opt} = 0opt=0,表示将位于 [l,r][l, r][l,r] 的之间的数字都开方。对于区间中每个 ai(lir),aiai

若 opt=1\mathrm{opt} = 1opt=1,表示询问位于 [l,r][l, r][l,r] 的所有数字的和。

输出格式

对于每次询问,输出一行一个数字表示答案。

样例
样例输入
4
1 2 2 3
0 1 3 1
1 1 4 4
0 1 2 2
1 1 2 4
样例输出
6
2
数据范围与提示

对于 100% 100\%100% 的数据,1≤n≤50000,−231≤others 1 \leq n \leq 50000, -2^{31} \leq \mathrm{others}1n50000,231​​others、ans≤231−1 \mathrm{ans} \leq 2^{31}-1ans231​​1。

 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define maxn 100010
using namespace std;
long long n,m,a[maxn],sum[maxn],block,pos[maxn];
bool v[maxn];
int findl(int id){return (id-1)*block+1;}
int findr(int id){return min(id*block,n);}
long long qread(){
    long long i=0,j=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')j=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){i=i*10+ch-'0';ch=getchar();}
    return i*j;
}
void reset(int id){
    if(v[id])return;
    v[id]=1;sum[id]=0;
    for(int i=findl(id);i<=findr(id);i++){
        a[i]=sqrt(a[i]);
        sum[id]+=a[i];
        if(a[i]>1)v[id]=0;
    }
}
void update(int l,int r){
    if(!v[pos[l]]){
        for(int i=l;i<=min(findr(pos[l]),r);i++){
            sum[pos[l]]-=a[i];
            a[i]=sqrt(a[i]);
            sum[pos[l]]+=a[i];
        }
        v[pos[l]]=1;
        for(int i=findl(pos[l]);i<=findr(pos[l]);i++)
            if(a[i]>1)v[pos[l]]=0;
    }
    if(pos[l]==pos[r])return;
    if(!v[pos[r]]){
        for(int i=findl(pos[r]);i<=r;i++){
            sum[pos[r]]-=a[i];
            a[i]=sqrt(a[i]);
            sum[pos[r]]+=a[i];
        }
        v[pos[r]]=1;
        for(int i=findl(pos[r]);i<=findr(pos[r]);i++)
            if(a[i]>0)v[pos[r]]=0;
    }
    for(int i=pos[l]+1;i<=pos[r]-1;i++)reset(i);
}
void query(int l,int r){
    long long res=0;
    for(int i=l;i<=min(findr(pos[l]),r);i++)res+=a[i];
    if(pos[l]!=pos[r])
        for(int i=findl(pos[r]);i<=r;i++)res+=a[i];
    for(int i=pos[l]+1;i<=pos[r]-1;i++)res+=sum[i];
    cout<<res<<endl;
}
int main(){
    n=qread();
    for(int i=1;i<=n;i++)a[i]=qread();
    block=sqrt(n);
    for(int i=1;i<=n;i++){
        pos[i]=(i-1)/block+1;
        sum[pos[i]]+=a[i];
    }
    for(int i=1;i<=n;i++){
        long long op=qread(),l=qread(),r=qread(),c=qread();
        if(op==0)update(l,r);
        else query(l,r);
    }
    return 0;
}

 

转载于:https://www.cnblogs.com/thmyl/p/8965930.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值