首页 > > 程序设计 > C/C++ >

LOJ#6280. 数列分块入门 4

2018-06-17 21:14:05来源：未知阅读 ()

内存限制：256 MiB时间限制：500 ms标准输入输出

题目类型：传统评测方式：文本比较

上传者： hzwer

提交提交记录统计讨论测试数据

给出一个长为

第一行输入一个数字

第二行输入 $aia_ia?i??，以空格隔开。$

接下来输入 $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] 的所有数字的和 mod(c+1)mod (c+1)mod(c+1)。$

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

样例输入

样例输出

1
4

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

$分块$

$维护块的值$

$其余的按照套路来$

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define int long long 
using namespace std;
const int MAXN=2*1e6+10;
const int INF=1e8+10;
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],sum[MAXN],tag[MAXN],block;
void IntervalAdd(int l,int r,int val)
{
    for(int i=l;i<=min(R[l],r);i++) a[i]+=val,sum[belong[i]]+=val;
    if(belong[l]!=belong[r])
        for(int i=L[r];i<=r;i++)
            a[i]+=val,sum[belong[i]]+=val;
    for(int i=belong[l]+1;i<=belong[r]-1;i++)
        tag[i]+=val;
}
int Query(int l,int r,int mod)
{
    int ans=0;
    for(int i=l;i<=min(R[l],r);i++) 
        ans+=a[i]+tag[belong[i]],ans=ans%mod;
    
    if(belong[l]!=belong[r])
        for(int i=L[r];i<=r;i++)
            ans+=a[i]+tag[belong[i]];
    for(int i=belong[l]+1;i<=belong[r]-1;i++)
        ans+=sum[i]+tag[i]*block,ans=ans%mod;
    return ans%mod;
}
main()
{
    #ifdef WIN32
    freopen("a.in","r",stdin);
   // freopen("b.out","w",stdout);
    #else
    #endif
    int N=read();block=sqrt(N);
    for(int i=1;i<=N;i++) belong[i]=(i-1)/block+1,L[i]=(belong[i]-1)*block+1,R[i]=belong[i]*block;
    for(int i=1;i<=N;i++) a[i]=read();
    for(int i=1;i<=N;i++) sum[belong[i]]+=a[i];
    for(int i=1;i<=N;i++)
    {
        int opt=read(),l=read(),r=read(),c=read();
        if(opt==0) IntervalAdd(l,r,c);
        else        
            printf("%lld\n",Query(l,r,c+1));
     } 
    return 0;
}