分析:

数据范围表示:c特别的小(c<20)

我们可以考虑nlogn*c^2的算法。

线段树维护区间信息:f[i]表示在[l,r]这段区间中选择i个数相乘的和。

因此,我们可以将区间看成一个点,在PushUp的时候用背包的方式更新父节点。(仔细观察发现这是卷积)

剩下的就是一些优化了...

附上代码:

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <iostream>
using namespace std;
#define N 50005
#define mod 19940417
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define ll long long
int C[N][21],a[N],n,Q;
char s[10];
struct node
{
ll f[21],tag,add;
int siz;
node ()
{
memset(f,0,sizeof(f));
siz=tag=add=0;
}
node operator +(const node &b)
{
node c;
c.f[0]=1;
for(int i=1;i<=20;i++)
{
for(int j=0;j<=i;j++)
{
c.f[i]+=f[j]*b.f[i-j]%mod;
c.f[i]%=mod;
}
}
c.siz=siz+b.siz;
return c;
}
void plus(ll x)
{
add+=x;
add%=mod;
for(int i=min(siz,20);i;i--)
{
ll y=x;
for(int j=1;j<=i;j++)
{
f[i]=(f[i]+y*f[i-j]%mod*C[siz-i+j][j]%mod)%mod;
y=y*x%mod;
}
}
}
void rev()
{
tag^=1;
add=(mod-add)%mod;
for(int i=min(siz,20);i;i--)
{
if(i&1)f[i]=(mod-f[i])%mod;
}
}
}tr[N<<2];
void PushUp(int rt)
{
tr[rt]=tr[rt<<1]+tr[rt<<1|1];
}
void build(int l,int r,int rt)
{
if(l==r)
{
tr[rt].f[1]=a[l];
tr[rt].f[0]=tr[rt].siz=1;
return ;
}
int m=(l+r)>>1;
build(lson);
build(rson);
PushUp(rt);
}
void PushDown(int rt)
{
if(tr[rt].tag)
{
tr[rt<<1].rev();
tr[rt<<1|1].rev();
tr[rt].tag=0;
}
if(tr[rt].add)
{
tr[rt<<1].plus(tr[rt].add);
tr[rt<<1|1].plus(tr[rt].add);
tr[rt].add=0;
}
}
void Update(int L,int R,int c,int l,int r,int rt)
{
if(L<=l&&r<=R)
{
tr[rt].plus(c);
return ;
}
PushDown(rt);
int m=(l+r)>>1;
if(m>=L)Update(L,R,c,lson);
if(m<R)Update(L,R,c,rson);
PushUp(rt);
}
void Update_rev(int L,int R,int l,int r,int rt)
{
if(L<=l&&r<=R)
{
tr[rt].rev();
return ;
}
PushDown(rt);
int m=(l+r)>>1;
if(m>=L)Update_rev(L,R,lson);
if(m<R)Update_rev(L,R,rson);
PushUp(rt);
}
node query(int L,int R,int l,int r,int rt)
{
if(L<=l&&r<=R)
{
return tr[rt];
}
PushDown(rt);
int m=(l+r)>>1;
if(m>=R)return query(L,R,lson);
if(m<L)return query(L,R,rson);
return query(L,R,lson)+query(L,R,rson);
}
int main()
{
scanf("%d%d",&n,&Q);
C[0][0]=1;
for(int i=1;i<=n;i++)
{
C[i][0]=1;
for(int j=1;j<=20;j++)
{
C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
}
}
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
build(1,n,1);
while(Q--)
{
int x,y,z;
scanf("%s%d%d",s,&x,&y);
if(s[0]=='I')
{
scanf("%d",&z);
Update(x,y,z,1,n,1);
}else if(s[0]=='R')
{
Update_rev(x,y,1,n,1);
}else
{
scanf("%d",&z);
printf("%lld\n",(query(x,y,1,n,1).f[z]+mod)%mod);
}
}
return 0;
}

  

最新文章

  1. MySQL: @variable vs. variable. Whats the difference?
  2. 金浪网关ESR-6400G 告警解决办法
  3. C++中的复制构造函数
  4. Codeforces Round #237 (Div. 2) B题模拟题
  5. 一些有用的javascript实例分析(一)
  6. shell-3
  7. ActiveMQ入门练习
  8. 19.C++-(=)赋值操作符、智能指针编写(详解)
  9. SSIS获得Excel行号(转自http://blog.csdn.net/zplume/article/details/19113911)
  10. loglog 函数的使用
  11. Visual C++ 6.0中关于for的简单问题
  12. python使用tcp实现一个简单的下载器
  13. 8个爽滑如丝的Windows小软件,不好用你拿王思葱砸死我
  14. PTA——时间转换
  15. python之绘制图形库turtle
  16. 基于 Laravel 的 文件管理
  17. Nginx+tomcat组合实现高并发场景的动静分离和负载均衡方案
  18. Git 环境设置(安装)
  19. Flask请求处理流程(request)[待续]
  20. firedac二进制序列和JSON序列的对比

热门文章

  1. web集群时session同步的3种方法[转]
  2. mysql统计类似SQL语句查询次数
  3. Nodejs经验谈
  4. Day12 前端html
  5. How 5 Natural Language Processing APIs Stack Up
  6. classes目录中没有class文件的一个原因
  7. 从JavaWeb危险字符过滤浅谈ESAPI使用
  8. jquery文本框内容实时监控
  9. leetCode刷题(找出数组里的两项相加等于定值)
  10. my views--软件工程、python