LOJ #6283. 數列分塊入門 7-分塊(區間乘法、區間加法、單點查詢)

#6283. 數列分塊入門 7

內存限制：256 MiB時間限制：500 ms標准輸入輸出

題目類型：傳統評測方式：文本比較

上傳者： hzwer

提交 提交記錄 統計 測試數據 討論

2

 

題目描述

給出一個長為 nn 的數列，以及 nn 個操作，操作涉及區間乘法，區間加法，單點詢問。

輸入格式

第一行輸入一個數字 nn。

第二行輸入 nn 個數字，第 ii 個數字為 a_iai​，以空格隔開。

接下來輸入 nn 行詢問，每行輸入四個數字 \mathrm{opt}opt、ll、rr、cc，以空格隔開。

若 \mathrm{opt} = 0opt=0，表示將位於 [l, r][l,r] 的之間的數字都加 cc。

若 \mathrm{opt} = 1opt=1，表示將位於 [l, r][l,r] 的之間的數字都乘 cc。

若 \mathrm{opt} = 2opt=2，表示詢問 a_rar​ 的值 \mathop{\mathrm{mod}} 10007mod10007（ll 和 cc 忽略）。

輸出格式

對於每次詢問，輸出一行一個數字表示答案。

樣例

樣例輸入

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% 的數據，1 \leq n \leq 100000, -2^{31} \leq \mathrm{others}1≤n≤100000,−231≤others、\mathrm{ans} \leq 2^{31}-1ans≤231−1。

 

這道題要開比1e5大，要不然過不去，真實測試。。。

 

代碼:

//#6283. 數列分塊入門 7-區間乘法，區間加法，單點查詢-要開1e6，mdzz
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e6+10;

int n,m;
int a[maxn],add[maxn],mul[maxn],pos[maxn];
const int mod=10007;

void update_add(int l,int r,int c)
{
    if(pos[l]==pos[r]){
        for(int i=(pos[l]-1)*m+1;i<=pos[l]*m;i++){
            a[i]=(a[i]*mul[pos[l]]%mod+add[pos[l]]+mod)%mod;
        }
        for(int i=l;i<=r;i++)
            a[i]=(a[i]+c+mod)%mod;
        mul[pos[l]]=1;add[pos[l]]=0;
    }
    else{
        for(int i=(pos[l]-1)*m+1;i<=pos[l]*m;i++){
            a[i]=(a[i]*mul[pos[l]]%mod+add[pos[l]]+mod)%mod;
        }
        for(int i=l;i<=pos[l]*m;i++)
            a[i]=(a[i]+c+mod)%mod;
        mul[pos[l]]=1;add[pos[l]]=0;
        for(int i=pos[l]+1;i<pos[r];i++){
            add[i]=(add[i]+c+mod)%mod;
        }
        for(int i=(pos[r]-1)*m+1;i<=min(pos[r]*m,n);i++){
            a[i]=(a[i]*mul[pos[r]]%mod+add[pos[r]]+mod)%mod;
        }
        for(int i=(pos[r]-1)*m+1;i<=r;i++)
            a[i]=(a[i]+c+mod)%mod;
        mul[pos[r]]=1;add[pos[r]]=0;
    }
}

void update_mul(int l,int r,int c)
{
    if(pos[l]==pos[r]){
        for(int i=(pos[l]-1)*m+1;i<=pos[l]*m;i++){
            a[i]=(a[i]*mul[pos[l]]%mod+add[pos[l]]+mod)%mod;
        }
        for(int i=l;i<=r;i++)
            a[i]=(a[i]*c+mod)%mod;
        mul[pos[l]]=1;add[pos[l]]=0;
    }
    else{
        for(int i=(pos[l]-1)*m+1;i<=pos[l]*m;i++){
            a[i]=(a[i]*mul[pos[l]]%mod+add[pos[l]]+mod)%mod;
        }
        for(int i=l;i<=pos[l]*m;i++)
            a[i]=(a[i]*c+mod)%mod;
        mul[pos[l]]=1;add[pos[l]]=0;
        for(int i=pos[l]+1;i<pos[r];i++){
            add[i]=(add[i]*c+mod)%mod;
            mul[i]=(mul[i]*c+mod)%mod;
        }
        for(int i=(pos[r]-1)*m+1;i<=min(pos[r]*m,n);i++){
            a[i]=(a[i]*mul[pos[r]]%mod+add[pos[r]]+mod)%mod;
        }
        for(int i=(pos[r]-1)*m+1;i<=r;i++)
            a[i]=(a[i]*c+mod)%mod;
        mul[pos[r]]=1;add[pos[r]]=0;
    }
}

int main()
{
    scanf("%d",&n);
    m=sqrt(n);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        a[i]%=mod;
        pos[i]=(i-1)/m+1;
        mul[pos[i]]=1;
        add[pos[i]]=0;
    }
    for(int i=1;i<=n;i++){
        int op,l,r,c;
        scanf("%d%d%d%d",&op,&l,&r,&c);
        c%=mod;
        if(op==0){
            update_add(l,r,c);
        }
        else if(op==1){
            update_mul(l,r,c);
        }
        else if(op==2){
            printf("%d\n",(a[r]*mul[pos[r]]+add[pos[r]])%mod);
        }
    }
}


/*
10
1 3 4 2 5 7 11 3 5 1
0 1 5 1
1 1 7 2
2 3 9 1
0 4 8 1
1 1 5 2
1 3 5 2
2 5 7 1
1 3 5 2
2 2 3 2
2 3 4 5

5
23
80
56
*/