数列分块入门1~3

数列分块入门1：

题目：

        给出一个长为nn 的数列，以及 nn 个操作，操作涉及区间加法，单点查值。

Input

第一行输入一个数字nn。 
第二行输入nn个数字，第ii个数字为aiai，以空格隔开。 
接下来输入nn行询问，每行输入四个数字opt,l,r,copt,l,r,c，以空格隔开。 
若opt=0opt=0，表示将位于[l,r][l,r]的之间的数字都加cc。 
若opt=1opt=1，表示询问arar 的值（ll和cc忽略）。 

Output

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

Samples

Input 

4
1 2 2 3
0 1 3 1
1 0 1 0
0 1 2 2
1 0 2 0

Output

2
5

解题思路:

        完整块数量不超过sqrt(n)，不完整块是长度小于block的块，对于完整块加值可以标记记录，对于不完整块进行暴力修改元素的值，这样复杂度最小为O(n3/2)。

代码：

#include<bits/stdc++.h> 
#define ll long long
using namespace std;
int n,a[1001100],blo[1001100],atag[1001100];
void update(int l,int r,int c){
    for(int i=l;i<=min(blo[l]*blo[0],r);i++)
        a[i]+=c;
    if(blo[l]!=blo[r])
        for(int i=(blo[r]-1)*blo[0]+1;i<=r;i++)
            a[i]+=c;
    for(int i=blo[l]+1;i<=blo[r]-1;i++)
        atag[i]+=c;
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) cin>>a[i];
    blo[0]=sqrt(n);
    for(int i=1;i<=n;i++) blo[i]=(i-1)/blo[0]+1;
    for(int i=1;i<=n;i++){
        int opt,l,r,c;
        cin>>opt>>l>>r>>c;
        if(!opt) update(l,r,c);
        else cout<<a[r]+atag[blo[r]]<<endl;
    }
	return 0;
}

Hint

n≤50000n≤50000,−231≤others,ans≤231−1−231≤others,ans≤231−1

数列分块入门2

题目：

        给出一个长为 nn 的数列，以及 nn 个操作，操作涉及区间加法，询问区间内小于某个值 xx 的元素个数。

Input

第一行输入一个数字 nn。 
第二行输入 nn 个数字，第 ii 个数字为 aiai，以空格隔开。 
接下来输入 nn 行询问，每行输入四个数字opt,l,r,copt,l,r,c，以空格隔开。 
若opt=0opt=0，表示将位于[l,r][l,r]的之间的数字都加cc。 
若opt=1opt=1，表示询问中[l,r][l,r]小于c2c2的数字的个数。 

Output

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

Samples

Input 

4
1 2 2 3
0 1 3 1
1 1 3 2
1 1 4 1
1 2 3 2

Output

3
0
2

解题思路：

        对于不完整块用暴力进行查找小于X的元素，对于完整块可以使用二分查找，对于整块相对顺序不会改变，不需要重新排序，对于不完整快需要暴力修改之后再重新排序。

代码：

#include<bits/stdc++.h> 
#define ll long long
using namespace std;
int rd(){
    int x=0,f=1;
	char ch=getchar();
    while(ch<'0'||ch>'9'){
		if(ch=='-')f=-1;ch=getchar();
	}
    while(ch>='0'&&ch<='9'){
		x=x*10+ch-'0';ch=getchar();
	}
    return x*f;
}
int n,blo,v[100001],tag[100001],bl[100001];
vector <int> vec[1001];
void reset(int x) {
    vec[x].clear();
    for(int i=(x-1)*blo+1;i<=min(x*blo,n);i++) 
        vec[x].push_back(v[i]);
    sort(vec[x].begin(),vec[x].end());
}
void add(int a,int b,int c) {
    for(int i=a;i<=min(bl[a]*blo,b);i++) v[i]+=c;
    reset(bl[a]);
    if(bl[a]!=bl[b]){
        for(int i=(bl[b]-1)*blo+1;i<=b;i++) v[i]+=c;
        reset(bl[b]);
    }
    for(int i=bl[a]+1;i<=bl[b]-1;i++) tag[i]+=c;
}

int query(int l,int r,int c) {
    int ans=0;
    for(int i=l;i<=min(bl[l]*blo,r);i++) 
        if(v[i]+tag[bl[l]]<c) ans++;
    if(bl[l]!=bl[r])
        for(int i=(bl[r]-1)*blo+1;i<=r;i++)
            if(v[i]+tag[bl[r]]<c) ans++;
    for(int i=bl[l]+1;i<=bl[r]-1;i++){
        int x=c-tag[i];
        ans+=lower_bound(vec[i].begin(),vec[i].end(),x)-vec[i].begin();
    }
    return ans;
} 
int main() {
    n=rd(),blo=sqrt(n);
    for(int i=1;i<=n;i++) v[i]=rd();
    for(int i=1;i<=n;i++){
        bl[i]=(i-1)/blo+1;
        vec[bl[i]].push_back(v[i]);
    }
    for(int i=1;i<=bl[n];i++) 
        sort(vec[i].begin(),vec[i].end());
    for(int i=1;i<=n;i++) {
        int opt=rd(),l=rd(),r=rd(),c=rd();
        if(!opt) add(l,r,c);
        else cout<<query(l,r,c*c)<<endl;
    }
    return 0;
}

Hint

n≤50000n≤50000, −231≤others,ans≤231−1

数列分块入门3

题目：

        给出一个长为 nn  的数列，以及 nn 个操作，操作涉及区间加法，询问区间内小于某个值 xx 的前驱（比其小的最大元素）。

Input

第一行输入一个数字 nn 。 

第二行输入 nn 个数字，第 ii 个数字为 aiai，以空格隔开。 
接下来输入 nn 行询问，每行输入四个数字opt,l,r,copt,l,r,c，以空格隔开。 
若opt=0opt=0，表示将位于[l,r][l,r]的之间的数字都加cc。 
若opt=1opt=1，表示询问[l,r][l,r]中cc的前驱的值（不存在则输出-1）。 

Output

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

Samples

Input 

4
1 2 2 3
0 1 3 1
1 1 4 4
0 1 2 2
1 1 2 4

Output

3
-1

解题思路：

        与数列分块入门二相似，用到排序加二分查找即可。

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int a[100010],b[100010],pos[100010],ans,flag,n,block;
vector<int>v[1010];
void reset(int x)
{
    v[pos[x]].clear();
    for(int i=(pos[x]-1)*block+1;i<=min(pos[x]*block,n);i++)
        v[pos[x]].push_back(a[i]);
    sort(v[pos[x]].begin(),v[pos[x]].end());
} 
void update(int l,int r,int c) {
    for(int i=l;i<=min(pos[l]*block,r);i++) a[i]+=c;
    reset(l);
    if(pos[l]!=pos[r]){
        for(int i=(pos[r]-1)*block+1;i<=r;i++) a[i]+=c;
        reset(r);
    }
    for(int i=pos[l]+1;i<=pos[r]-1;i++)
        b[i]+=c;
}
void query(int l,int r,int x)
{
    flag=0;
    ans=-10001010;
    for(int i=l;i<=min(pos[l]*block,r);i++)
        if(a[i]+b[pos[l]]<x){
            ans=max(a[i]+b[pos[l]],ans);
            flag=1;    
        } 
    if(pos[l]!=pos[r])
        for(int i=(pos[r]-1)*block+1;i<=r;i++)
            if(a[i]+b[pos[r]]<x){
                ans=max(ans,a[i]+b[pos[r]]);
                flag=1;
            }
    for(int i=pos[l]+1;i<=pos[r]-1;i++){
        int c=x-b[i];
        int k=lower_bound(v[i].begin(),v[i].end(),c)-v[i].begin();
        if(k!=0&&v[i][k-1]+b[i]<x) {
            flag=1;
            ans=max(ans,v[i][k-1]+b[i]);
        }
    }
    if(flag) cout<<ans<<endl;
    else cout<<"-1"<<endl;
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) cin>>a[i];
    block=sqrt(n); 
    for(int i=1;i<=n;i++)
        pos[i]=(i-1)/block+1,v[pos[i]].push_back(a[i]);
    for(int i=1;i<=pos[n];i++) sort(v[i].begin(),v[i].end());
    int opt,l,r,c;
    for(int i=0;i<n;i++){
        cin>>opt>>l>>r>>c;
        if(opt==0) update(l,r,c);
        else query(l,r,c);        
    } 
    return 0;
}

Hint

n≤100000n≤100000, −231≤others,ans≤231−1