CodeForces 551E GukiZ and GukiZiana 【分块】

题目：点击打开链接

题意：对于给定数列，操作一是把[l,r]区间内的数加x，操作二是求出一个值，该值为数x出现的最右位置减最左位置，若不存在x输出为-1。

思路：对于这种迷之询问，用分块大法。将n个数分成sqrt(n)块，初始化将a数组赋值给d数组，将d数组分块从小到大排序，新建一个数组p储存每一块所加的和，更新时，如果l,r都在同一块，暴力更新a并再排序d，若不在同一块，对于l所在块，暴力找出所有本来位置在a数组l和l之后的数加上x再给d排序，对于r所在块，暴力找出所有本来位置在a数组r和r之前的数加上x再给d排序，对于中间的块，则p数组加上x。二分查询时，将x减去p数组的对应的值，对于d从左往右每一块进行二分查找该数的最左端位置，找到后再继续二分查找每一块的该数的最右端位置。若找不到最左端的数则输出-1，否则输出两位置之差。

代码：

#include <stdio.h>
#include <cmath>
#include <algorithm>
#include <queue>
#include <set>
#include <vector>
#include <map>
#define PR pair<int,int>
#define MP make_pair
#define fi first
#define se second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define sqr(x) ((x)*(x))
#define LL long long
const LL INF = 1e18;
const int inf=0x3f3f3f3f;
const int M=18;
const int N=500100;
const int MOD=100000;
const double eps=1e-10;
const double pi=acos(-1.0);
using namespace std;

int n,m,belong[N],l[N],r[N],block,num,op,x,y,z;
LL p[N];
struct node
{
    LL v;
    int id;
    bool operator < (const node t) const
    {
        if(v!=t.v) return v<t.v;
        return id<t.id;
    }
}a[N],d[N];
void build()
{
    int i,j;
    block=sqrt(n);
    num=n/block;if(n%block) num++;
    for(i=1;i<=num;i++)
        l[i]=(i-1)*block+1,r[i]=i*block;
    r[num]=n;
    for(i=1;i<=n;i++)
        belong[i]=(i-1)/block+1,d[i]=a[i];
    for(i=1;i<=num;i++)
        sort(d+l[i],d+r[i]+1);
}
void update(int L,int R,int z)
{
    int i,ll=belong[L],rr=belong[R];
    if(ll==rr)
    {
        for(i=L;i<=R;i++) a[i].v+=z;
        for(i=l[ll];i<=r[ll];i++) d[i]=a[i];
        sort(d+l[ll],d+r[ll]+1);
        return;
    }
    for(i=L;i<=r[ll];i++) a[i].v+=z;
    for(i=l[ll];i<=r[ll];i++) d[i]=a[i];
    sort(d+l[ll],d+r[ll]+1);
    for(i=ll+1;i<rr;i++) p[i]+=z;
    for(i=l[rr];i<=R;i++) a[i].v+=z;
    for(i=l[rr];i<=r[rr];i++) d[i]=a[i];
    sort(d+l[rr],d+r[rr]+1);
}
int bisearch1(int L,int R,LL x)
{
    int l=L,r=R;
    while(l<r)
    {
        int m=l+(r-l)/2;
        if(d[m].v<x) l=m+1;
        else r=m;
    }
    return d[r].v==x?d[r].id:-1;
}
int bisearch2(int L,int R,LL x)
{
    int l=L,r=R;
    while(l<r)
    {
        int m=l+(r-l+1)/2;
        if(d[m].v<=x)l=m;
        else r=m-1;
    }
    return d[l].v==x?d[l].id:-1;
}
int query(int x)
{
    int i,j,lpos=-1,rpos=-1;
    for(i=1;i<=num;i++)
    {
        if(lpos==-1) lpos=bisearch1(l[i],r[i],x-p[i]);
        rpos=max(rpos,bisearch2(l[i],r[i],x-p[i]));
    }
    if(lpos==-1||rpos==-1) return -1;
    else return rpos-lpos;
}

int main()
{
    int i;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)scanf("%I64d",&a[i].v),a[i].id=i;
    build();
    while(m--)
    {
        scanf("%d",&op);
        if(op==1) scanf("%d%d%d",&x,&y,&z),update(x,y,z);
        else scanf("%d",&x),printf("%d\n",query(x));
    }
}