#数列分块入门 2

这道题他们说不用分块就用平衡树，所以暴力拆解就用分块。

做了两道分块已经渐渐明白分块的一些模板在这里就放一下吧（用vector比用数组md方便多了）

 1 void /*修改操作*//*查询操作*/ (int l,int r.....,int val)/*修改的要素*/
 2 {
 3     /*总是先修改（查询）角区间，就是那种不是整块整块的*/
 4     for(int i=l;i<=(id[l],n)/*注意结束的区间可能就是 l~r 所以取min*/;i++) 
 5     {
 6         /*暴力枚举修改（查询）*/
 7     } 
 8     if(/* l 和 r 不在同一个区间 就把 r 所在的角区间枚举修改（查询）*/)
 9     {
10         /*暴力~*/
11         /*但是能二分或者别的什么能搞出来，那肯定换一下啦*/
12     }
13     /*最后就整块整块的枚举啦*/ 
14 }

这道题就是维护每个块内的元素成单调的，这样最后查询的时候就可以用二分（lower_bound）来进行查询会快很多,而剩下的就只有暴力、暴力再暴力。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 int n;
 4 int block;
 5 int id[50009];
 6 int num[50009];
 7 int tag[10009];
 8 vector<int>v[10009];
 9 inline long long kd()
10 {
11     long long x=0,f=1;char ch=getchar();
12     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
13     while(ch>='0'&&ch<='9'){x=x*10+(ch^48);ch=getchar();}
14     return x*f;
15 }
16 inline void reset(int now)
17 {
18     int r=min(id[now]*block,n);
19     v[id[now]].clear();/*清空*/
20     for(int i=(id[now]-1)*block+1;i<=r;i++)v[id[i]].push_back(num[i]);/*挼进去*/
21     sort(v[id[now]].begin(),v[id[now]].end());/*排个序*/
22 }
23 inline void add(int l,int r,int c)
24 {
25     int minn=min(r,id[l]*block);
26     for(int i=l;i<=minn;i++)num[i]+=c;/*暴力修改角区间*/
27     reset(l);/*同下面reset的注释*/
28     if(id[l]!=id[r])/*如果 l 和 r 不在个区间之内*/
29     {
30         for(int i=(id[r]-1)*block+1;i<=r;i++)num[i]+=c;/*同样暴力*/
31         reset(r);/*角区间进行修改之后就打乱的顺序，所以重新排序*/
32     }
33     for(int i=id[l]+1;i<=id[r]-1;i++)tag[i]+=c;/*直接用tag来修改整块区间*/
34 }
35 int query(int l,int r,int c)
36 {
37     int ans=0;
38     int minn=min(r,id[l]*block);
39     for(int i=l;i<=minn;i++)
40         if(num[i]+tag[id[i]]<c)ans++;/*暴力枚举角区间*/
41     if(id[l]!=id[r])
42     {
43         for(int i=(id[r]-1)*block+1;i<=r;i++)
44         {
45             if(num[i]+tag[id[i]]<c)ans++;/*再次暴力枚举*/
46         }
47     }
48     for(int i=id[l]+1;i<=id[r]-1;i++)
49     {
50         int x=c-tag[i];
51         ans+=lower_bound(v[i].begin(),v[i].end(),x)-v[i].begin();/*整块的就直接二分查找*/
52     }
53     return ans;
54 }
55 int main()
56 {
57     n=kd();
58     for(int i=1;i<=n;i++)num[i]=kd();
59     block=(int)sqrt(n);
60     for(int i=1;i<=n;i++)
61     {
62         id[i]=(i-1)/block+1;
63         v[id[i]].push_back(num[i]);
64     }
65     for(int i=1;i<=id[n];i++)
66         sort(v[i].begin(),v[i].end());
67     for(int i=1;i<=n;i++)
68     {
69         int opt=kd();
70         int l=kd(),r=kd(),c=kd();
71         if(opt==0)add(l,r,c);
72         if(opt==1)cout<<query(l,r,c*c)<<endl;
73     }
74     return 0;
75 }

猜你喜欢