Luogu 3373(线段树标记混合下传)
2017-11-05 16:20
316 查看
传送门
题意:nlogn实现区间加、区间乘、区间求和(都要求取模)
题解:
只要记住一句话:(x+a)*b=x*b+a*b,先乘后加,乘法标记管两个(sum,add),加法标记管一个(sum)。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define root 1,1,n
#define lson rt<<1,l,mid
#define rson rt<<1|1,mid+1,r
typedef long long ll;
const int N=1e5+4;
int n,m;
ll a
,mod,sum[N<<2],add[N<<2],mul[N<<2];
inline ll read() {
ll x=0,f=1;char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
return x*f;
}
void build(int rt,int l,int r) {
add[rt]=0,mul[rt]=1;
if (l==r) {sum[rt]=a[l];return ;}
int mid=l+r>>1;
build(lson);
build(rson);
(sum[rt]=sum[rt<<1]+sum[rt<<1|1])%=mod;
}
inline void pushdown(int rt,int len) {
int lf=len-(len>>1),rg=len>>1;
if (mul[rt]^1) {
(add[rt<<1]*=mul[rt])%=mod,(mul[rt<<1]*=mul[rt])%=mod,(sum[rt<<1]*=mul[rt])%=mod;
(add[rt<<1|1]*=mul[rt])%=mod,(mul[rt<<1|1]*=mul[rt])%=mod,(sum[rt<<1|1]*=mul[rt])%=mod;
mul[rt]=1;
}
if (add[rt]) {
(add[rt<<1]+=add[rt])%=mod,(sum[rt<<1]+=lf*add[rt]%mod)%=mod;
(add[rt<<1|1]+=add[rt])%=mod,(sum[rt<<1|1]+=rg*add[rt]%mod)%=mod;
add[rt]=0;
}
}
void modify(int rt,int l,int r,int L,int R,ll v,int opt) {
if (L<=l&&r<=R) {
if (opt&1) (add[rt]*=v)%=mod,(mul[rt]*=v)%=mod,(sum[rt]*=v)%=mod;
else (add[rt]+=v)%=mod,(sum[rt]+=(r-l+1)*v%mod)%=mod;
return ;
}
pushdown(rt,r-l+1);
int mid=l+r>>1;
if (L<=mid) modify(lson,L,R,v,opt);
if (mid<R) modify(rson,L,R,v,opt);
(sum[rt]=sum[rt<<1]+sum[rt<<1|1])%=mod;
}
ll query(int rt,int l,int r,int L,int R) {
if (L<=l&&r<=R) return sum[rt];
pushdown(rt,r-l+1);
int mid=l+r>>1;
ll ret=0;
if (L<=mid) (ret+=query(lson,L,R))%=mod;
if (mid<R) (ret+=query(rson,L,R))%=mod;
return ret;
}
int main() {
// freopen("P3373.in","r",stdin);
n=read(),m=read(),mod=read();
for (register int i=1;i<=n;++i) a[i]=read()%mod;
build(root);
for (register int t=0;t<m;++t) {
int opt=read();
switch (opt) {
case 1:{
int x=read(),y=read();
ll v=read();
modify(root,x,y,v,1);
break;
}
case 2:{
int x=read(),y=read();
ll v=read();
modify(root,x,y,v,2);
break;
}
case 3:{
int x=read(),y=read();
printf("%lld\n",query(root,x,y));
break;
}
}
}
return 0;
}
题意:nlogn实现区间加、区间乘、区间求和(都要求取模)
题解:
只要记住一句话:(x+a)*b=x*b+a*b,先乘后加,乘法标记管两个(sum,add),加法标记管一个(sum)。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define root 1,1,n
#define lson rt<<1,l,mid
#define rson rt<<1|1,mid+1,r
typedef long long ll;
const int N=1e5+4;
int n,m;
ll a
,mod,sum[N<<2],add[N<<2],mul[N<<2];
inline ll read() {
ll x=0,f=1;char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
return x*f;
}
void build(int rt,int l,int r) {
add[rt]=0,mul[rt]=1;
if (l==r) {sum[rt]=a[l];return ;}
int mid=l+r>>1;
build(lson);
build(rson);
(sum[rt]=sum[rt<<1]+sum[rt<<1|1])%=mod;
}
inline void pushdown(int rt,int len) {
int lf=len-(len>>1),rg=len>>1;
if (mul[rt]^1) {
(add[rt<<1]*=mul[rt])%=mod,(mul[rt<<1]*=mul[rt])%=mod,(sum[rt<<1]*=mul[rt])%=mod;
(add[rt<<1|1]*=mul[rt])%=mod,(mul[rt<<1|1]*=mul[rt])%=mod,(sum[rt<<1|1]*=mul[rt])%=mod;
mul[rt]=1;
}
if (add[rt]) {
(add[rt<<1]+=add[rt])%=mod,(sum[rt<<1]+=lf*add[rt]%mod)%=mod;
(add[rt<<1|1]+=add[rt])%=mod,(sum[rt<<1|1]+=rg*add[rt]%mod)%=mod;
add[rt]=0;
}
}
void modify(int rt,int l,int r,int L,int R,ll v,int opt) {
if (L<=l&&r<=R) {
if (opt&1) (add[rt]*=v)%=mod,(mul[rt]*=v)%=mod,(sum[rt]*=v)%=mod;
else (add[rt]+=v)%=mod,(sum[rt]+=(r-l+1)*v%mod)%=mod;
return ;
}
pushdown(rt,r-l+1);
int mid=l+r>>1;
if (L<=mid) modify(lson,L,R,v,opt);
if (mid<R) modify(rson,L,R,v,opt);
(sum[rt]=sum[rt<<1]+sum[rt<<1|1])%=mod;
}
ll query(int rt,int l,int r,int L,int R) {
if (L<=l&&r<=R) return sum[rt];
pushdown(rt,r-l+1);
int mid=l+r>>1;
ll ret=0;
if (L<=mid) (ret+=query(lson,L,R))%=mod;
if (mid<R) (ret+=query(rson,L,R))%=mod;
return ret;
}
int main() {
// freopen("P3373.in","r",stdin);
n=read(),m=read(),mod=read();
for (register int i=1;i<=n;++i) a[i]=read()%mod;
build(root);
for (register int t=0;t<m;++t) {
int opt=read();
switch (opt) {
case 1:{
int x=read(),y=read();
ll v=read();
modify(root,x,y,v,1);
break;
}
case 2:{
int x=read(),y=read();
ll v=read();
modify(root,x,y,v,2);
break;
}
case 3:{
int x=read(),y=read();
printf("%lld\n",query(root,x,y));
break;
}
}
}
return 0;
}
相关文章推荐
- Luogu 3373 又乘又加的线段树
- Luogu_3373 (双 tag 线段树模板)
- [luogu3373][模板]线段树 2(线段树)
- BZOJ3038 && BZOJ3211 上帝造题的七分钟2 && 花神游历各国 (线段树 + 开方标记)
- [ACM] poj 3468 A Simple Problem with Integers(线段树,成段更新,懒惰标记)
- 郁闷的出纳员 (splay的区间标记模板,删除区间,add标记,类似线段树)
- 线段树单标记下穿【模板】
- HDU 1698 Just a Hook 线段树 区间更新 惰性标记
- POJ 3468 线段树成段更新区间求和 + 延迟标记详解 A Simple Problem with Integers
- 浅析__线段树延迟标记
- Luogu 3369 我用线段树骗了一道平衡树题……
- HDU 1698-Just a Hook(线段树_区间更新+lazy标记)
- 线段树区间更新延迟标记
- 线段树 (区间修改 区间查询 延迟标记)
- 线段树_懒标记
- 线段树小引申(寻找最近的没被标记过的点)
- bzoj 1798 段乘段加线段树 双懒标记
- HDOJ 4027 Can you answer these queries?(线段树+区间标记)
- BZOJ_5249_Luogu_P4364_[2018多省省队联测]_IIIDX_九省联考2018_JLOI2018_线段树
- uva11992 Fast Matrix Operations 线段树打标记