1798: [Ahoi2009]Seq 维护序列seq
1798: [Ahoi2009]Seq 维护序列seq
Time Limit: 30 Sec Memory Limit: 64 MB
Submit: 2930 Solved: 1087
[Submit][Status]
Description
Input
Output
Sample Input
1 2 3 4 5 6 7
5
1 2 5 5
3 2 4
2 3 7 9
3 1 3
3 4 7
Sample Output
35
8
HINT
【样例说明】
初始时数列为(1,2,3,4,5,6,7)。
经过第1次操作后,数列为(1,10,15,20,25,6,7)。
对第2次操作,和为10+15+20=45,模43的结果是2。
经过第3次操作后,数列为(1,10,24,29,34,15,16}
对第4次操作,和为1+10+24=35,模43的结果是35。
对第5次操作,和为29+34+15+16=94,模43的结果是8。
测试数据规模如下表所示
数据编号 1 2 3 4 5 6 7 8 9 10
N= 10 1000 1000 10000 60000 70000 80000 90000 100000 100000
M= 10 1000 1000 10000 60000 70000 80000 90000 100000 100000
Source
题解:涉及到区间加法和乘法的线段树呵呵呵。。。虽然不难但是我已经被其折腾了一年了。。。。还好AC了
个人觉得重点在于对于一个线段,当需要改变这个线段的一部分时,这个线段上的总和值该如何维护的问题,直到今天我方才想出了一个机智的方法——既然直接维护比较讨厌,那么何不在下面方便维护的点直接维护,然后将各个子线段的变化量累计起来,不就是此线段的变化量了么。。。呵呵呵(phile:今天才明白? HansBug:TT)
type
vet=record
a0,a1:int64;
end;
var
i,j,k,l,m,n,a2,a3,a4:longint;
p:int64;
d1,d2,d:vet;
a,c:array[..] of int64;
b:array[..] of vet;
function min(x,y:longint):longint;inline;
begin
if x<y then min:=x else min:=y;
end;
function max(x,y:longint):longint;inline;
begin
if x>y then max:=x else max:=y;
end;
function merge(d1,d2:vet):vet;inline;
var d3:vet;
begin
d3:=d1;
d3.a0:=d3.a0 mod p;
d3.a1:=d3.a1 mod p;
d2.a0:=d2.a0 mod p;
d2.a1:=d2.a1 mod p;
d3.a0:=(d3.a0*d2.a0) mod p;
d3.a1:=((d3.a1*d2.a0) mod p+d2.a1) mod p;
exit(d3);
end;
procedure built(z,x,y:longint);inline;
begin
if (x=y) then
a[z]:=c[x] mod p
else
begin
built(z*,x,(x+y) div );
built(z*+,(x+y) div +,y);
a[z]:=(a[z*]+a[z*+]) mod p;
end;
b[z].a0:=;b[z].a1:=;
end;
procedure ext(z,x,y:longint);inline;
begin
a[z]:=((a[z]*b[z].a0) mod p+b[z].a1*(y-x+)) mod p;
b[z*]:=merge(b[z*],b[z]);
b[z*+]:=merge(b[z*+],b[z]);
b[z].a0:=;b[z].a1:=;
end;
function op(z,x,y,l,r:longint;d:vet):int64;inline;
var a2,a3,a4:int64;
begin
if l>r then exit();
ext(z,x,y);a2:=a[z];
if (x=l) and (y=r) then
begin
b[z]:=d;
exit((a2*(b[z].a0-)) mod p+(b[z].a1*(r-l+)) mod p);
end
else
begin
a3:=op(z*,x,(x+y) div ,l,min((x+y) div ,r),d) mod p;
a4:=op(z*+,(x+y) div +,y,max((x+y) div +,l),r,d) mod p;
a[z]:=(a[z]+(a3+a4) mod p) mod p;exit((a3+a4) mod p);
end;
end;
function sum(z,x,y,l,r:longint;d:vet):int64;inline;
var d1,d2:vet;
begin
if l>r then exit();
d1:=b[z];d1:=merge(d1,d);;
if (x=l) and (y=r) then
exit(((d1.a0*a[z]) mod p+(d1.a1*(r-l+)) mod p) mod p)
else
exit((sum(z*,x,(x+y) div ,l,min((x+y) div ,r),d1)+sum(z*+,(x+y) div +,y,max((x+y) div +,l),r,d1)) mod p);
end;
procedure showoff(z,x,y,l:longint);inline;
begin
writeln('':l*,z,'(',x,',',y,') = Tag=(',b[z].a0,',',b[z].a1,') ',a[z]);
if x<y then
begin
showoff(z*,x,(x+y) div ,l+);
showoff(z*+,(x+y) div +,y,l+);
end;
end;
begin
readln(n,p);
for i:= to n do read(c[i]);
readln;
built(,,n);
readln(m);
for i:= to m do
begin
read(j);
case j of
:begin
readln(a2,a3,a4);
d1.a0:=a4;d1.a1:=;
op(,,n,a2,a3,d1);
end;
:begin
readln(a2,a3,a4);
d1.a0:=;d1.a1:=a4;
op(,,n,a2,a3,d1);
end;
:begin
readln(a2,a3);
d1.a0:=;d1.a1:=;
writeln((sum(,,n,a2,a3,d1)+p) mod p);
end;
end;
end;
readln;
end.
1798: [Ahoi2009]Seq 维护序列seq的更多相关文章
- BZOJ 1798: [Ahoi2009]Seq 维护序列seq( 线段树 )
线段树.. 打个 mul , add 的标记就好了.. 这个速度好像还挺快的...( 相比我其他代码 = = ) 好像是#35.. ---------------------------------- ...
- bzoj 1798: [Ahoi2009]Seq 维护序列seq (线段树 ,多重标记下放)
1798: [Ahoi2009]Seq 维护序列seq Time Limit: 30 Sec Memory Limit: 64 MBSubmit: 7773 Solved: 2792[Submit ...
- bzoj 1798: [Ahoi2009]Seq 维护序列seq 线段树 区间乘法区间加法 区间求和
1798: [Ahoi2009]Seq 维护序列seq Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeO ...
- Bzoj 1798: [Ahoi2009]Seq 维护序列seq(线段树区间操作)
1798: [Ahoi2009]Seq 维护序列seq Time Limit: 30 Sec Memory Limit: 64 MB Description 老师交给小可可一个维护数列的任务,现在小可 ...
- BZOJ1798: [Ahoi2009]Seq 维护序列seq[线段树]
1798: [Ahoi2009]Seq 维护序列seq Time Limit: 30 Sec Memory Limit: 64 MBSubmit: 5504 Solved: 1937[Submit ...
- 【BZOJ】1798: [Ahoi2009]Seq 维护序列seq
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1798 大概就是维护两个标记的线段树模板题. 设定优先级,先乘后加(只是相对的),$push ...
- 【BZOJ】1798: [Ahoi2009]Seq 维护序列seq(线段树)
http://www.lydsy.com/JudgeOnline/problem.php?id=1798 之前写了个快速乘..........................20多s...... 还好 ...
- bzoj 1798 [Ahoi2009]Seq 维护序列seq(线段树+传标)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1798 [题意] 给定一个序列,要求提供区间乘/加,以及区间求和的操作 [思路] 线段树 ...
- AC日记——[Ahoi2009]Seq 维护序列seq bzoj 1798
1798 思路: 维护两个标记: 乘:m 和 加:a 先下放乘,再下放加: 下放乘的时候要把子节点的加一块乘了: 开long long: 来,上代码: #include <cstdio> ...
随机推荐
- LINQ to Sql系列一 增,删,改
首先,我们来了解一下LINQ to sql 是什么,我想从以下几点来做说明. 1,LINQ是从.net framework 3.0开始引入的的语言特性,主要是用来对集合数据进行操作的.2,LINQ t ...
- Android学习笔记总结
第一步: Android(1) - 在 Windows 下搭建 Android 开发环境,以及 Hello World 程序 搭建 Android 的开发环境,以及写一个简单的示例程序 · 在 Win ...
- Sql Server 查询多行并一行
干货 CREATE TABLE #benefit_code21 (id INT, number nvarchar(MAX), pname ), collegeID INT, applicationda ...
- JS插件库
http://www.cnblogs.com/xiaoyao2011/category/327551.html
- bootstrap 按钮 文本 浮动 隐藏
bootstrap 按钮 文本 浮动 隐藏 <!DOCTYPE html> <html lang="en"> <head> <meta c ...
- 深入了解Android中的AsyncTask
AsyncTask,即异步任务,是Android给我们提供的一个处理异步任务的类.通过此类,可以实现UI线程和后台线程进行通讯,后台线程执行异步任务,并把结果返回给UI线程. 我们知道,Androi ...
- zoj3823--构造
题目大意: 在n*n(n<=512)的网格上,从边界某个点出发,经过每个点一次且回到边界上,构造出一种方案使拐弯的数量至少为n*(n-1)-1次. 构造方法:我们可以手算出n=2~6时的方案. ...
- 抽象类 abstract 和 接口 interface 类的区别
在看一些框架的优秀改良时,都会设计一层base层,并且 base里面一般都是 abstract 类,然后 就找了为什么做的原因.发现: PHP5支持抽象类和抽象方法.抽象类不能直接被实例化,你必须先继 ...
- quagga源码学习--BGP协议路由更新
BGP的核心就是交换路由,所以关键的部分还是在路由的更新与撤销上面,这之间包含了冗长的属性,community等等处理过程,不做详述. bgp_read函数是路由更新的事件处理函数,在收到BGP_MS ...
- iOS核心笔记—源代码管理工具-GIT
源代码管理工具-GIT 一. git 概述 1. git 简介? 什么是git? > git是一款开源的分布式版本控制工具 > 在世界上所有的分布式版本控制工具中,git是最快.最简单.最 ...