分块+二分专题

洛谷P2801 教主的魔法

区间修改\区间查询大于C的个数

#include<bits/stdc++.h>
using namespace std;
const int MAX=2*1e6;
const int eps=1e-8;
int a[MAX],aa[MAX],b[MAX],c[MAX],cc[MAX][2],block,cnt;
inline void read(int &x)
{
    x=0;char s=getchar();
    while(s<'0'||s>'9'){s=getchar();}
    while(s<='9'&&s>='0'){x=x*10+s-'0';s=getchar();}
}
void sort(int x)
{
    int L=cc[x][0],R=cc[x][1];
    for(int k=L;k<=R;k++) aa[k]=a[k];
    sort(aa+L,aa+R+1);
}
void add(int l,int r,int w)
{
    bool bl=l==cc[c[l]][0],br=r==cc[c[r]][1];
    if(c[l]==c[r])
    {
        for(int k=l;k<=r;k++) a[k]+=w; sort(c[l]);
    }
    else
    {
        if(!bl){for(int k=l;k<=cc[c[l]][1];k++) a[k]+=w; sort(c[l]);}
        if(!br){for(int k=cc[c[r]][0];k<=r;k++) a[k]+=w; sort(c[r]);}
        for(int k=c[l]+1-bl;k<c[r]+br;k++) b[k]+=w;
    }
}
int query(int l,int r,int w)
{
    int ans=0;
    bool bl=l==cc[c[l]][0],br=r==cc[c[r]][1];
    if(c[l]==c[r])
    {
        for(int k=l,v=w-b[c[l]];k<=r;k++) ans+=a[k]>=v;
    }
    else
    {
        if(!bl)for(int k=l,v=w-b[c[l]];k<=cc[c[l]][1];k++) ans+=a[k]>=v;
        if(!br)for(int k=cc[c[r]][0],v=w-b[c[r]];k<=r;k++) ans+=a[k]>=v;
        for(int k=c[l]+1-bl;k<c[r]+br;k++)
        {
            int L=cc[k][0],R=cc[k][1],v=w-b[k],mid;
            while(L<=R)
            {
                mid=(L+R)>>1;
                aa[mid]>=v?R=mid-1:L=mid+1;
            }
            ans+=cc[k][1]-L+1;
        }
    }
    return ans;
}
int main()
{
    int n,m,x,y,z;
    char s[5];
    read(n);read(m);
    memset(a,0,sizeof(a));
    block=(int)(sqrt(n)+eps),cnt=1;
    for(int k=1;k<=n;k++) c[k]=k%block==0?cnt++:cnt;
    for(int k=1;k<=n;k++) read(a[k]);
    for(int k=1;k<=cnt;k++) {cc[k][0]=k*block-block+1;cc[k][1]=k*block;}cc[cnt][1]=n;
    for(int k=1;k<=cnt;k++) sort(k);
    for(int k=0;k<m;k++)
    {
        scanf("%s",s);
        read(x);read(y);read(z);
        if(s[0]=='M') add(x,y,z);
        else printf("%d\n",query(x,y,z));
    }
    return 0;
}

qustoj 二与安生

单点修改\区间查询次大值

#include<bits/stdc++.h>
using namespace std;
const int MAX=1e6+5;
const double eps=1e-6;
int n,m,block,cnt,a[MAX],c[MAX],cc[MAX][2],b[MAX][2];
void update(int x)
{
    int L=cc[x][0],R=cc[x][1];
    b[x][0]=b[x][1]=b[x][2]-1;
    for(int i=L;i<=R;++i)
    {
        if(a[i]>a[b[x][0]])
        {
            b[x][1]=b[x][0];
            b[x][0]=i;
        }
        else if(a[i]>a[b[x][1]])
        {
            b[x][1]=i;
        }
    }
    if(b[x][0]>b[x][1]) swap(b[x][0],b[x][1]);
}
int query(int l,int r)
{
    if(r-l<1) return -1;
    int m1,m2,m3,m4;
    m1=m2=m3=m4=-1;
    bool bl=l==cc[c[l]][0],br=r==cc[c[r]][1];
    if(c[l]==c[r])
    {
        for(int k=l;k<=r;k++)
        if(a[k]>m1) {m2=m1;m1=a[k];m4=m3;m3=k;}
        else if(a[k]>m2) {m2=a[k];m4=k;}
    }
    else
    {
        if(!bl)for(int k=l;k<=cc[c[l]][1];k++)
        {
            if(a[k]>m1) {m2=m1;m1=a[k];m4=m3;m3=k;}
            else if(a[k]>m2) {m2=a[k];m4=k;}
        }
        for(int k=c[l]+1-bl;k<c[r]+br;k++)
        {
            if(a[b[k][0]]>m1) {m2=m1;m1=a[b[k][0]];m4=m3;m3=b[k][0];}
            else if(a[b[k][0]]>m2) {m2=a[b[k][0]];m4=b[k][0];}
            if(a[b[k][1]]>m1) {m2=m1;m1=a[b[k][1]];m4=m3;m3=b[k][1];}
            else if(a[b[k][1]]>m2) {m2=a[b[k][1]];m4=b[k][1];}
        }
        if(!br)for(int k=cc[c[r]][0];k<=r;k++)
        {
            if(a[k]>m1) {m2=m1;m1=a[k];m4=m3;m3=k;}
            else if(a[k]>m2) {m2=a[k];m4=k;}
        }
    }
    return m4;
}
int main()
{
//    freopen("in.txt","r",stdin);
//    freopen("out.txt","w",stdout);
    scanf("%d%d",&n,&m);
    block=(int)(sqrt(n)+eps);
    cnt=1;
    for(int i=1;i<=n;++i) c[i]=i%block==0?cnt++:cnt;
    for(int i=1;i<cnt;++i) {cc[i][0]=i*block-block+1;cc[i][1]=i*block;}
    cc[cnt-1][1]=n;
    for(int i=1;i<=n;++i) scanf("%d",&a[i]);
    for(int i=1;i<cnt;++i) update(i);
    int f,x,y;
    for(int i=0;i<m;++i)
    {
        scanf("%d%d%d",&f,&x,&y);
        if(f==1) {a[x]=y;update(c[x]);}
        else printf("%d\n",query(x,y));
    }
    return 0;
}