分块——入门7

这里写图片描述

如果只是加法操作的话就和入门1没啥区别了,现在最大的问题就是乘法咋搞??
首先我们要知道的是,乘法的优先级是高于加法,那么我们在做乘法的时候就要十分的小心了,所以我们在修改乘法标记的时候要把加法标记也作修改。
例如:
若当前的一个块乘以m1后加上a1,这时进行一个乘m2的操作,则原来的标记变成m1m2,a1m2
若当前的一个块乘以m1后加上a1,这时进行一个加a2的操作,则原来的标记变成m1,a1+a2

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
//By Bibi
/// .-~~~~~~~~~-._ _.-~~~~~~~~~-.
/// __.' ~. .~ `.__
/// .'// \./ \\`.
/// .'// | \\`.
/// .'// .-~"""""""~~~~-._ | _,-~~~~"""""""~-. \\`.
/// .'//.-" `-. | .-' "-.\\`.
/// .'//______.============-.. \ | / ..-============.______\\`.
/// .'______________________________\|/______________________________`.
#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define dep(i,a,b) for(int i=a;i>=b;--i)
using namespace std;
const int MAXN=100005;
const int mod=10007;
int read()
{
int sum=0,flag=1;
char c;
for(;c<'0'||c>'9';c=getchar())if(c=='-') flag=-1;
for(;c>='0'&&c<='9';c=getchar())sum=(sum<<1)+(sum<<3)+c-'0';
return sum*flag;
}
int n,block,num,l[MAXN],r[MAXN],L,R,c,op,a[MAXN],p[MAXN],be[MAXN],pp[MAXN];
void build()
{
block=sqrt(n);
num=n/block;
if(n%block) num++;
rep(i,1,n) be[i]=(i-1)/block+1;
rep(i,1,num)
{
l[i]=(i-1)*block+1;
r[i]=i*block;
}
r[num]=n;
}
void push(int ll,int rr)
{
rep(i,ll,rr) a[i]=(a[i]*pp[be[ll]]+p[be[ll]])%mod;
p[be[ll]]=0,pp[be[ll]]=1;
}

void add()
{
push(l[be[L]],r[be[L]]);
if(be[L]==be[R]) rep(i,L,R) a[i]+=c,a[i]%=mod;
else
{
rep(i,L,r[be[L]]) a[i]+=c,a[i]%=mod;
push(l[be[R]],r[be[R]]);
rep(i,l[be[R]],R) a[i]+=c,a[i]%=mod;
rep(i,be[L]+1,be[R]-1) p[i]+=c,p[i]%=mod;
}
}
void cheng()
{
push(l[be[L]],r[be[L]]);
if(be[L]==be[R]) rep(i,L,R) a[i]*=c,a[i]%=mod;
else
{
rep(i,L,r[be[L]]) a[i]*=c,a[i]%=mod;
push(l[be[R]],r[be[R]]);
rep(i,l[be[R]],R) a[i]*=c,a[i]%=mod;
rep(i,be[L]+1,be[R]-1) pp[i]*=c,p[i]*=c,pp[i]%=mod,p[i]%=mod;
}
}
void init()
{
n=read();
rep(i,1,n) a[i]=read(),pp[i]=1;
build();
rep(i,1,n)
{
op=read();L=read();R=read();c=read();
if(!op) add();
else if(op==1) cheng();
else printf("%d\n",(a[R]*pp[be[R]]+p[be[R]])%mod);
}
}
int main()
{
init();
return 0;
}