[bzoj]2962序列操作

标签: 线段树


题目链接

题意

给你一串序列,要你维护三个操作:

1.区间加法

2.区间取相反数

3.区间内任意选k个数相乘的积

题解

第三个操作看起来一脸懵逼啊。

其实是可以合并的。

$ c[o].s[i]=\sum_{j=0}^{20}c[lc].s[j]×c[rc].s[i-j]\(
跟\)C_mn=\sum_{i=0}n C_ni×C_{m-n}{n-i} $这个等式是一个道理的吧。

然后想怎么维护加和取反。

取反比较容易,把取奇数个的答案变成相反数。

加法就稍微复杂一点。

假如p是加上的数,s[i]是区间内取k个数的答案。

\[s[i]=\sum_{j=0}^k C_{n-j}^i×s[j]×p^{k-j}
\]

Code

(不知道为什么常数很大,bzoj上43s才过)

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
#include<stack>
#include<set>
#include<map>
using namespace std;
#define ll long long
#define REP(i,a,b) for(int i=(a),_end_=(b);i<=_end_;i++)
#define DREP(i,a,b) for(int i=(a),_end_=(b);i>=_end_;i--)
#define EREP(i,a) for(int i=start[(a)];i;i=e[i].next)
inline int read()
{
int sum=0,p=1;char ch=getchar();
while(!(('0'<=ch && ch<='9') || ch=='-'))ch=getchar();
if(ch=='-')p=-1,ch=getchar();
while('0'<=ch && ch<='9')sum=sum*10+ch-48,ch=getchar();
return sum*p;
} const int mod=19940417;
const int maxn=1e5+20;
struct node {
int s[21];
int add,rev;
void Merge(node a,node b)
{
REP(i,0,20)
{
s[i]=0;
REP(j,0,i)s[i]=(s[i]+(ll)a.s[j]*b.s[i-j])%mod;
}
}
}; node c[maxn*4]; int n;
int a[maxn],C[maxn][25]; #define lc (o<<1)
#define rc (o<<1 | 1)
#define left lc,l,mid
#define right rc,mid+1,r void Reverse(int o,int l,int r)
{
c[o].rev^=1;
c[o].add*=-1;
REP(i,0,20)
if(i & 1)c[o].s[i]*=-1;
} void Add(int o,int l,int r,ll add)
{
c[o].add=(c[o].add+add)%mod;
DREP(i,min(20,r-l+1),0)
{
ll X=add;
DREP(j,i-1,0)
{
c[o].s[i]=(c[o].s[i]+(ll)c[o].s[j]*X%mod*C[r-l+1-j][i-j])%mod;
X=(ll)X*add%mod;
}
}
} void pushdown(int o,int l,int r)
{
int mid=(l+r)>>1;
if(c[o].rev)
{
c[o].rev=0;
Reverse(lc,l,mid);
Reverse(rc,mid+1,r);
}
if(c[o].add)
{
Add(lc,l,mid,c[o].add);
Add(rc,mid+1,r,c[o].add);
c[o].add=0;
}
} void make_tree(int o,int l,int r)
{
if(l==r)
{
c[o].s[0]=1;c[o].s[1]=a[l];
c[o].add=c[o].rev=0;
return;
}
int mid=(l+r)>>1;
make_tree(left);
make_tree(right);
c[o].Merge(c[lc],c[rc]);
} int q; void init()
{
n=read();q=read();
REP(i,1,n)a[i]=read();
C[0][0]=1;
REP(i,1,n)
{
C[i][0]=1;
REP(j,1,i)C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
}
make_tree(1,1,n);
} void updatea(int ql,int qr,ll x,int o,int l,int r)
{
if(ql<=l && r<=qr)
{
Add(o,l,r,x);
return;
}
int mid=(l+r)>>1;
pushdown(o,l,r);
if(ql<=mid )updatea(ql,qr,x,left);
if(qr>mid)updatea(ql,qr,x,right);
c[o].Merge(c[lc],c[rc]);
} void updater(int ql,int qr,int o,int l,int r)
{
if(ql<=l && r<=qr)
{
Reverse(o,l,r);
return;
}
int mid=(l+r)>>1;
pushdown(o,l,r);
if(ql<=mid)updater(ql,qr,left);
if(qr>mid)updater(ql,qr,right);
c[o].Merge(c[lc],c[rc]);
} node query(int ql,int qr,int o,int l,int r)
{
if(ql<=l && r<=qr)return c[o];
int mid=(l+r)>>1;
pushdown(o,l,r);
if(ql>mid)return query(ql,qr,right);
else if(qr<=mid)return query(ql,qr,left);
else
{
node a=query(ql,mid,left),b=query(mid+1,qr,right);
node c;c.Merge(a,b);
return c;
}
} void doing()
{
REP(i,1,q)
{
char ch;
//cin>>ch;
scanf("\n%c",&ch);
if(ch=='I')
{
int l=read(),r=read();ll x=read();
updatea(l,r,x,1,1,n);
}
else if(ch=='R')
{
int l=read(),r=read();
updater(l,r,1,1,n);
}else
{
int l=read(),r=read(),k=read();
node x=query(l,r,1,1,n);
printf("%d\n",(x.s[k]+mod)%mod);
}
}
} int main()
{
freopen("sequence.in","r",stdin);
freopen("sequence.out","w",stdout);
init();
doing();
return 0;
}

[bzoj]2962序列操作的更多相关文章

  1. bzoj 2962 序列操作

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

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

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

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

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

  4. bzoj 1858 序列操作

    bzoj 1858 序列操作 带有随机多个区间单值覆盖的区间操作题,可考虑用珂朵莉树解决. #include<bits/stdc++.h> using namespace std; #de ...

  5. 【BZOJ-2962】序列操作 线段树 + 区间卷积

    2962: 序列操作 Time Limit: 50 Sec  Memory Limit: 256 MBSubmit: 678  Solved: 246[Submit][Status][Discuss] ...

  6. BZOJ 2962

    2962: 序列操作 Time Limit: 50 Sec  Memory Limit: 256 MBSubmit: 618  Solved: 225[Submit][Status][Discuss] ...

  7. bzoj2962 序列操作

    2962: 序列操作 Time Limit: 50 Sec  Memory Limit: 256 MBSubmit: 1145  Solved: 378[Submit][Status][Discuss ...

  8. (WAWAWAWAWAWA) BZOJ 1858: [Scoi2010]序列操作

    二次联通门 : BZOJ 1858: [Scoi2010]序列操作 /* BZOJ 1858: [Scoi2010]序列操作 已经... 没有什么好怕的的了... 16K的代码... 调个MMP啊.. ...

  9. bzoj 1858: [Scoi2010]序列操作

    1858: [Scoi2010]序列操作 Time Limit: 10 Sec  Memory Limit: 64 MB 线段树,对于每个区间需要分别维护左右和中间的1和0连续个数,并在op=4时特殊 ...

随机推荐

  1. fsockopen — 打开一个网络连接或者一个Unix套接字连接

    fsockopen (PHP 4, PHP 5, PHP 7) 说明 resource fsockopen ( string $hostname [, int $port = -1 [, int &a ...

  2. HTML meta refresh 刷新与跳转(重定向)页面

    下面为各位整理了一些HTML meta refresh 刷新与跳转(重定向)页面的例子吧,后面本站长自己也补充了一些js页面刷新与跳转例子吧. refresh 属性值  --  刷新与跳转(重定向)页 ...

  3. 《并行程序设计导论》——MPI(Microsoft MPI)(1):Hello

    =================================版权声明================================= 版权声明:原创文章 禁止转载  请通过右侧公告中的“联系邮 ...

  4. iOS 注册页面 报错

    1.SpringBoard 无法启动应用程序(错误:-3) 方法:重启模拟器 2.

  5. vue中什么样的数据可以是在视图中显示

    1. Vue中不可以添加不存在的属性,因为不存在的属性是没有getter和setter的. <div id="app"> {{msg.a}} {{msg.b}} < ...

  6. python_tornado_session用户验证

    什么是session? -- Django中带有session,tornado中自己写 -- 逻辑整理 用户请求过来,验证通过,随机生成一个字符串当作value返回给浏览器, 在服务器中用户信息与随机 ...

  7. 函数式编程--使用lambda表达式

    前面一篇博客我们已经说到了,lambda表达式允许使用更简洁的代码来创建只有一个抽象方法的接口的实例.现在我们来写一段java的命令者模式来自己研究下lambda表达式的语法. 这里重复下命令者模式: ...

  8. linkin大话数据结构--Map

    Map 映射关系,也有人称为字典,Map集合里存在两组值,一组是key,一组是value.Map里的key不允许重复.通过key总能找到唯一的value与之对应.Map里的key集存储方式和对应的Se ...

  9. android dialog弹出的情况下监听返回键

    view = LayoutInflater.from(getActivity()).inflate( R.layout.dialog_tips, null); title2 = (TextView) ...

  10. Hibernate (三)

    1 一对多的单向 示例:一个已经存在的学生,新建一个班级,然后将该学生加入到该班级之下 设置inverse="false" <?xml version="1.0&q ...