带 sin, cos 的线段树 - 牛客
链接:https://www.nowcoder.com/acm/contest/160/D
来源:牛客网
题目描述
给出一个长度为n的整数序列a1,a2,...,an,进行m次操作,操作分为两类。
操作1:给出l,r,v,将al,al+1,...,ar分别加上v;
操作2:给出l,r,询问
输入描述:
第一行一个整数n
接下来一行n个整数表示a1,a2,...,an
接下来一行一个整数m
接下来m行,每行表示一个操作,操作1表示为1 l r v,操作2表示为2 l r
保证1≤n,m,ai,v≤200000;1≤l≤r≤n,v是整数
输出描述:
对每个操作2,输出一行,表示答案,四舍五入保留一位小数
保证答案的绝对值大于0.1,且答案的准确值的小数点后第二位不是4或5
数据随机生成(n,m人工指定,其余整数在数据范围内均匀选取),并去除不满足条件的操作2
输入
4
1 2 3 4
5
2 2 4
1 1 3 1
2 2 4
1 2 4 2
2 1 3
输出
0.3
-1.4
-0.3
题意 : 给你 n 个数字,第一种操作是将一个区间内每一个数字加上同一个数字,第二种操作是求一个区间内每一个数 sin 的累加和
思路分析 :对于每个区间维护一下 cos 和 sin 的值,当一个区间要加上一个数字时,此时再重新计算 sin的值时 , sin(a + x) = sin(a)*cos(x) + cos(a)*sin(x) ,一个区间内的所有值都可以这样计算,因此就会用到区间内的 sin 总和 以及 cos 的总和
这个题有个很坑的地方,就是大量的地方用到 sin 与 cos 函数,若输入的是一个整形数强制转变为浮点数后再用 sin函数,cos函数则会超时,不强制转换会快一半的时间!!!
代码示例 :
#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn = 2e5+5;
#define lson k<<1
#define rson k<<1|1 int n, m;
struct node
{
int l, r;
ll lazy;
double rs, rc;
}t[maxn<<2]; void pushdown(int k){
double x1 = t[lson].rs*cos(t[k].lazy) + t[lson].rc*sin(t[k].lazy);
double x2 = t[lson].rc*cos(t[k].lazy) - t[lson].rs*sin(t[k].lazy);
t[lson].rs = x1, t[lson].rc = x2; x1 = t[rson].rs*cos(t[k].lazy) + t[rson].rc*sin(t[k].lazy);
x2 = t[rson].rc*cos(t[k].lazy) - t[rson].rs*sin(t[k].lazy);
t[rson].rs = x1, t[rson].rc = x2; t[lson].lazy += t[k].lazy;
t[rson].lazy += t[k].lazy;
t[k].lazy = 0;
}
int x;
void build(int l, int r, int k){
t[k].l = l, t[k].r = r;
t[k].rc = t[k].rs = 0.0;
t[k].lazy = 0;
if (l == r) {
scanf("%d", &x);
t[k].rs = sin(x), t[k].rc = cos(x);
return;
}
int m = (l+r) >> 1;
build(l, m, lson);
build(m+1, r, rson); t[k].rs = t[lson].rs+t[rson].rs;
t[k].rc = t[lson].rc+t[rson].rc;
} void update(int l, int r, ll v, int k){
if (l <= t[k].l && t[k].r <= r){
double x1 = t[k].rs*cos(v)+t[k].rc*sin(v);
double x2 = t[k].rc*cos(v)-t[k].rs*sin(v);
t[k].rs = x1, t[k].rc = x2;
t[k].lazy += v;
return;
}
if (t[k].lazy) pushdown(k);
int m = (t[k].l+t[k].r) >> 1;
if (l <= m) update(l, r, v, lson);
if (r > m) update(l, r, v, rson);
t[k].rs = t[lson].rs+t[rson].rs;
t[k].rc = t[lson].rc+t[rson].rc;
} double sum;
void query(int l, int r, int k){
if (l <= t[k].l && t[k].r <= r){
sum += t[k].rs;
return;
}
if (t[k].lazy) pushdown(k);
int m = (t[k].l + t[k].r) >> 1;
if (l <= m) query(l, r, lson);
if (r > m) query(l, r, rson);
} int main() {
int pt, l, r;
ll v; cin >> n;
build(1, n, 1);
cin >> m;
while(m--){
scanf("%d%d%d", &pt, &l, &r);
if (pt == 1) {
scanf("%lld", &v);
update(l, r, v, 1);
}
else {
sum = 0;
query(l, r, 1);
printf("%.1lf\n", sum);
}
}
return 0;
}
带 sin, cos 的线段树 - 牛客的更多相关文章
- BZOJ 4129 树上带修莫队+线段树
思路: 可以先做做BZOJ3585 是序列上的mex 考虑莫队的转移 如果当前数字出现过 线段树上把它置成1 对于询问 二分ans 线段树上查 0到ans的和 是不是ans+1 本题就是把它搞到了序列 ...
- 笛卡尔树--牛客第四场(sequence)
思路: O(n)建一颗笛卡尔树,再O(n)dfs向上合并答案就行了. #define IOS ios_base::sync_with_stdio(0); cin.tie(0); #include &l ...
- 线段树(单标记+离散化+扫描线+双标记)+zkw线段树+权值线段树+主席树及一些例题
“队列进出图上的方向 线段树区间修改求出总量 可持久留下的迹象 我们 俯身欣赏” ----<膜你抄> 线段树很早就会写了,但一直没有总结,所以偶尔重写又会懵逼,所以还是要总结一下. ...
- BZOJ 3878 [AHOI&JSOI2014]奇怪的计算器 (线段树)
题面:BZOJ传送门 洛谷传送门 线段树好题 题目保证$a$一定是正整数,容易发现计算结果是单调的 我们把询问离线,并按照从小到大排序 某次操作可能导致某些位置达到边界$L/R$ 根据单调性的结论 这 ...
- 「ZJOI2019」线段树
传送门 Description 线段树的核心是懒标记,下面是一个带懒标记的线段树的伪代码,其中 tag 数组为懒标记: 其中函数\(Lson(Node)\)表示\(Node\)的左儿子,\(Rson( ...
- @loj - 3043@「ZJOI2019」线段树
目录 @description@ @solution@ @accepted code@ @details@ @description@ 九条可怜是一个喜欢数据结构的女孩子,在常见的数据结构中,可怜最喜 ...
- 2019牛客第八场多校 E_Explorer 可撤销并查集(栈)+线段树
目录 题意: 分析: @(2019牛客暑期多校训练营(第八场)E_Explorer) 题意: 链接 题目类似:CF366D,Gym101652T 本题给你\(n(100000)\)个点\(m(1000 ...
- 牛客练习赛28 B数据结构(线段树)
链接:https://www.nowcoder.com/acm/contest/200/B来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言5242 ...
- 牛客练习赛28-B(线段树,区间更新)
牛客练习赛28 - B 传送门 题目 qn姐姐最好了~ qn姐姐给你了一个长度为n的序列还有m次操作让你玩, 1 l r 询问区间[l,r]内的元素和 2 l r 询问区间[l,r]内的 ...
随机推荐
- VMware 注册码
VMware 12 Pro 永久许可证激活密钥 5A02H-AU243-TZJ49-GTC7K-3C61NVF5XA-FNDDJ-085GZ-4NXZ9-N20E6UC5MR-8NE16-H81WY- ...
- SVN:符号
问号:未纳入SVN控制之下的新文件感叹号:已经修改过还没上传的文件红叉:准备删除尚未提交绿勾:从服务器上取下来没修改过的文件黄色感叹号:发生冲突蓝色加号:准备加到服务器上还没提交
- window 系统下修改`CMD`的编码格式的方法,`CHCP` 的 使用
CHCP的使用 CHCP是一个计算机指令,能够显示或设置活动代码页编号. 一般上是在命令提示框中使用,用来查询和修改命令提示框的编码格式 具体使用方法 查看活动代码页编号 方式1: >>& ...
- H3C DHCP中继显示及维护
- java 一个类加载器的高级问题分析
编写一个能打印出自己的类加载器名称和当前类加载器的父子结构关系链的MyServlet,正常发布后,看到打印结果为WebAppClassloader. 把MyServlet.class文件打jar包,放 ...
- 2019-9-22-dotnet-core-导出-COM-组件
title author date CreateTime categories dotnet core 导出 COM 组件 lindexi 2019-09-22 20:25:38 +0800 2019 ...
- JUnit 单元测试断言推荐 AssertJ
文章转自:http://sgq0085.iteye.com/blog/2030609 前言 由于JUnit的Assert是公认的烂API,所以不推荐使用,目前推荐使用的是AssertJ. Assert ...
- <QluOJ2018NewCode>计算几何(寄蒜几盒)
题目描述 现在有一个圆圈,圆圈上有若干个点,请判断能否在若干个点中选择三个点两两相连组成一个等边三角形? 这若干个点在圆圈上按顺时针顺序分布. 如果可以的话输出"Yes"(不含引号 ...
- 【34.54%】【codeforces 675E】Trains and Statistic
time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard o ...
- HDU - 4289 Control (Dinic)
You, the head of Department of Security, recently received a top-secret information that a group of ...