【题目描述】

给出一个长为 n 的数列，以及 n 个操作，操作涉及区间开方，区间求和。

【输入格式】

第一行输入一个数字 n。

第二行输入 n 个数字，第 i 个数字为 ai，以空格隔开。

接下来输入 n 行询问，每行输入四个数字 opt、l、r、c，以空格隔开。

若 opt=0，表示将位于 [l,r] 的之间的数字都开方，对于区间中的每个 ai(l<=i<=r)，ai->|√ai|

若 opt=1，表示询问位于 [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% 的数据，1<=n<=50000,-2^31<=other,ans<=2^31-1。

【源代码】

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<utility>
#include<stack>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<bitset>
#define EPS 1e-9
#define PI acos(-1.0)
#define INF 0x3f3f3f3f
#define LL long long
const int MOD = 1E9+7;
const int N = 1000000+5;
const int dx[] = {-1,1,0,0,-1,-1,1,1};
const int dy[] = {0,0,-1,1,-1,1,-1,1};
using namespace std;
LL block,sum;//block为块的长度，sum为块的个数
LL a[N];
LL pos[N];//pos记录第i个元素在第几个块中
LL ans[N];//维护整块和
bool flag[N];//标记整块是否全为1
void init(LL n){block=sqrt(n);//块的长度sum=n/block;//块个数if(n%block)sum++;for(LL i=1;i<=n;i++)//第i个元素在第几块中pos[i]=(i-1)/block+1;
}
void work(LL L,LL R,LL x){//维护整块和LL start=0;for(LL i=L;i<=R;i++)start+=a[i];ans[x]=start;
}
void update(LL L,LL R){for(LL i=L;i<=min(pos[L]*block,R);i++){//左边的边角料ans[pos[L]]-=a[i];//开方前元素和减去第i个元素a[i]=sqrt(a[i]);//开方ans[pos[L]]+=a[i];//开方后元素和加上第i个元素}if(pos[L]!=pos[R]){//存在右区间才遍历，防止重复计算for(LL i=(pos[R]-1)*block+1;i<=R;i++){//右边的边角料ans[pos[R]]-=a[i];//开方前元素和减去第i个元素a[i]=sqrt(a[i]);//开方ans[pos[R]]+=a[i];//开方后元素和加上第i个元素}}for(LL i=pos[L]+1;i<=pos[R]-1;i++){//中间的整块if(flag[i])//全为1，跳过continue;else{//不全为1，暴力处理flag[i]=true;//先假设为1，如果不是后面再改为0for(int j=(i-1)*block+1;j<=i*block;j++){//块中的所有元素ans[i]-=a[j];//开方前元素和减去第j个元素a[j]=sqrt(a[j]);//开方ans[i]+=a[j];//开方后元素和加上第j个元素if(a[j]!=1)//不全为1flag[i]=false;}}}
}
LL query(LL L,LL R){LL res=0;for(LL i=L;i<=min(pos[L]*block,R);i++){//左边的边角料res+=a[i];}if(pos[L]!=pos[R]){//存在右区间才遍历，防止重复计算for(LL i=(pos[R]-1)*block+1;i<=R;i++){//右边的边角料res+=a[i];}}for(LL i=pos[L]+1;i<=pos[R]-1;i++){//中间的整块进行二分查找res+=ans[i];}return res;
}
int main(){LL n;scanf("%lld",&n);for(LL i=1;i<=n;i++)scanf("%lld",&a[i]);init(n);for(LL i=1;i<=sum;i++)//传入块的左右边界与编号work((i-1)*block+1,i*block,i);for(LL i=1;i<=n;i++){LL op;LL left,right,x;scanf("%lld",&op);scanf("%lld%lld%lld",&left,&right,&x);if(op==0)update(left,right);else if(op==1)printf("%lld\n",query(left,right));}return 0;
}

数列分块入门 5（LibreOj-6281）相关推荐

「分块」数列分块入门1 – 9
ACM模板放暑假了,回归!!! 自己不会写暴力,而且好久没写代码了,于是学学分块的优雅暴力~ 「分块入门-LibreOJ」「分块」数列分块入门1 – 9 by hzwer 数列简单分块问题实际上有 ...
数列分块入门 (1 ~ 7)
分块 6277. 数列分块入门 1 分块思想我们把每m个元素分成一块,所以我们总共的块数就是n/mn / mn/m块,一般情况下我们取m=nm = \sqrt{n}m=n.对于区间加操作,我们可以 ...
LOJ——#6277. 数列分块入门 1
~~推荐播客~~ 「分块」数列分块入门1 – 9 by hzwer 浅谈基础根号算法--分块博主蒟蒻,有缘人可直接观摩以上大佬的博客... #6277. 数列分块入门 1 题目大意: 给出一个长为 ...
【分块入门】LOJ 数列分块入门 1 - 9 (学习更新……)
dl题解 _「分块」数列分块入门1 – 9 by hzwer LOJ #6277. 数列分块入门 1 题意:给出一个长为n的数列,以及n个操作,操作涉及区间加法,单点查值. 时间限制:100ms 分块 ...
#6279. 数列分块入门 3（区间修改，查询权值前驱）
#6279. 数列分块入门 3 这是使用hzwer建议的set写的分块代码,set自动排序,支持二分查找,但是常数较大,比我下面写的用vector实现的分块慢了三倍,空间大了10倍. #include ...
LOJ 数列分块入门6
LOJ 数列分块入门6 题目: 题目题解: 我都不懂这题为什么要用分块... ... 直接vector就好了... 但是如果有区间修改的话就不行了.所以这题是启示我们也可以动态分块.具体就是每次插入 ...
LOJ #6280. 数列分块入门 4-分块(区间加法、区间求和)
#6280. 数列分块入门 4 内存限制:256 MiB时间限制:500 ms标准输入输出题目类型:传统评测方式:文本比较上传者: hzwer 提交提交记录统计测试数据讨论题目描述给出一个长为 ...
数列分块入门（套题）（loj6277,loj6278,loj6279,loj6280,loj6281,loj6282,loj6283,loj6284,loj6285）
前言 zjoi考差了,码一些分块题缓解一下心情壹数列分块入门 1[loj6277] 题目大意:区间加,单点查直接分块,区间加时完全覆盖的块打tag,边界块暴力重构块大小设为n\sqrt nn ...
「分块系列」数列分块入门3 解题报告
数列分块入门3 题意概括区间加法,区间求前驱. 写在前面这题的方法与分块2方法极其类似,建议自行解决. 正题和上一题类似,但是二分不是用来计数的,而是用来求小于c的最大值的.然后对于不完整快,将 ...
#6277. 数列分块入门 1
题目链接:https://loj.ac/problem/6277 学习博客:http://hzwer.com/8053.html #6277. 数列分块入门 1 内存限制:256 MiB时间限制:10 ...

数列分块入门 5（LibreOj-6281）

【题目描述】

【输入格式】

【输出格式】

【样例】

【数据范围与提示】

【源代码】

数列分块入门 5（LibreOj-6281）相关推荐

最新文章

热门文章