Loj 6283. 数列分块入门 7

题目描述

给出一个长为 nnn 的数列，以及 nnn 个操作，操作涉及区间乘法，区间加法，单点询问。

输入格式

第一行输入一个数字 nnn。

第二行输入 nnn 个数字，第 i 个数字为 aia_iai​，以空格隔开。

接下来输入 nnn 行询问，每行输入四个数字 opt\mathrm{opt}opt、lll、rrr、ccc，以空格隔开。

若 opt=0\mathrm{opt} = 0opt=0，表示将位于 [l,r][l, r][l,r] 的之间的数字都加 ccc。

若 opt=1\mathrm{opt} = 1opt=1，表示将位于 [l,r][l, r][l,r] 的之间的数字都乘 ccc。

若 opt=2\mathrm{opt} = 2opt=2，表示询问 ara_rar​ 的值 mod 10007mod \ 10007mod 10007（lll 和 ccc 忽略）。

输出格式

对于每次询问，输出一行一个数字表示答案。

样例

样例输入

7
1 2 2 3 9 3 2
0 1 3 1
2 1 3 1
1 1 4 4
0 1 7 2
1 2 6 4
1 1 6 5
2 2 6 4

样例输出

3
100

数据范围与提示

对于 100% 100\%100% 的数据，1≤n≤100000,−231≤others 1 \leq n \leq 100000, -2^{31} \leq \mathrm{others}1≤n≤100000,−231≤others、ans≤231−1 \mathrm{ans} \leq 2^{31}-1ans≤231−1。

思路： 把加法和乘法都看成 x*a+b，乘法的时候b是0，加法的时候a是1

#include<bits/stdc++.h>
using namespace std;

#define F(i,a,b) for(int i=a;i<=b;i++) 
#define D(i,a,b) for(int i=a;i>=b;i--)
#define ms(i,a)  memset(a,i,sizeof(a)) 
#define LL       long long 
#define st(x)    ((x-1)*B+1)
#define ed(x)    min(x*B,n)
#define bl(x)    ((x-1)/B+1) 

int inline read(){
    int x=0,w=0; char c=getchar(); 
    while (c<'0' || c>'9') w+=c=='-',c=getchar(); 
    while (c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar(); 
    return w? -x: x;  
}

int const maxn=100003; 
int const mod=10007;  

int n,B;
LL add[350],mul[350],a[maxn];  

void update(int l,int r,int M,int A){
    int x=bl(l); 
    int y=bl(r); 
    if(x==y){
        F(i,st(x),ed(x)) a[i]=(a[i]*mul[x]+add[x]) % mod;  
        F(i,l,r) a[i]=(a[i]*M+A) % mod;   
        mul[x]=1; add[x]=0;  
    }else {
        F(i,st(x),ed(x)) a[i]=(a[i]*mul[x]+add[x]) % mod;  
        F(i,st(y),ed(y)) a[i]=(a[i]*mul[y]+add[y]) % mod; 
        F(i,l,ed(x)) a[i]=(a[i]*M+A) % mod; 
        F(i,st(y),r) a[i]=(a[i]*M+A) % mod;    
        mul[x]=mul[y]=1; add[x]=add[y]=0;  
        F(i,x+1,y-1) mul[i]=(mul[i]*M)% mod,add[i]=(add[i]*M+A) % mod;  
    }
}

int query(int k){
    int x=bl(k); 
    F(i,st(x),ed(x)) a[i]=(a[i]*mul[x]+add[x]) % mod;
    mul[x]=1; add[x]=0;  
    return a[k];  
}

int main(){
    n=read(); 
    B=(int) sqrt(n); 
    F(i,1,n) a[i]=read()% mod; 
    F(i,1,bl(n)) mul[i]=1; 
    F(i,1,n){
        int x=read(); 
        int l=read(); 
        int r=read(); 
        int c=read(); 
        if(x==0) update(l,r,1,c); 
        else if(x==1) update(l,r,c,0); 
        else printf("%d\n",query(r));  
    }
    return 0; 
}