另外开一个数组维护每一个块内的总和。

给区间加值是,残余的块一个一个点更新,整个的块一次性更新

查询的时候也是,残余的块一个一个点加,整个的块一次性加

  1. #include<map>
  2. #include<set>
  3. #include<ctime>
  4. #include<cmath>
  5. #include<stack>
  6. #include<queue>
  7. #include<string>
  8. #include<vector>
  9. #include<cstdio>
  10. #include<cstdlib>
  11. #include<cstring>
  12. #include<iostream>
  13. #include<algorithm>
  14. #define first fi
  15. #define second se
  16. #define lowbit(x) (x & (-x))
  17.  
  18. typedef unsigned long long int ull;
  19. typedef long long int ll;
  20. const double pi = 4.0*atan(1.0);
  21. const int inf = 0x3f3f3f3f;
  22. const int maxn = ;
  23. const int maxm = ;
  24. const int mod = ;
  25. using namespace std;
  26.  
  27. int n, m, tol, T;
  28. int block;
  29. ll a[maxn];
  30. ll b[maxn];
  31. ll add[maxn];
  32. int belong[maxn];
  33.  
  34. void init() {
  35. memset(a, , sizeof a);
  36. memset(b, , sizeof b);
  37. memset(add, , sizeof add);
  38. memset(belong, , sizeof belong);
  39. }
  40.  
  41. int L(int x) {
  42. return (x-)*block+;
  43. }
  44.  
  45. int R(int x) {
  46. return min(n, x*block);
  47. }
  48.  
  49. void update(int l, int r, int c) {
  50. for(int i=l; i<=min(r, R(belong[l])); i++) {
  51. a[i] += c;
  52. b[belong[i]] += c;
  53. }
  54. if(belong[l] == belong[r]) return ;
  55. if(belong[l] != belong[r]) {
  56. for(int i=L(belong[r]); i<=r; i++) {
  57. a[i] += c;
  58. b[belong[i]] += c;
  59. }
  60. }
  61. for(int i=belong[l]+; i<belong[r]; i++) {
  62. add[i] += c;
  63. b[i] += c * (R(i) - L(i) + );
  64. }
  65. }
  66.  
  67. ll query(int l, int r, int c) {
  68. ll ans = ;
  69. for(int i=l; i<=min(r, R(belong[l])); i++) {
  70. ans += a[i] + add[belong[i]];
  71. ans %= c;
  72. }
  73. if(belong[l] == belong[r]) return ans;
  74. if(belong[l] != belong[r]) {
  75. for(int i=L(belong[r]); i<=r; i++) {
  76. ans += a[i] + add[belong[i]];
  77. ans %= c;
  78. }
  79. }
  80. for(int i=belong[l]+; i<belong[r]; i++) {
  81. ans += b[i];
  82. ans %= c;
  83. }
  84. return ans % c;
  85. }
  86.  
  87. int main() {
  88. while(~scanf("%d", &n)) {
  89. init();
  90. block = sqrt(n);
  91. for(int i=; i<=n; i++) {
  92. scanf("%lld", &a[i]);
  93. belong[i] = (i-) / block + ;
  94. b[belong[i]] += a[i];
  95. }
  96. m = n;
  97. while(m--) {
  98. int op, l, r, c;
  99. scanf("%d%d%d%d", &op, &l, &r, &c);
  100. if(op == ) {
  101. update(l, r, c);
  102. } else {
  103. ll ans = query(l, r, c+);
  104. printf("%lld\n", ans);
  105. }
  106. }
  107. }
  108. return ;
  109. }

LOJ#6280. 数列分块入门 4的更多相关文章

  1. LOJ #6280. 数列分块入门 4-分块(区间加法、区间求和)

    #6280. 数列分块入门 4 内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论   题目描述 给出一个 ...

  2. LOJ——#6277. 数列分块入门 1

    ~~推荐播客~~ 「分块」数列分块入门1 – 9 by hzwer 浅谈基础根号算法——分块 博主蒟蒻,有缘人可直接观摩以上大佬的博客... #6277. 数列分块入门 1 题目大意: 给出一个长为 ...

  3. LOJ #6285. 数列分块入门 9-分块(查询区间的最小众数)

    #6285. 数列分块入门 9 内存限制:256 MiB时间限制:1500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 2   题目描述 给 ...

  4. LOJ #6284. 数列分块入门 8-分块(区间查询等于一个数c的元素,并将这个区间的所有元素改为c)

    #6284. 数列分块入门 8 内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 2   题目描述 给出 ...

  5. LOJ #6283. 数列分块入门 7-分块(区间乘法、区间加法、单点查询)

    #6283. 数列分块入门 7 内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 2   题目描述 给出 ...

  6. LOJ #6282. 数列分块入门 6-分块(单点插入、单点查询、数据随机生成)

    #6282. 数列分块入门 6 内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 1   题目描述 给出 ...

  7. LOJ #6281. 数列分块入门 5-分块(区间开方、区间求和)

    #6281. 数列分块入门 5 内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 5   题目描述 给出 ...

  8. LOJ #6279. 数列分块入门 3-分块(区间加法、查询区间内小于某个值x的前驱(比其小的最大元素))

    #6279. 数列分块入门 3 内存限制:256 MiB时间限制:1500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 3   题目描述 给 ...

  9. LOJ #6278. 数列分块入门 2-分块(区间加法、查询区间内小于某个值x的元素个数)

    #6278. 数列分块入门 2 内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 6   题目描述 给出 ...

随机推荐

  1. CodeIgniter框架对数据库查询结果进行统计

    假设有一个user表,如果要查询符合条件sex=male的记录数量,有下面几种方法: 方法一:先取回所有符合条件的记录,再count $res = $this->db->query(&qu ...

  2. CMake--模块的使用和自定义模块

    1.链接外部库 如果程序中使用了外部库,事先并不知道它的头文件和链接库的位置,就要给出头文件和链接库的查找方法,并将他们链接到程序中. FIND_PACKAGE(<name> [major ...

  3. 微信小程序错误码参考大全

    开发过程中,会遇到很多微信返回的状态码,鬼知道代表什么意思,现在好了,整理总结了一份状态码,方便大家. 转载:http://www.yiyongtong.com/archives/view-1856- ...

  4. Flutter的Appbar

    actions → List一个 Widget 列表,代表 Toolbar 中所显示的菜单,对于常用的菜单,通常使用 IconButton 来表示:对于不常用的菜单通常使用PopupMenuButto ...

  5. MyBatis基础:MyBatis入门(1)

    1. MyBatis简介 MyBatis 是支持定制化 SQL.存储过程以及高级映射的优秀的持久层框架. MyBatis 避免了几乎所有的 JDBC 代码和手动设置参数以及获取结果集. MyBatis ...

  6. Delphi窗体之间互相调用的简单问题

    问题是这样的,我的程序主窗口Form1上面有一个数据连接(ADOCONNECTION1)和ADOQUERY,然后还有一些数据感知组件用于浏览用的,我打算点击From1中的一个“修改数据”按钮,就弹出F ...

  7. ArrayList 初学小结!

    package good.com; import java.util.ArrayList;//导入 ArrayList 包 调用动态数组! public class ArrayListList { / ...

  8. 想要配置文件生效 需要通过添加到web.xml加载到内存中

    想要配置文件生效 需要通过添加到web.xml加载到内存中

  9. jsp页面给字体加颜色

    jsp页面给字体加颜色<span style="color:red">要加颜色的部分</span>

  10. 51-node-1649齐头并进(最短路)

    题意:中文题,没啥坑点: 解题思路:这道题一开始以为要跑两个最短路,后来发现不用,因为如果给定了铁路的线路,那么,公路一定是n个节点无向图的补图,所以,铁路和公路之间一定有一个是可以直接从1到n的,我 ...