数列分块入门 1-9 未完待续

最近又去研究暴力美学。。

深深的感到暴力之美。

 

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

思路：把数列分成根号n块 对于一个区间l r 先暴力处理掉连边不包含整个块的，然后对包含整个块的直接打上标记

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define sca(x) scanf("%d",&x)
#define per(i,j,k) for(int i=j;i>=k;i--)
#define inf 0x3f3f3f3f
#define LL long long
#define N 550005
#define inf 0x3f3f3f3f

int len;
int a[N],b[N],tag[N];
int read()
{
    char c;
    int x=0,f=1;
    c=getchar();
    while(c<'0'||c>'9')
    {
        if(c=='-')f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9')
    {
        x=x*10+c-'0';
        c=getchar();
    }
    return x*f;
}
void add(int l,int r,int k)
{
    for(int i=l;i<=min(b[l]*len,r);i++)a[i]+=k;
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1;i<=r;i++)a[i]+=k;
    }
    for(int i=b[l]+1;i<=b[r]-1;i++)tag[i]+=k;
}

int main()
{
    int n;
    int op,l,r,k;
    n=read();
    len=sqrt(n);
    rep(i,1,n)a[i]=read();
    rep(i,1,n)b[i]=(i-1)/len+1;
    rep(i,1,n)
    {
        op=read(),l=read(),r=read(),k=read();
        if(op==0)
        {
            add(l,r,k);
        }
        else
        {
            printf("%d\n",a[r]+tag[b[r]]);
        }
    }
    return 0;
}

题目2：给出一个长为 n的数列，以及 n 个操作，操作涉及区间加法，询问区间内小于某个值 x 的元素个数

思路：依旧是分块对于每个块内的数先排好序，不在整个块内的依然是暴力，对于在块内的就二分查询

注意不完整的块更新完后要重新排序。

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define sca(x) scanf("%d",&x)
#define per(i,j,k) for(int i=j;i>=k;i--)
#define inf 0x3f3f3f3f
#define LL long long
#define N 550005
#define pb(x) push_back(x)
#define inf 0x3f3f3f3f

int len;
int a[N],b[N],tag[N];
vector<int>v[25000];
int read()
{
    char c;int x=0,f=1;
    c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

void modify(int x)
{
    v[x].clear();
    for(int i=(x-1)*len+1;i<=x*len;i++)v[x].pb(a[i]);
    sort(v[x].begin(),v[x].end());
}

void add(int l,int r,int k)
{
    for(int i=l;i<=min(b[l]*len,r);i++)a[i]+=k;
    modify(b[l]);//l所在的整个块重新排序
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1;i<=r;i++)a[i]+=k;
        modify(b[r]);//r所在的块重新排序
    }
    for(int i=b[l]+1;i<=b[r]-1;i++)tag[i]+=k;
}

int query(int l,int r,int k)
{
    int ans=0;
    for(int i=l;i<=min(b[l]*len,r);i++)
    {
        if(a[i]+tag[b[l]]<k)ans++;
    }
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1;i<=r;i++)
        {
            if(a[i]+tag[b[r]]<k)ans++;
        }
    }
    for(int i=b[l]+1;i<=b[r]-1;i++)
    {
        int tmp=k-tag[i];
        ans+=lower_bound(v[i].begin(),v[i].end(),tmp)-v[i].begin();
    }
    return ans;
}
int main()
{
    int n;
    n=read();
    int op,l,r,c;
    len=sqrt(n);
    rep(i,1,n)a[i]=read();
    rep(i,1,n)b[i]=(i-1)/len+1;
    rep(i,1,n)v[b[i]].pb(a[i]);
    rep(i,1,b[n])sort(v[i].begin(),v[i].end());//先对每个块排序
    rep(i,1,n)
    {
        op=read(),l=read(),r=read(),c=read();
        if(op==0)add(l,r,c);
        else printf("%d\n",query(l,r,c*c));
    }
}

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

思路：与上一题基本上相似，这个要熟悉lower bound 的用法，测试数据里边应该是没有负数的  所以ans的值直接设位-1.

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define sca(x) scanf("%d",&x)
#define per(i,j,k) for(int i=j;i>=k;i--)
#define inf 0x3f3f3f3f
#define LL long long
#define N 550005
#define pb(x) push_back(x)
#define inf 0x3f3f3f3f

int len;
int a[N],b[N],tag[N];
vector<int>v[25000];
int read()
{
    char c;int x=0,f=1;
    c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

void modify(int x)
{
    v[x].clear();
    for(int i=(x-1)*len+1;i<=x*len;i++)v[x].pb(a[i]);
    sort(v[x].begin(),v[x].end());
}

void add(int l,int r,int k)
{
    for(int i=l;i<=min(b[l]*len,r);i++)a[i]+=k;
    modify(b[l]);
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1;i<=r;i++)a[i]+=k;
        modify(b[r]);
    }
    for(int i=b[l]+1;i<=b[r]-1;i++)tag[i]+=k;
}

int query(int l,int r,int k)
{
    int ans=-1;
    for(int i=l;i<=min(b[l]*len,r);i++)
    {
        if(a[i]+tag[b[l]]<k)ans=max(ans,a[i]+tag[b[l]]);
    }
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1;i<=r;i++)
        {
            if(a[i]+tag[b[r]]<k)ans=max(ans,a[i]+tag[b[r]]);
        }
    }
    for(int i=b[l]+1;i<=b[r]-1;i++)
    {
        int tmp=k-tag[i];
        int pos=lower_bound(v[i].begin(),v[i].end(),tmp)-v[i].begin();
        if(pos>=1)ans=max(ans,v[i][pos-1]+tag[i]);
    }
    return ans;
}
int main()
{
    int n;
    n=read();
    int op,l,r,c;
    len=sqrt(n);
    rep(i,1,n)a[i]=read();
    rep(i,1,n)b[i]=(i-1)/len+1;
    rep(i,1,n)v[b[i]].pb(a[i]);
    rep(i,1,b[n])sort(v[i].begin(),v[i].end());
    rep(i,1,n)
    {
        op=read(),l=read(),r=read(),c=read();
        if(op==0)add(l,r,c);
        else printf("%d\n",query(l,r,c));
    }
}