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:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...
随机推荐
- C#Windows窗体界面设计_01_绘制三角函数_附强制类型转换
binzhouweichao@163.com 今天开始学习C#windows窗体界面设计. 首先说一下类型转换. 参考http://www.csharpwin.com/csharpspace/6848 ...
- ElasticSearch+Springboot实际应用:索引同步建设,搜索过程
1.介绍 springboot框架,众多自动化的部署和约定配置,造成了springboot的着手麻烦,熟练后可以快速快捷进行开发,常用作快捷开发的java底层框架.各位看官都是大神,自行体会. ...
- Python写的东西在CMD下打印中文
以前遇到的问题是在IDLE中可以print出中文,在cmd却是乱码. 后来想明白,IDLE中默认编码是unicode,而cmd中是ANSI,即cp963,即GBK 所以这样就能输出中文了: s = “ ...
- Codeforces Round #329 (Div. 2) A. 2Char 暴力
A. 2Char Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/contest/593/problem/A De ...
- 在Zend Studio中为ThinkPHP添加代码自动提示功能
身边很多朋友都使用ThinkPHP或CodeIgniter等开发框架为自己的项目提高开发效率. 在得力于这些优秀框架良好的设计结构的同时,也头疼于代码的自动完成提示功能没有纯PHP网站那么完善了.经常 ...
- iOS开发——网络编程Swift篇&(六)异步Post方式
异步Post方式 // MARK: - 异步Post方式 func asynchronousPost() { //创建NSURL对象 var url:NSURL! = NSURL(string: &q ...
- Metadata Lock原理1
https://www.percona.com/blog/2013/02/01/implications-of-metadata-locking-changes-in-mysql-5-5/ impli ...
- XCode帮助文档离线下载解决办法
1.菜单栏Xcode->Preferences选择Documentation,在线下载 2.离线下载(用迅雷即可下载) 在上述在线下载列表中,点击某一列,下拉框可看见 info,可得到其网络所在 ...
- 打开Excel时提示"您尝试打开的文件**.xls的格式与文件扩展名指定的格式不一致"
问题描述: 系统安装了WPS时,Analyzer导出excel时候,会提示"您尝试打开的文件**.xls的格式与文件扩展名指定的格式不一致",这是Excel的安全问题, ...
- centos安装zendstudio centos系统
查看centos系统32或64位命令(位数):在终端中执行“getconf LONG_BIT”命令并回车键确定 安装zend studio : 1. 首先安装好Java环境yum -y install ...