UVaLive6443(线段树)
要点
- 题意……题意往往是个大坎Orz:输入操作
p
则在区间\([x_1,x_2]\)插入一个三次函数,t
则先查询区间\([x_1,x_2]\)的函数值的和,然后按题目要求得到新的\(x_1\)、\(x_2\)并插入一个三次函数。都是整形。 - 解决方法:三次函数只有四个部分,栽4棵代表不同次幂函数的线段树即可,这样每次只维护某一特定次幂的函数是很容易做到的。
- 感觉这题很简单了对吧,然而每个月都要来那么一次红红的一板WA ……我调了很久才过的(菜鸡才会掉的)坑:
- 预处理的sum不应在询问答案时再使用,因为题目意思是没有函数的点函数值算0。所以要在线段树的sum的维护的过程中就使用这个数组。
- 一些手(智)生(障)的错误:把标准同步关了却还用printf、减法取模不加mod直接%mod、线段树更新写错了(这可还行……)。
- 休息并出门(调出bug的神奇做法)之后回来一眼看出\(x[i]\)是int范围的,要
x[i] = (x[i] % mod + mod) % mod;
,原来写的x[i] = (x[i] + mod) % mod
。
压着时限过的
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int X = 1e6 + 5, mod = 1e9 + 7;
const ll INF = 1e18;
int Test, N, l, r;
ll x[4];
ll sum[4][X];
char s[2];
struct SegmentTree {
#define ls(p) p << 1
#define rs(p) p << 1 | 1
struct Node {
int l, r;
int sum, tag;
}t[X << 2];
void Deal(int s, int p, int id) {
t[s].sum = ((ll)t[s].sum + (ll)t[p].tag * (sum[id][t[s].r] - sum[id][t[s].l - 1] + mod) % mod) % mod;
t[s].tag = ((ll)t[s].tag + (ll)t[p].tag) % mod;
}
void Push_down(int p, int id) {
if (t[p].tag) {
Deal(ls(p), p, id), Deal(rs(p), p, id);
t[p].tag = 0;
}
}
void Build(int l, int r, int p) {
t[p].l = l, t[p].r = r, t[p].tag = t[p].sum = 0;
if (l == r) return;
int mid = (l + r) >> 1;
Build(l, mid, ls(p));
Build(mid + 1, r, rs(p));
}
void Update(int l, int r, int p, ll k, int id) {
if (l <= t[p].l && t[p].r <= r) {
t[p].sum = ((ll)t[p].sum + k * (sum[id][t[p].r] - sum[id][t[p].l - 1] + mod) % mod) % mod;
t[p].tag = ((ll)t[p].tag + k) % mod;
return;
}
Push_down(p, id);
int mid = (t[p].l + t[p].r) >> 1;
if (l <= mid) Update(l, r, ls(p), k, id);
if (mid < r) Update(l, r, rs(p), k, id);
t[p].sum = ((ll)t[ls(p)].sum + t[rs(p)].sum) % mod;
}
ll Query(int l, int r, int p, int id) {
if (l <= t[p].l && t[p].r <= r) return t[p].sum;
Push_down(p, id);
int mid = (t[p].l + t[p].r) >> 1;
if (l > mid) return Query(l, r, rs(p), id);
if (r <= mid) return Query(l, r, ls(p), id);
return ((ll)Query(l, r, ls(p), id) + Query(l, r, rs(p), id)) % mod;
}
}T[4];
void Pre(int k) {
for (int i = 0; i <= X - 5; i++) {
ll t = 1LL;
for (int j = 1; j <= k; j++)
t = t * i % mod;
sum[k][i + 1] = (t + sum[k][i]) % mod;
}
}
int main() {
ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0);
for (int i = 0; i < 4; i++)
Pre(i);
cin >> Test;
for (int kase = 1; kase <= Test; kase++) {
cout << "Case #" << kase << ":\n";
for (int i = 0; i < 4; i++)
T[i].Build(1, X - 4, 1);
for (cin >> N; N; N--) {
cin >> s >> l >> r;
l++, r++;
for (int i = 3; ~i; --i)
cin >> x[i], x[i] = (x[i] % mod + mod) % mod;
if (s[0] == 'p') {
for (int i = 3; ~i; --i)
T[i].Update(l, r, 1, x[i], i);
} else {
ll ans = 0LL;
for (int i = 3; ~i; --i) {
ans = (ans + T[i].Query(l, r, 1, i)) % mod;
}
cout << ans << '\n';
l--, r--;
l = ans * l % 1000000;
r = ans * r % 1000000;
if (l > r) swap(l, r);
l++, r++;
for (int i = 3; ~i; --i)
T[i].Update(l, r, 1, x[i], i);
}
}
}
return 0;
}
UVaLive6443(线段树)的更多相关文章
- bzoj3932--可持久化线段树
题目大意: 最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分.超级计算机中的 任务用三元组(Si,Ei,Pi)描述,(Si,Ei,Pi)表示任务从第Si秒开始,在第 ...
- codevs 1082 线段树练习 3(区间维护)
codevs 1082 线段树练习 3 时间限制: 3 s 空间限制: 128000 KB 题目等级 : 大师 Master 题目描述 Description 给你N个数,有两种操作: 1:给区 ...
- codevs 1576 最长上升子序列的线段树优化
题目:codevs 1576 最长严格上升子序列 链接:http://codevs.cn/problem/1576/ 优化的地方是 1到i-1 中最大的 f[j]值,并且A[j]<A[i] .根 ...
- codevs 1080 线段树点修改
先来介绍一下线段树. 线段树是一个把线段,或者说一个区间储存在二叉树中.如图所示的就是一棵线段树,它维护一个区间的和. 蓝色数字的是线段树的节点在数组中的位置,它表示的区间已经在图上标出,它的值就是这 ...
- codevs 1082 线段树区间求和
codevs 1082 线段树练习3 链接:http://codevs.cn/problem/1082/ sumv是维护求和的线段树,addv是标记这歌节点所在区间还需要加上的值. 我的线段树写法在运 ...
- PYOJ 44. 【HNSDFZ2016 #6】可持久化线段树
#44. [HNSDFZ2016 #6]可持久化线段树 统计 描述 提交 自定义测试 题目描述 现有一序列 AA.您需要写一棵可持久化线段树,以实现如下操作: A v p x:对于版本v的序列,给 A ...
- CF719E(线段树+矩阵快速幂)
题意:给你一个数列a,a[i]表示斐波那契数列的下标为a[i],求区间对应斐波那契数列数字的和,还要求能够维护对区间内所有下标加d的操作 分析:线段树 线段树的每个节点表示(f[i],f[i-1])这 ...
- 【BZOJ-3779】重组病毒 LinkCutTree + 线段树 + DFS序
3779: 重组病毒 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 224 Solved: 95[Submit][Status][Discuss] ...
- 【BZOJ-3673&3674】可持久化并查集 可持久化线段树 + 并查集
3673: 可持久化并查集 by zky Time Limit: 5 Sec Memory Limit: 128 MBSubmit: 1878 Solved: 846[Submit][Status ...
随机推荐
- BZOJ 1614 [Usaco2007 Jan]Telephone Lines架设电话线:spfa + 二分【路径中最大边长最小】
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1614 题意: 给你一个无向图,n个点,m条边. 你需要找出一条从1到n的路径,使得这条路径 ...
- 如何改变cmd窗口大小
按下Win+R(windows徽标和R一起按),再输入cmd. 输入mode ,查看当前屏幕的大小,来为后面准备(大概估计数值大小). 再输入mode con:cols=400 lines=400 ...
- centos7搭建mysql-5.7.22主从复制
mysql7.7.22主从复制 本项目是根据真实环境搭建编写出文档,文档中的目录也是根据自己公司环境所创建.公司原来是一台服务器搭建的数据库(5.7.22),由于业务的扩展需要搭建一台从服务器,减轻主 ...
- OpenCV——PS 滤镜算法之平面坐标到极坐标的变换
// define head function #ifndef PS_ALGORITHM_H_INCLUDED #define PS_ALGORITHM_H_INCLUDED #include < ...
- 【Shell】基础正则表示法及grep用法
——<鸟哥的私房菜> 正规表示法就是处理字串的方法,他是以行为单位来进行字串的处理行为:正规表示法透过一些特殊符号的辅助,可以让使用者轻易的达到『搜寻/删除/取代』某特定字串的处理程序:只 ...
- Highcharts小数保留两位方法
tooltip: { valueSuffix: '%', formatter:function(){ return''+this.series.name+' '+Highcharts.numberFo ...
- Godot-3D教程-02.3D性能和局限性
介绍 Introduction Godot遵循表现与性能平衡信条.在这个表现的世界中,它们总是有许多约定俗成的东西,主要是在用执行速度换取可用性与扩展性方面.下面是一些实际的例子: 高效地渲染对象是比 ...
- 百度地图API的第一次接触
因为项目的需求,第一次接触了百度API. 第一步:引用百度地图API的脚本 如果在局域网环境中,要把地图文件和js文件都要下载下来 <script type="text/javascr ...
- HL7 V2 分隔符
Delimiter Characters Character Purpose 0x0D Marks the end of each segment. | Composite delimiter. ^ ...
- ZigBee自组网地址分配与路由协议概述
1. ZigBee简介 ZigBee是基于IEEE802.15.4标准的低功耗局域网协议.根据国际标准规定,ZigBee技术是一种短距离.低功耗的无线通信技术. ZigBee协议从下到上分别为物理层( ...