学了分块,感觉这玩意好难啊,怎么听起来这么简单?【】【】分块!
先推荐一个东西:loj 分块全家桶!
首先,把一整个数组劈成 \(\sqrt n\) 块是最优的!(当然如果你想写一个 \(114514\) 块的分块也没问题但他不优啊!)
长这样:
这样它的复杂度是:
- 预处理:\(O(n\sqrt n+q)\)
- 在线处理:\(O(q\sqrt n+n)\)
分块其实就是三层的树,每个非叶子结点的节点有 \(\sqrt n\) 个子节点。
像这样:
(第一层:整个大区间,第二层:每个块,第三层:每个位置)
然后呢?
没了。
你问咋处理?每个块的处理,两边的“散块”就暴力啊!
分块的思路很简单。
但某些毒瘤题的代码不做评价。
T1
模板。只放代码注释不放解析。
本题代码
#include<stdio.h>
#include<bits/stdc++.h>
#define N 1000010
#define MOD 998244353
#define esp 1e-8
#define INF 999999999999999999
#define LL long long
#define rep(i,a,b,g) for(LL i=a;i<=b;i+=g)
#define rem(i,a,b,g) for(LL i=a;i>=b;i-=g)
#define repn(i,a,b,g) for(LL i=a;i<b;i+=g)
#define remn(i,a,b,g) for(LL i=a;i>b;i-=g)
#define pll pair<LL,LL>
#define mkp(x,y) make_pair(x,y)
#define i128 __int128
#define lowbit(x) ((x)&(-(x)))
using namespace std;
void read(i128 &x)
{
i128 f=1;
x=0;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
x*=f;
}
void write(i128 x)
{
if(x>=10)write(x/10);
putchar(x%10+'0');
}
LL n,a[50010],lz[50010],op,l,r,c,kc;
//kc:块长(根号n),lz:类似lazytag,给整个块的标记
LL q(LL x)
{
return lz[x/kc]+a[x];
}
void ud(LL l,LL r,LL c)
{
if(l/kc==r/kc)
{
rep(i,l,r,1)
{
a[i]+=c;
}
}
else
{
rep(i,l,(l/kc+1)*kc-1,1)a[i]+=c;
rep(i,r/kc*kc,r,1)a[i]+=c;
//两边的散块
repn(i,l/kc+1,r/kc,1)lz[i]+=c;
//中间的整块
}
}
int main()
{
cin>>n;
kc=sqrt(n);
rep(i,1,n,1)cin>>a[i];
rep(i,1,n,1)
{
cin>>op>>l>>r>>c;
if(!op)
{
ud(l,r,c);
}
else
{
cout<<q(r)<<endl;
}
}
return 0;
}
T2
这道题目的基础是咋求一个数 \(c\) 在 \(\left [l,r\right ]\) 的排名。
咋办?肯定二分啊!排个序!
void px(LL x)
{
rep(i,max(1LL,x*kc),min(n,kc*(x+1)-1),1)b[i]=a[i];
sort(b+max(1LL,x*kc),b+min(n+1,(x+1)*kc));
}
为什么排序用
b[i]
而不是用a[i]
?
两边的散块:你这么说,我不存在?
你排序,和原来不一样了,散块表示:你【】【】!
剩下的有手就行。
有个坑点,记得及时排序。
本题代码
#include<stdio.h>
#include<bits/stdc++.h>
#define N 1000010
#define MOD 998244353
#define esp 1e-8
#define INF 999999999999999999
#define LL long long
#define rep(i,a,b,g) for(LL i=a;i<=b;i+=g)
#define rem(i,a,b,g) for(LL i=a;i>=b;i-=g)
#define repn(i,a,b,g) for(LL i=a;i<b;i+=g)
#define remn(i,a,b,g) for(LL i=a;i>b;i-=g)
#define pll pair<LL,LL>
#define mkp(x,y) make_pair(x,y)
#define i128 LL
#define lowbit(x) ((x)&(-(x)))
using namespace std;
void read(i128 &x)
{
i128 f=1;
x=0;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
x*=f;
}
void write(i128 x)
{
if(x>=10)write(x/10);
putchar(x%10+'0');
}
LL n,a[50010],b[50010],lz[50010],op,l,r,c,kc;
LL q(LL l,LL r,LL c)
{
LL sum=0;
if(l/kc==r/kc)
{
rep(i,l,r,1)
{
if(a[i]+lz[i/kc]<c)sum++;
}
}
else
{
rep(i,l,min(n,(l/kc+1)*kc-1),1)
{
if(a[i]+lz[i/kc]<c)sum++;
}
rep(i,r/kc*kc,r,1)
{
if(a[i]+lz[i/kc]<c)sum++;
}
repn(i,l/kc+1,r/kc,1)
{
LL l=i*kc-1,r=(i+1)*kc,mid;
while(l+1<r)
{
mid=l+r>>1;
if(b[mid]+lz[i]<c)l=mid;
else r=mid;
}
sum+=l-i*kc+1;
}
}
return sum;
}
void px(LL x)
{
rep(i,max(1LL,x*kc),min(n,kc*(x+1)-1),1)b[i]=a[i];
sort(b+max(1LL,x*kc),b+min(n+1,(x+1)*kc));
}
void ud(LL l,LL r,LL c)
{
if(l/kc==r/kc)
{
rep(i,l,r,1)
{
a[i]+=c;
}
px(l/kc);
}
else
{
rep(i,l,min(n,(l/kc+1)*kc-1),1)a[i]+=c;
rep(i,r/kc*kc,r,1)a[i]+=c;
px(l/kc);
px(r/kc);
repn(i,l/kc+1,r/kc,1)lz[i]+=c;
}
}
int main()
{
cin>>n;
kc=sqrt(n);
rep(i,1,n,1)
{
read(a[i]);
}
rep(i,0,n/kc,1)
{
px(i);
}
rep(i,1,n,1)
{
read(op);
read(l);
read(r);
read(c);
if(!op)
{
ud(l,r,c);
}
else
{
write(q(l,r,c*c));
puts("");
}
}
return 0;
}
T3
和上一题几乎一样。
就是维护排好序的序列。
对了我上题写的太石山了就重新写了一遍&改了自己的模板。
本题代码
#include<stdio.h>
#include<bits/stdc++.h>
#define N 100010
#define MOD 998244353
#define esp 1e-8
#define INF 999999999999999999
#define LL long long
#define rep(i,a,b,g) for(LL i=a;i<=b;i+=g)
#define rem(i,a,b,g) for(LL i=a;i>=b;i-=g)
#define repn(i,a,b,g) for(LL i=a;i<b;i+=g)
#define remn(i,a,b,g) for(LL i=a;i>b;i-=g)
#define pll pair<LL,LL>
#define mkp(x,y) make_pair(x,y)
#define i128 LL
#define lowbit(x) ((x)&(-(x)))
using namespace std;
void read(i128 &x)
{
i128 f=1;
x=0;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
x*=f;
}
void writing(i128 x)
{
if(x>=10)writing(x/10);
putchar(x%10+'0');
}
void write(i128 x)
{
if(x<0)
{
cout<<'-';
x=-x;
}
writing(x);
}
LL n,a[N],b[N],add[N],st[N],ed[N],pos[N],op,l,r,c,x,num;
vector<LL>sum[N];
void build()
{
x=sqrt(n);
num=n/x;
if(n%x)num++;
rep(i,1,num,1)
{
st[i]=(i-1)*x+1;
ed[i]=x*i;
}
ed[num]=n;
rep(i,1,n,1)
{
pos[i]=(i-1)/x+1;
sum[pos[i]].push_back(a[i]);
}
rep(i,1,num,1)sort(sum[i].begin(),sum[i].end());
}
void change(LL l,LL r,LL k)
{
LL p=pos[l],q=pos[r];
if(p==q)
{
rep(i,l,r,1)
{
a[i]+=k;
}
sum[p].clear();
rep(i,st[p],ed[p],1)
{
sum[p].push_back(a[i]);
}
sort(sum[p].begin(),sum[p].end());
return;
}
repn(i,p+1,q,1)
{
add[i]+=k;
}
rep(i,l,ed[p],1)
{
a[i]+=k;
}
sum[p].clear();
rep(i,st[p],ed[p],1)
{
sum[p].push_back(a[i]);
}
sort(sum[p].begin(),sum[p].end());
rep(i,st[q],r,1)
{
a[i]+=k;
}
sum[q].clear();
rep(i,st[q],ed[q],1)
{
sum[q].push_back(a[i]);
}
sort(sum[q].begin(),sum[q].end());
}
LL ask(LL l,LL r,LL k)
{
LL ans=-1,p=pos[l],q=pos[r];
if(p==q)
{
rep(i,l,r,1)if(a[i]+add[p]<k)ans=max(ans,a[i]+add[p]);
return ans;
}
repn(i,p+1,q,1)
{
LL tt=lower_bound(sum[i].begin(),sum[i].end(),k-add[i])-sum[i].begin();
if(tt==0)continue;
ans=max(ans,add[i]+sum[i][tt-1]);
}
rep(i,l,ed[p],1)if(a[i]+add[p]<k)ans=max(ans,a[i]+add[p]);
rep(i,st[q],r,1)if(a[i]+add[q]<k)ans=max(ans,a[i]+add[q]);
return ans;
}
int main()
{
cin>>n;
rep(i,1,n,1)
{
read(a[i]);
}
build();
rep(i,1,n,1)
{
read(op);
read(l);
read(r);
read(c);
if(!op)
{
change(l,r,c);
}
else
{
write(ask(l,r,c));
puts("");
}
}
return 0;
}
T4
水题,维护一段块的和。
本题代码
#include<stdio.h>
#include<bits/stdc++.h>
#define N 1000010
#define MOD 998244353
#define esp 1e-8
#define INF 999999999999999999
#define LL long long
#define rep(i,a,b,g) for(LL i=a;i<=b;i+=g)
#define rem(i,a,b,g) for(LL i=a;i>=b;i-=g)
#define repn(i,a,b,g) for(LL i=a;i<b;i+=g)
#define remn(i,a,b,g) for(LL i=a;i>b;i-=g)
#define pll pair<LL,LL>
#define mkp(x,y) make_pair(x,y)
#define i128 LL
#define lowbit(x) ((x)&(-(x)))
using namespace std;
void read(i128 &x)
{
i128 f=1;
x=0;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
x*=f;
}
void write(i128 x)
{
if(x>=10)write(x/10);
putchar(x%10+'0');
}
LL n,a[50010],lz[50010],ss[50010],op,l,r,c,kc;
LL q(LL l,LL r,LL c)
{
LL sum=0;
if(l/kc==r/kc)
{
rep(i,l,r,1)
{
sum+=a[i]%c+lz[i/kc]%c;
sum%=c;
}
}
else
{
rep(i,l,(l/kc+1)*kc-1,1)
{
sum+=a[i]%c+lz[i/kc]%c;
sum%=c;
}
rep(i,r/kc*kc,r,1)
{
sum+=a[i]%c+lz[i/kc]%c;
sum%=c;
}
repn(i,l/kc+1,r/kc,1)
{
sum+=ss[i];
sum%=c;
}
}
return sum;
}
void ud(LL l,LL r,LL c)
{
if(l/kc==r/kc)
{
rep(i,l,r,1)
{
a[i]+=c;
ss[i/kc]+=c;
}
}
else
{
rep(i,l,(l/kc+1)*kc-1,1)
{
a[i]+=c;
ss[i/kc]+=c;
}
rep(i,r/kc*kc,r,1)
{
a[i]+=c;
ss[i/kc]+=c;
}
repn(i,l/kc+1,r/kc,1)
{
lz[i]+=c;
ss[i]+=c*kc;
}
}
}
int main()
{
cin>>n;
kc=sqrt(n);
rep(i,1,n,1)read(a[i]);
rep(i,1,n,1)ss[i/kc]+=a[i];
rep(i,1,n,1)
{
read(op);
read(l);
read(r);
read(c);
if(!op)
{
ud(l,r,c);
}
else
{
write(q(l,r,c+1));
puts("");
}
}
return 0;
}
T5
啊?根号可以用分块吗?——刚开题的时候的我beeeeeeeeeee like
直到我想起一道题,但是忘了是哪道。
你想啊,就这么点数据范围(也许 \(a_i\) 范围上到 LONG_LONG_MAX
都可以做!),根号不是几下就废了吗?(变成 \(1\) 或者 \(0\))
没算是几下,应该是不超过 \(10\) 下,够了。
所有部分,包括两头散块和中间整块,都可以用一个 sol
函数解决。
这个函数干嘛呢?
把需要处理的块内部还可以开方的开方。文章来源:https://www.toymoban.com/news/detail-839329.html
好了就是这么简单。文章来源地址https://www.toymoban.com/news/detail-839329.html
本题代码
#include<stdio.h>
#include<bits/stdc++.h>
#define N 1000010
#define MOD 998244353
#define esp 1e-8
#define INF 999999999999999999
#define LL long long
#define rep(i,a,b,g) for(LL i=a;i<=b;i+=g)
#define rem(i,a,b,g) for(LL i=a;i>=b;i-=g)
#define repn(i,a,b,g) for(LL i=a;i<b;i+=g)
#define remn(i,a,b,g) for(LL i=a;i>b;i-=g)
#define pll pair<LL,LL>
#define mkp(x,y) make_pair(x,y)
#define i128 LL
#define lowbit(x) ((x)&(-(x)))
using namespace std;
void read(i128 &x)
{
i128 f=1;
x=0;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
x*=f;
}
void writing(i128 x)
{
if(x>=10)writing(x/10);
putchar(x%10+'0');
}
void write(i128 x)
{
if(x<0)
{
cout<<'-';
x=-x;
}
writing(x);
}
LL n,a[50010],ss[50010],op,l,r,c,kc;
vector<LL>b[50010];
LL q(LL l,LL r)
{
LL sum=0;
if(l/kc==r/kc)
{
rep(i,l,r,1)
{
sum+=a[i];
}
}
else
{
rep(i,l,(l/kc+1)*kc-1,1)
{
sum+=a[i];
}
rep(i,r/kc*kc,r,1)
{
sum+=a[i];
}
repn(i,l/kc+1,r/kc,1)
{
sum+=ss[i];
}
}
return sum;
}
void sol(LL l,LL r)
{
LL po=l/kc;
repn(i,0,b[po].size(),1)
{
LL j=b[po][i];
if(j>=l&&j<=r)
{
ss[po]-=a[j];
a[j]=sqrt(a[j]);
ss[po]+=a[j];
if(a[j]<=1)
{
b[po].erase(b[po].begin()+i--);
}
}
}
}
void ud(LL l,LL r)
{
if(l/kc==r/kc)
{
sol(l,r);
}
else
{
sol(l,(l/kc+1)*kc-1);
sol(r/kc*kc,r);
repn(i,l/kc+1,r/kc,1)
{
sol(i*kc,(i+1)*kc-1);
}
}
}
int main()
{
cin>>n;
kc=sqrt(n);
rep(i,1,n,1)read(a[i]);
rep(i,1,n,1)
{
b[i/kc].push_back(i);
ss[i/kc]+=a[i];
}
rep(i,1,n,1)
{
read(op);
read(l);
read(r);
read(c);
if(!op)
{
ud(l,r);
}
else
{
write(q(l,r));
puts("");
}
}
return 0;
}
到了这里,关于分块学习笔记的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!