loj #6283. 數列分塊入門 7

#6283. 數列分塊入門 7

題目描述

給出一個長爲 nnn 的數列，以及 nnn 個操做，操做涉及區間乘法，區間加法，單點詢問。

輸入格式

第一行輸入一個數字 nnn。

第二行輸入 nnn 個數字，第 i 個數字爲 aia_ia​i​​，以空格隔開。

接下來輸入 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_ra​r​​ 的值 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,−2​31​​≤others、ans≤231−1 \mathrm{ans} \leq 2^{31}-1ans≤2​31​​−1$。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define maxn 100010
#define mod 10007
using namespace std;
int n,block,a[maxn],pos[maxn],add[maxn],mul[maxn];
int findl(int id){return (id-1)*block+1;}
int findr(int id){return min(id*block,n);}
void reset(int id){
    int l=findl(id),r=findr(id);
    for(int i=l;i<=r;i++){
        a[i]=1LL*a[i]*mul[id]%mod;
        a[i]=(a[i]+add[id])%mod;
    }
    mul[id]=1;
    add[id]=0;
}
void modify1(int l,int r,int c){
    reset(pos[l]);
    for(int i=l;i<=min(findr(pos[l]),r);i++)(a[i]+=c)%=mod;
    if(pos[l]==pos[r])return;
    reset(pos[r]);
    for(int i=findl(pos[r]);i<=r;i++)(a[i]+=c)%=mod;
    for(int i=pos[l]+1;i<=pos[r]-1;i++)
        (add[i]+=c)%=mod;
}
void modify2(int l,int r,int c){
    reset(pos[l]);
    for(int i=l;i<=min(findr(pos[l]),r);i++)a[i]=1LL*a[i]*c%mod;
    if(pos[l]==pos[r])return;
    reset(pos[r]);
    for(int i=findl(pos[r]);i<=r;i++)a[i]=1LL*a[i]*c%mod;
    for(int i=pos[l]+1;i<=pos[r]-1;i++){
        add[i]=1LL*add[i]*c%mod;
        mul[i]=1LL*mul[i]*c%mod;
    }
}
int query(int x){
    int id=0,suf=0;
    while(1){
        suf+=block;
        if(suf>=x){
            if(suf>n)id=pos[n];
            else id=pos[suf];
            break;
        }
    }
    return (1LL*a[x]*mul[id]%mod+add[id])%mod;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    block=sqrt(n);
    for(int i=1;i<=n;i++)pos[i]=(i-1)/block+1,mul[i]=1;
    int op,l,r,c;
    for(int i=1;i<=n;i++){
        scanf("%d%d%d%d",&op,&l,&r,&c);
        if(op==0)modify1(l,r,c);
        else if(op==1)modify2(l,r,c);
        else printf("%d\n",query(r));
    }
    return 0;
}