BZOJ 2962
2962: 序列操作
Time Limit: 50 Sec Memory Limit: 256 MB
Submit: 618 Solved: 225
[Submit][Status][Discuss]
Description
有一个长度为n的序列,有三个操作1.I a b c表示将[a,b]这一段区间的元素集体增加c,2.R a b表示将[a,b]区间内所有元素变成相反数,3.Q a b c表示询问[a,b]这一段区间中选择c个数相乘的所有方案的和mod 19940417的值。
Input
第一行两个数n,q表示序列长度和操作个数。
第二行n个非负整数,表示序列。
接下来q行每行输入一个操作I a b c或者 R a b或者Q a b c意义如题目描述。
Output
对于每个询问,输出选出c个数相乘的所有方案的和mod19940417的值。
Sample Input
1 2 3 4 5
I 2 3 1
Q 2 4 2
R 1 5
I 1 3 -1
Q 1 5 1
Sample Output
19940397
样例说明
做完第一个操作序列变为1 3 4 4 5。
第一次询问结果为3*4+3*4+4*4=40。
做完R操作变成-1 -3 -4 -4 -5。
做完I操作变为-2 -4 -5 -4 -5。
第二次询问结果为-2-4-5-4-5=-20。
HINT
100%的数据n<=50000,q<=50000,初始序列的元素的绝对值<=109,I a b c中保证[a,b]是一个合法区间,|c|<=109,R a b保证[a,b]是个合法的区间。Q a b c中保证[a,b]是个合法的区间1<=c<=min(b-a+1,20)。
Source

#include<cstdio>
#include<cstring>
#include<iostream>
#define N 60050
#define mod 19940417LL
#define ll long long
#define lc k<<1
#define rc k<<1|1
#define mid (l+r>>1)
using namespace std;
struct use{
ll fi[],delta,rev;
void init(){memset(fi,,sizeof(fi));delta=rev=;}
}tr[N<<];
ll ai[N],cc[N][],mi[N];
inline char letter(){
for(register char ch=getchar();;ch=getchar()) if(ch>='A'&&ch<='Z') return ch;
}
inline ll read(){
register ll x=,f=;
register char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
return x*f;
}
use updata(use x,use y){
use ci;ci.init();
for(int i=;i<=;i++){
ci.fi[i]=(x.fi[i]+y.fi[i])%mod;
for(int j=;j<i;j++)
ci.fi[i]=(ci.fi[i]+x.fi[j]*y.fi[i-j]%mod)%mod;
}
return ci;
}
void build(int k,int l,int r){
if(l==r){
tr[k].fi[]=ai[l];
return ;
}
build(lc,l,mid);build(rc,mid+,r);
tr[k]=updata(tr[lc],tr[rc]);
}
void paint(int k,int l,int r,ll x){
ll xx;
if(!x){
tr[k].rev^=;
tr[k].delta=((mod-tr[k].delta)%mod+mod)%mod;
for(int j=;j<=min(r-l+,);j+=) tr[k].fi[j]=((mod-tr[k].fi[j])%mod+mod)%mod;
}
else{
x=(x%mod+mod)%mod;
tr[k].delta=((tr[k].delta+x)%mod+mod)%mod;
for(int q,j=min(r-l+,);j;j--){
for(xx=x,q=j-;q;q--){
tr[k].fi[j]=(tr[k].fi[j]+cc[r-l+-q][j-q]*xx%mod*tr[k].fi[q]%mod)%mod;
xx=x*xx%mod;
}
tr[k].fi[j]=(tr[k].fi[j]+cc[r-l+][j]*xx%mod)%mod;
} }
}
void pushdown(int k,int l,int r){
if(tr[k].rev) paint(lc,l,mid,),paint(rc,mid+,r,);
if(tr[k].delta) paint(lc,l,mid,tr[k].delta),paint(rc,mid+,r,tr[k].delta);
tr[k].rev=tr[k].delta=;
}
void ins(int k,int l,int r,int x,int y,int v){
if(x<=l&&r<=y){
paint(k,l,r,v);
return ;
}
pushdown(k,l,r);
if(x<=mid) ins(lc,l,mid,x,y,v);
if(y>mid) ins(rc,mid+,r,x,y,v);
tr[k]=updata(tr[lc],tr[rc]);
}
void rev(int k,int l,int r,int x,int y){
if(x<=l&&r<=y){
paint(k,l,r,);
return ;
}
pushdown(k,l,r);
if(x<=mid) rev(lc,l,mid,x,y);
if(y>mid) rev(rc,mid+,r,x,y);
tr[k]=updata(tr[lc],tr[rc]);
}
use query(int k,int l,int r,int x,int y){
if(x<=l&&r<=y) return tr[k];
pushdown(k,l,r);
use x1,x2;x1.init();x2.init();
bool f1(),f2();
if(x<=mid) x1=query(lc,l,mid,x,y),f1=;
if(y>mid) x2=query(rc,mid+,r,x,y),f2=;
if(f1&&f2) return updata(x1,x2);
return f1?x1:x2;
}
int main(){
ll n,q,a,b,c;char ch;
n=read();q=read();
for(int i=;i<=n;i++) cc[i][]=;
for(int i=;i<=n;i++){
for(int j=;j<=;j++){
cc[i][j]=(cc[i-][j-]+cc[i-][j])%mod;
}
}
for(int i=;i<=n;i++){
ai[i]=read();ai[i]=(ai[i]%mod+mod)%mod;
}
build(,,n);
for(int i=;i<=q;i++){
if((ch=letter())=='I'){
a=read();b=read();c=read();
c=(c%mod+mod)%mod;
if(!c) continue;
ins(,,n,a,b,c);
}
else if(ch=='R'){
a=read();b=read();
rev(,,n,a,b);
}
else{
a=read();b=read();c=read();
use ci=query(,,n,a,b);
printf("%lld\n",ci.fi[c]);
}
}
return ;
}
BZOJ 2962的更多相关文章
- [bzoj]2962序列操作
[bzoj]2962序列操作 标签: 线段树 题目链接 题意 给你一串序列,要你维护三个操作: 1.区间加法 2.区间取相反数 3.区间内任意选k个数相乘的积 题解 第三个操作看起来一脸懵逼啊. 其实 ...
- bzoj 2962 序列操作
2962: 序列操作 Time Limit: 50 Sec Memory Limit: 256 MB[Submit][Status][Discuss] Description 有一个长度为n的序列, ...
- bzoj 2962 序列操作——线段树(卷积?)
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2962 如果 _,_,_,…… 变成了 (_+k),(_+k),(_+k),…… ,计算就是在 ...
- bzoj 2962 序列操作 —— 线段树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2962 维护 sum[i] 表示选 i 个的乘积和,合并两个子树就枚举两边选多少,乘起来即可: ...
- BZOJ 2127: happiness [最小割]
2127: happiness Time Limit: 51 Sec Memory Limit: 259 MBSubmit: 1815 Solved: 878[Submit][Status][Di ...
- BZOJ 3275: Number
3275: Number Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 874 Solved: 371[Submit][Status][Discus ...
- BZOJ 2879: [Noi2012]美食节
2879: [Noi2012]美食节 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1834 Solved: 969[Submit][Status] ...
- bzoj 4610 Ceiling Functi
bzoj 4610 Ceiling Functi Description bzoj上的描述有问题 给出\(n\)个长度为\(k\)的数列,将每个数列构成一个二叉搜索树,问有多少颗形态不同的树. Inp ...
- BZOJ 题目整理
bzoj 500题纪念 总结一发题目吧,挑几道题整理一下,(方便拖板子) 1039:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...
随机推荐
- thymeleaf中的Literals
Literals即为文字 一.Text literals:文本文字 文本文字只是字符串指定的单引号之间.他们可以包含任何字符,但你应避免任何单引号里面\ ' <p> Now you are ...
- hibernate[版本四]知识总结
1.hibernate是orm对象关系映射,是对jdbc的封装 2.hibernate版helloworld 2.1导入jar <dependencies> <dependency& ...
- C#全角半角转换函数
Code#region 全角半角转换 /// <summary> /// 转全角的函数(SBC case) /// </summary> /// <param name= ...
- Oracle超出最大连接数问题及解决
用过Oracle的应该都熟悉如何查看和设置Oracle数据库的最大连接数.这里就再啰嗦一遍. 查看当前的连接数,可以用select count(*) from v$process;设置的最大连接数(默 ...
- Oracle DataGuard搭建(二)
三.配置备库 创建catalog数据库 用dbca创建数据库,用oracle自带模板,不用em,启用归档,同一管理密码oracle,global name:dbcat1.node249.gewara, ...
- 又来折腾Linux
硬盘坏了之后就没装过Linux了,因为弄了一个很老的台式机的80G并口硬盘,根本不够用的,一直懒得理. 前段时间实验室的老机子得报销了,但是里面的东西还可以拆下来,所以又拆下了两个硬盘,这样就有三个8 ...
- visual studio 2012进行C语言开发[图文]
现在大家计算机大概都脱离XP了,so,之前蛮多可以用的编译器,可能放在我们现在的Win7,win8下面会出现一些比如不兼容了之类的问题.其实,用微软强大的IDE-visual studio系列,也是可 ...
- Android学习笔记(20)————利用ListView制作带竖线的多彩表格
http://blog.csdn.net/conowen/article/details/7421805 /********************************************** ...
- dubbo-RPC学习(二) http protocol
dubbo中的协议默认支持很多种,从简单的开始,先解读HttpProtocol吧. dubbo中的http-RPC基于spring web支持的http invoker,web容器默认使用jetty. ...
- 网络IPC:套接字之建立连接
如果处理的是面向连接的网络服务(SOCK_STREAM或SOCK_SEQPACKET),在开始交换数据以前,需要在请求服务的进程套接字(客户端)和提供服务的进程套接字(服务器)之间建立一个连接.客户端 ...