分块-教主的魔法
https://www.luogu.org/problem/show?pid=2801#sub
思路就是分块嘛;
第一次分块题啊,那讲的细一点吧;
首先堆这些数分块;
这里我的长度是sqrt(n);
分块的长度根据题目而定;
往往推出一个公式,然后根据这个公式求出最优的长度;
这个就是数学范畴;
我的strcut
l,r是区间范围,len是长度(在这一题里无用),tag是加法标记
struct kuai{
int l,r,len,tag;
}c[1010];
分块的思想拉,就是区间内包含整个块的就整个块地操作,
不是整个快的最多就两个部分,操作区间的左右两端,直接暴力;
我的主程序
a[]是原始数列,不改变;
b[]是分块序列,所有操作都在b[]进行;
int a[1000001],b[1000001];
int ll,n,m,x,y,z,ans;
char cc;
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%d",&a[i]),b[i]=a[i];
fenkuai();
while(m--){
cin>>cc;
scanf("%d%d%d",&x,&y,&z);
if(cc=='A')outit(x,y,z);else change(x,y,z);
}
}
这种复杂的题目,主程序轻巧一点,而且一定要写好一个模块就调试一遍,千万要小心,不然全打好了调试得要死;
我的分块函数
ll就是快的数量
void fenkuai(){
int k=sqrt(n);
int l=1;
while(l<=n){
c[++ll].l=l;
c[ll].r=min(n,l+k-1);
c[ll].len=c[ll].r-c[ll].l+1;
sort(b+c[ll].l,b+c[ll].r+1);
l+=k;
}
}
我们将每个块都排序;这样就方便二分找答案;
findans的函数
就是把x,y分为左中右三个区间
void outit(int x,int y,int z){
int ans=0;
for(int i=1;i<=ll;i++)
if(x<=c[i].l&&c[i].r<=y)ans+=er(i,z);else
if(c[i].l<=x&&x<=c[i].r)ans+=BL(i,x,c[i].r,z);else
if(c[i].l<=y&&y<=c[i].r)ans+=BL(i,c[i].l,y,z);
printf("%d\n",ans);
}
二分
int er(int num,int z){
int l=c[num].l,r=c[num].r,ans=1e9;
while(r>=l){
int mid=(l+r)/2;
if(b[mid]+c[num].tag>=z){//要加上 c[num].tag
ans=min(ans,mid);r=mid-1;
}else l=mid+1;
}
if(ans==1e9)return 0;//这个特判容易忘
return c[num].r-ans+1;
}
暴力
暴力部分要在a数组里找
int BL(int num,int x,int y,int z){
int ans=0;
for(int i=x;i<=y;i++)if(a[i]+c[num].tag>=z)ans++;
return ans;
}
然后就是改变部分
void change(int x,int y,int z){
for(int i=1;i<=ll;i++)
if(x<=c[i].l&&c[i].r<=y)c[i].tag+=z;else
if(c[i].l<=x&&x<=c[i].r)BLchange(i,x,c[i].r,z);else
if(c[i].l<=y&&y<=c[i].r)BLchange(i,c[i].l,y,z);
}
void BLchange(int num,int x,int y,int z){
for(int i=c[num].l;i<=c[num].r;i++){
b[i]=a[i]+c[num].tag;
if(x<=i&&i<=y)b[i]+=z;
}
sort(b+c[num].l,b+c[num].r+1);
}
总得来说,复杂的题目要小心,思路明确,打好一个模块调试一个;
代码
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define Ll unsigned long long
using namespace std;
struct kuai{
int l,r,len,tag;
}c[1010];
int a[1000001],b[1000001];
int ll,n,m,x,y,z,ans;
char cc;
void fenkuai(){
int k=sqrt(n);
int l=1;
while(l<=n){
c[++ll].l=l;
c[ll].r=min(n,l+k-1);
c[ll].len=c[ll].r-c[ll].l+1;
sort(b+c[ll].l,b+c[ll].r+1);
l+=k;
}
}
int er(int num,int z){
int l=c[num].l,r=c[num].r,ans=1e9;
while(r>=l){
int mid=(l+r)/2;
if(b[mid]+c[num].tag>=z){//要加上 c[num].tag
ans=min(ans,mid);r=mid-1;
}else l=mid+1;
}
if(ans==1e9)return 0;//这个特判容易忘
return c[num].r-ans+1;
}
int BL(int num,int x,int y,int z){
int ans=0;
for(int i=x;i<=y;i++)if(a[i]+c[num].tag>=z)ans++;
return ans;
}
void outit(int x,int y,int z){
int ans=0;
for(int i=1;i<=ll;i++)
if(x<=c[i].l&&c[i].r<=y)ans+=er(i,z);else
if(c[i].l<=x&&x<=c[i].r)ans+=BL(i,x,c[i].r,z);else
if(c[i].l<=y&&y<=c[i].r)ans+=BL(i,c[i].l,y,z);
printf("%d\n",ans);
}
void BLchange(int num,int x,int y,int z){
for(int i=c[num].l;i<=c[num].r;i++){
b[i]=a[i]+c[num].tag;
if(x<=i&&i<=y)b[i]+=z;
}
sort(b+c[num].l,b+c[num].r+1);
}
void change(int x,int y,int z){
for(int i=1;i<=ll;i++)
if(x<=c[i].l&&c[i].r<=y)c[i].tag+=z;else
if(c[i].l<=x&&x<=c[i].r)BLchange(i,x,c[i].r,z);else
if(c[i].l<=y&&y<=c[i].r)BLchange(i,c[i].l,y,z);
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%d",&a[i]),b[i]=a[i];
fenkuai();
while(m--){
cin>>cc;
scanf("%d%d%d",&x,&y,&z);
if(cc=='A')outit(x,y,z);else change(x,y,z);
}
}