[GZOI2016] 亚索的量子实验【分块】
第二题 亚索的粒子实验
【问题描述】
亚索是一名伟大的科学家,他最近在做一个粒子的实验,粒子初始有一定的能量,实验过程中倘若第i个粒子被注入k能量,那该粒子就会增加k能量,同时由于辐射作用,第2i,3i,4i……(即i的倍数的粒子)也会增加k能量。
他有一台特殊的机器,能够为一段连续的粒子注入一定的能量,此时辐射同样存在,实验过程中,由于他需要取出某些粒子观察一下数据变化,但这十分困扰他,因此,他找到了聪明的你们,希望你们能在他需要数据的时候告诉他相应粒子能量的大小。
【输入格式】
输入文件名为yasuo.in。
输入第 1 行包含 1 个正整数 n ,表示 n 个粒子。
接下来第 2 行有n个数,第i个数ai表示第i个粒子的初始能量值。
第 3 行包括一个数q,表示操作数,包括询问和操作。
接下来 q 行,每行格式如下:
1 i,表示询问第i个粒子当前的能量值。
2 l r d,表示对[l,r]区间中的所有粒子注入d能量。
【输出格式】
输出文件名为yasuo.out。
对于每一个询问输出一行,表示相应粒子的能量值。
【输入输出样例1】
yasuo.in |
yasuo.out |
3 1 2 3 2 2 1 3 5 1 2 |
12 |
【输入输出样例2】
yasuo.in |
yasuo.out |
6 1 2 1 4 5 6 5 2 2 4 2 1 3 1 4 2 3 5 1 1 5 |
3 8 6 |
【数据规模与约定】
对于 30%的数据 0<=n<=100,0<=q<=100.
对于 70%的数据 0<=q<=10^5.
对于 100%的数据 0<=n<=10^5,0<=q<=5*10^5,0<=ai<=10^6
1<=i<=n,1<=l<=r<=n,0<=d<=10^6
半年前做这道题是一脸懵逼,半年后依然懵逼= =,撸了半个多小时线段树+调了1个多小时仍然是只过样例,看了std才发现是分块!
具体实在没啥好说的了,直接看代码好了,代码实在是太清晰了。。。
#include <cstdio>
#include <cmath> const int maxn = 100005; int n, div[maxn][135], idx[maxn], opr, t1, t2, p, lmt, q, a[maxn], x, y;
long long f[maxn], g[maxn], t3, ans; int main(void) {
freopen("yasuo.in", "r", stdin);
freopen("yasuo.out", "w", stdout);
scanf("%d", &n);
p = (int)sqrt((float)n + 0.5f);
for (int i = 1; i <= n; ++i) {
scanf("%d", a + i);
}
div[1][idx[1]++] = 1;
for (int i = 2; i <= n; ++i) {
lmt = (int)sqrt((float)i + 0.5f);
for (int j = 1; j <= lmt; ++j) {
if (i % j == 0) {
div[i][idx[i]++] = j;
if (j != i / j) {
div[i][idx[i]++] = i / j;
}
}
}
} scanf("%d", &q);
while (q--) {
scanf("%d", &opr);
if (opr == 1) {
scanf("%d", &t1);
ans = (long long)a[t1];
for (int i = 0; i < idx[t1]; ++i) {
t2 = div[t1][i];
ans += f[t2] + g[t2 / p];
}
printf("%I64d\n", ans);
}
else {
scanf("%d%d%I64d", &t1, &t2, &t3);
x = t1 / p;
y = t2 / p;
if (x == y) {
for (int i = t1; i <= t2; ++i) {
f[i] += t3;
}
continue;
}
for (int i = t1; i < (x + 1) * p; ++i) {
f[i] += t3;
}
for (int i = x + 1; i < y; ++i) {
g[i] += t3;
}
for (int i = y * p; i <= t2; ++i) {
f[i] += t3;
}
}
}
return 0;
}
[GZOI2016] 亚索的量子实验【分块】的更多相关文章
- 量子:基于ERP块对的两步量子直接通信
学习论文: 题目:Two-step quantum direct communication protocol using the Einstein-Podolsky-Rosen pair block ...
- Vertica集群扩容实验过程记录
需求: 将3个节点的Vertica集群扩容,额外增加3个节点,即扩展到6个节点的Vertica集群. 实验环境: RHEL 6.5 + Vertica 7.2.2-2 步骤: 1.三节点Vertica ...
- 【图像处理】第三次实验:JPEG图像压缩
1.任务说明 将LENA图像用JPEG方式压缩. 2.算法原理 JPEG(Joint Photographic Experts Group)是一个由ISO和IEC两个组织机构联合组成的一个专家组,负责 ...
- 第一部分:IBM量子体验
(一)量子世界 今天的计算机使用标准的(或用物理学家的话来说,“经典的”)计算模型来执行计算与处理信息,此计算模型要追溯到图灵(Turing)和冯・诺伊曼(Von Neumann)时期.在此模型 ...
- 20162318 实验三《 敏捷开发与XP实践》实验报告
北京电子科技学院(BESTI) 实 验 报 告 课程:程序设计与数据结构 班级:1623班 姓名:张泰毓 指导老师:娄老师.王老师 实验日期:2017年5月12日 实验密级:非密级 实验器材:带Lin ...
- 20172302 《Java软件结构与数据结构》实验三:查找与排序实验报告
课程:<Java软件结构与数据结构> 班级: 1723 姓名: 侯泽洋 学号:20172302 实验教师:王志强老师 实验日期:2018年11月19日 必修/选修: 必修 实验内容 (1) ...
- 20172309 《Java软件结构与数据结构》实验三报告
课程:<程序设计与数据结构(下)> 班级:1723 姓名: 王志伟 学号:20172309 实验教师:王志强老师 实验日期:2018年11月2日 必修/选修: 必修 实验内容: 实验一: ...
- Burpsuit分块传输插件绕WAF原理和技巧(转)
0x00 原理 给服务器发送payload数据包,使得waf无法识别出payload,当apache,tomcat等web容器能正常解析其内容.如图一所示 0x02 实验环境 本机win10+x ...
- 20172301 《Java软件结构与数据结构》实验三报告
20172301 <Java软件结构与数据结构>实验三报告 课程:<Java软件结构与数据结构> 班级: 1723 姓名: 郭恺 学号:20172301 实验教师:王志强老师 ...
随机推荐
- javascript array-like object
http://www.nfriedly.com/techblog/2009/06/advanced-javascript-objects-arrays-and-array-like-objects/ ...
- Hadoop在window上运行 user=Administrator, access=WRITE, inode="hadoop"
win7下eclipse中错误的详细描述如下: org.apache.hadoop.security.AccessControlException: org.apache.hadoop.securit ...
- Visual Studio Visual assistant注释也做拼写检查怎么办
1 打开Visual Assistant 2 在Advanced中找到Underlines,取消勾选"Underline spelling errors in comments and ...
- android事件分发(二)
非常早之前写过一篇android事件分发的博客,主要写的是它是怎样分发的,具体非常多原理的东西都没有涉及到.今天就从源代码看android怎样控制它的分发机制. 鉴于手机屏幕的限制,所以android ...
- 理解Android ANR的触发原理(转)
一.概述 ANR(Application Not responding),是指应用程序未响应,Android系统对于一些事件需要在一定的时间范围内完成,如果超过预定时间能未能得到有效响应或者响应时间过 ...
- HBase单机环境搭建
在搭建HBase单机环境之前,首先你要保证你已经搭建好Java环境: $ java -version java version "1.8.0_51" Java(TM) SE Run ...
- Linux Grub系统加密、破密、修复
一.在重新启动系统时候按任意键进入 grub界面 Grub加密 一.title前的密码 修改grub.conf 这种加密只是在用户要进入grub界面的时候提示要输入密码,但是可以正常进入系统,有没有 ...
- Redis和Memcache性能测试对比
Redis和Memcache在写入性能上面差别不大,读取性能上面尤其是批量读取性能上面Memcache全面胜出,当然Redis也有自己的优点:比如数据持久化.支持更多的数据结构(Set List ZS ...
- RabbitMQ常用命令、管理界面
1.运行CMD,cd切换到RabbitMQ安装目录sbin下E:\Program Files\RabbitMQ Server\rabbitmq_server-3.7.2\sbin 执行 rabbitm ...
- CAS和ABA
1 CAS compare and swap的缩写,详见乐观锁和悲观锁. 2 ABA 就是说,我获取的旧值是A,然后被人修改成了B,但是又被人修改成了A,我就认为并没有修改,更新内存. 解决办法,给每 ...