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

5 5
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

40
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

 
题解:
线段树训练的好题(纯线段树维护的代码试炼),这题自己敲出来,就说明你的线段树方面的代码能力已经过关了。
 
AC代码:
#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的更多相关文章

  1. [bzoj]2962序列操作

    [bzoj]2962序列操作 标签: 线段树 题目链接 题意 给你一串序列,要你维护三个操作: 1.区间加法 2.区间取相反数 3.区间内任意选k个数相乘的积 题解 第三个操作看起来一脸懵逼啊. 其实 ...

  2. bzoj 2962 序列操作

    2962: 序列操作 Time Limit: 50 Sec  Memory Limit: 256 MB[Submit][Status][Discuss] Description 有一个长度为n的序列, ...

  3. bzoj 2962 序列操作——线段树(卷积?)

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2962 如果 _,_,_,…… 变成了 (_+k),(_+k),(_+k),…… ,计算就是在 ...

  4. bzoj 2962 序列操作 —— 线段树

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2962 维护 sum[i] 表示选 i 个的乘积和,合并两个子树就枚举两边选多少,乘起来即可: ...

  5. BZOJ 2127: happiness [最小割]

    2127: happiness Time Limit: 51 Sec  Memory Limit: 259 MBSubmit: 1815  Solved: 878[Submit][Status][Di ...

  6. BZOJ 3275: Number

    3275: Number Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 874  Solved: 371[Submit][Status][Discus ...

  7. BZOJ 2879: [Noi2012]美食节

    2879: [Noi2012]美食节 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 1834  Solved: 969[Submit][Status] ...

  8. bzoj 4610 Ceiling Functi

    bzoj 4610 Ceiling Functi Description bzoj上的描述有问题 给出\(n\)个长度为\(k\)的数列,将每个数列构成一个二叉搜索树,问有多少颗形态不同的树. Inp ...

  9. BZOJ 题目整理

    bzoj 500题纪念 总结一发题目吧,挑几道题整理一下,(方便拖板子) 1039:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...

随机推荐

  1. poj 1847 Tram

    http://poj.org/problem?id=1847 这道题题意不太容易理解,n个车站,起点a,终点b:问从起点到终点需要转换开关的最少次数 开始的那个点不需要转换开关 数据: 3 2 1// ...

  2. spring注解使用

    一.各种注解方式 1.@Autowired注解(不推荐使用,建议使用@Resource) @Autowired可以对成员变量.方法和构造函数进行标注,来完成自动装配的工作.@Autowired的标注位 ...

  3. 转载 从Http到Https

    转载自 http://www.cnblogs.com/silin6/p/5928503.html HTTP 当你在浏览器输入一个网址 (例如 http://tasaid.com)的时候,浏览器发起一个 ...

  4. 解决iPhone上select时常失去焦点,随意跳到下一个输入框,影响用户操作

    window.addEventListener('load', function() { FastClick.attach(document.body); }, false); //300s延迟,解决 ...

  5. MSSQL索引优化

    转自:http://blog.itpub.net/16436858/viewspace-589275/ http://www.cnblogs.com/jams742003/archive/2011/1 ...

  6. 【转】WPF颜色相关操作

        using System.Windows.Media; 1.String转换成Color Color color = (Color)ColorConverter.ConvertFromStri ...

  7. php 建立类POST/GET 的HTTP请求

    1.第一种利用fsock的方式来建立类POST的请求. <?php $srv_ip = '192.168.1.5';//你的目标服务地址. $srv_port = 80;//端口 $url = ...

  8. 基础数据结构 之 队列(python实现)

    队也是编程开发中常见的一种数据结构.栈和队可用来模拟函数的递归过程.队的特点为先入先出,主要操作包括入队和出队.入队时需判断队是否已满,出队时需判断队是否为空.下面给出一个队的python实现的例子: ...

  9. 关于webapi 返回的类型的笔记

    经过测试发现使用IE浏览器返回的数据是json,而使用Firefox和Chrome返回的则为xml,经研究发现IE在发生http请求时请求头accpet节点相比Firefox和Chrome缺少&quo ...

  10. 从零开始学android开发-通过WebService进行网络编程,使用工具类轻松实现

    相信大家在平常的开发中,对网络的操作用到HTTP协议比较多,通过我们使用Get或者Post的方法调用一个数据接口,然后服务器给我们返回JSON格式的数据,我们解析JSON数据然后展现给用户,相信很多人 ...