LOJ#6283. 數列分塊入門 7

自為風月馬前卒發表於2018-02-19
記憶體限制:256 MiB時間限制:500 ms標準輸入輸出
題目型別:傳統評測方式:文字比較
上傳者: hzwer

題目描述

給出一個長為 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}1n100000,231​​others、ans≤231−1 \mathrm{ans} \leq 2^{31}-1ans231​​1。

顯示分類標籤

 

對於兩種標記,分別進行維護。

按照運算規律(這個可以自己推式子),我們先考慮乘法,再考慮加法

對於加的操作,相當於直接對於每個塊的加法標記加

對於乘的操作,我們需要對乘法標記和加法標記都乘上對應的值

零散塊不好記錄,直接暴力修改

 

 

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int MAXN=1e5+10,mod=10007;
inline char nc()
{
    static char buf[MAXN],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,MAXN,stdin)),p1==p2?EOF:*p1++;
}
inline int read()
{
    char c=nc();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=nc();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=nc();}
    return x*f;
}
int a[MAXN],belong[MAXN],L[MAXN],R[MAXN],block,N;
int mul[MAXN],add[MAXN];
void reset(int x)
{
    for(int i=L[x*block];i<=min(N,R[x*block]);i++)
        a[i]=(a[i]*mul[x]+add[x])%mod;
    mul[x]=1;add[x]=0;
}
void Add(int l,int r,int val)
{
    reset(belong[l]);
    for(int i=l;i<=min(r,R[l]);i++) a[i]+=val,a[i]%=mod;
    if(belong[l]!=belong[r])
    {
        reset(belong[r]);
        for(int i=L[r];i<=r;i++)
             a[i]+=val,a[i]%=mod;
    }
    for(int i=belong[l]+1;i<=belong[r]-1;i++)
        add[i]+=val,add[i]%=mod;
}
void Mul(int l,int r,int val)
{
    reset(belong[l]);
    for(int i=l;i<=min(r,R[l]);i++) a[i]*=val,a[i]%=mod;
    if(belong[l]!=belong[r])
    {
        reset(belong[r]);
        for(int i=L[r];i<=r;i++)
            a[i]*=val,a[i]%=mod;
    }
        
    for(int i=belong[l]+1;i<=belong[r]-1;i++)
        mul[i]*=val,add[i]*=val,mul[i]%=mod,add[i]%=mod;
}
int main()
{
    #ifdef WIN32
    freopen("a.in","r",stdin);
    #else
    #endif
    N=read();block=sqrt(N);
    fill(mul,mul+N+1,1);
    for(int i=1;i<=N;i++) a[i]=read()%mod,belong[i]=(i-1)/block+1;
    for(int i=1;i<=N;i++) L[i]=(belong[i]-1)*block+1,R[i]=belong[i]*block;
    for(int i=1;i<=N;i++)
    {
        int opt=read(),l=read(),r=read(),c=read()%mod;
        if(opt==0) 
            Add(l,r,c);
        else if(opt==1) 
            Mul(l,r,c);
        else if(opt==2) 
            printf("%d\n",(a[r]*mul[belong[r]]+add[belong[r]])%mod);
    }
    return 0;
}

 

 

 

 

相關文章