分块学习笔记

学了分块，感觉这玩意好难啊，怎么听起来这么简单？【】【】分块！

先推荐一个东西：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

直到我想起一道题，但是忘了是哪道。

1 hour later...

（把脑袋往墙上撞的声音）就这么点数据范围（也许 \(a_i\) 范围上到 LONG_LONG_MAX 都可以做！），根号不是几下就废了吗？（变成 \(1\) 或者 \(0\)）

没算是几下，一定是不超过 \(10\) 下，够了。

所有部分，包括两头散块和中间整块，都可以用一个 sol 函数解决。

这个函数干嘛呢？

把需要处理的块内部还可以开方的开方。

好了就是这么简单。

本题代码

#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;
}

T6

我【】【】【】【】【】【】【】【】【】【】【】【】【】【】【】分块跑链表玩你【】【】【】【】【】【】—— 刚开题的我 beeeeeeeeeeee like

结果口胡了个算法后来发现是正解？

是这样的，如果我们每个块内用其它数据结构，是能够支持其它不一样的操作的，比如这题是 vector，每次插入时先找到位置所在的块，再暴力插入，把块内的其它元素直接向后移动一位。

当然你想用链表也是可以的。

等会，万一复杂度炸了咋办？比如某个块你疯狂往里面插入。

那就重构我们的分块啊！

每 \(\sqrt n\) 次插入后，重新把数列分一下块（像最开始那样），重构需要的复杂度为 \(O(n)\)，重构的次数最多 \(\sqrt n\)，所以复杂度没有问题。

结果，脑子：哇真的可以！手：怎么玩？！

好吧，菜就多练（对自己说）。

对了我换了个机子结果把原来的模板忘了所以叕重构了代码。

本题代码

#include<stdio.h>
#include<bits/stdc++.h>
#define N 200010
#define MOD 998244353
#define esp 1e-8
#define INF 999999999999999999
#define LL int
#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)))
#define lc (u<<1)
#define rc (u<<1|1)
using namespace std;
i128 read()
{
	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;
	return x;
}
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,m,bl,a[100005],s[200005],tp;
vector<LL>v[1005];
pll q(LL b)
{
    LL x=1;
    while(b>v[x].size())
    {   	
        b-=v[x].size();
		x++;
	}
    return mkp(x,b-1);
}
void reb()
{
    tp=0;
    rep(i,1,m,1)
	{
        for(vector<LL>::iterator j=v[i].begin();j!=v[i].end();j++)s[++tp]=*j;
        v[i].clear();
    }
    LL blo=sqrt(tp);
    rep(i,1,tp,1)v[(i-1)/blo+1].push_back(s[i]);
    m=(tp-1)/blo+1;
}
void ins(LL a,LL b)
{
    pll t=q(a);
    v[t.first].insert(v[t.first].begin()+t.second,b);
    if(v[t.first].size()>20*bl)reb();
}
int main()
{
	cin>>n;
	bl=sqrt(n);
    rep(i,1,n,1)cin>>a[i];
    rep(i,1,n,1)v[(i-1)/bl+1].push_back(a[i]);
    m=(n-1)/bl+1;//最开始这里写成(n-1)*bl+1了，直接爆炸！
    rep(i,1,n,1)
    {
        LL op,a,b,c;
        cin>>op>>a>>b>>c;
        if(!op)ins(a,b);
        else
        {
            pll t=q(b);
            cout<<v[t.first][t.second]<<endl; 
        }
    }
	return 0;
}

T9

教练讲了这题所以我先做 T9 诶嘿

这题可以用莫队，但是他的双倍经验 P4168 不行。

这里我的思路是 P4168 一篇题解的，感谢其作者。

解析我觉得他说的很好了我就不打了（

本题代码

#include<stdio.h>
#include<bits/stdc++.h>
#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)))
#define lc (u<<1)
#define rc (u<<1|1)
using namespace std;
i128 read()
{
	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;
	return x;
}
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);
}
const LL N=100010,M=100010,K=1010;
LL n,m,tot,len,l,r,ans;
LL a[N],b[N],p[N],s[K][N],L[K],R[K],pp[K][K],bok[N];
LL cnt;
void lsh()
{
	len=sqrt(n);
	memset(L,0x3f,sizeof(L));
	rep(i,1,n,1)
	{
		p[i]=(i-1)/len+1;
		L[p[i]]=min(L[p[i]],i);
		R[p[i]]=max(R[p[i]],i);
		a[i]=read();
		b[i]=a[i];	
	}
	sort(b+1,b+1+n);
	tot=unique(b+1,b+1+n)-b-1;
	rep(i,1,n,1)a[i]=lower_bound(b+1,b+1+tot,a[i])-b;
}
void ycl()
{
	rep(i,1,n,1)rep(j,p[i],p[n],1)s[j][a[i]]++;
	rep(i,1,p[n],1)
	{
		memset(bok,0,sizeof(bok));
		rep(j,i,p[n],1)
		{
			pp[i][j]=pp[i][j-1];
			rep(k,L[j],R[j],1)
			{
				bok[a[k]]++;
				if((bok[a[k]]>bok[pp[i][j]])||(bok[a[k]]==bok[pp[i][j]]&&a[k]<pp[i][j]))pp[i][j]=a[k];
			}
		}
	}
	memset(bok,0,sizeof(bok));
}
LL q(LL l,LL r)
{
	if(p[r]-p[l]<=2)
	{
		LL res=0;
		rep(i,l,r,1)
		{
			bok[a[i]]++;
			if(bok[a[i]]>bok[res]||(bok[a[i]]==bok[res]&&a[i]<res))res=a[i];
		}
		rep(i,l,r,1)bok[a[i]]=0;
		return res;
	}
	LL res=0;
	rep(i,l,R[p[l]],1)if(!bok[a[i]])bok[a[i]]+=s[p[r]-1][a[i]]-s[p[l]][a[i]];
	rep(i,L[p[r]],r,1)if(!bok[a[i]])bok[a[i]]+=s[p[r]-1][a[i]]-s[p[l]][a[i]];
	rep(i,l,R[p[l]],1)
	{
		bok[a[i]]++;
		if(bok[a[i]]>bok[res]||(bok[a[i]]==bok[res]&&a[i]<res))res=a[i];
	}
	rep(i,L[p[r]],r,1)
	{
		bok[a[i]]++;
		if(bok[a[i]]>bok[res]||(bok[a[i]]==bok[res]&&a[i]<res))res=a[i];
	}
	LL k=pp[p[l]+1][p[r]-1];
	LL lxl=s[p[r]-1][k]-s[p[l]][k];
	rep(i,l,R[p[l]],1)lxl+=(a[i]==k);
	rep(i,L[p[r]],r,1)lxl+=(a[i]==k);
	if(lxl>bok[res]||(lxl==bok[res]&&k<res))res=k;
	rep(i,l,R[p[l]],1)bok[a[i]]=0;
	rep(i,L[p[r]],r,1)bok[a[i]]=0;
	return res;
}
int main()
{
	cin>>n;
	lsh();
	ycl();
	m=n;
	rep(i,1,m,1)
	{
		l=read();
		r=read();
		if(l>r)swap(l,r);
		cout<<(ans=b[q(l,r)])<<endl;
	}
	return 0;
}

T7

感谢 @Lixiang_is_potato 和她的大号 @MinimumSpanningTree（请勿在没有事的情况下向其大号发送任何私信或 @）提供的思路，我要给她发 —— 奖 —— 状！

想不到 yeeeeeeee 点怎么做

突然看到了她发的求助帖：捞帖和原帖。

~~她改完错以后还把原贴删了，我好感动啊~~

帮她改了错后突然发现这玩意好简单…

I have a problem（可以借鉴这题 tag 的思路），

I have a 分块！

A————————C！

好了不闹了正经的。

搞俩 tag：T1 那种加的，和乘的。

处理乘修改的时候，乘标记要乘，记得把加的也乘了。

加法只有加标记要加。

然后对于散块，我们可以不管直接暴力更新！

同学的代码，简洁、明了

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
typedef long long ll;
const int N=100100,M=1100;
const ll MOD=10007;
int n,op,x,y,block,nk,st[M],ed[M],p[N];
ll a[N],sum[M],add[M],c;
void build()
{
	block=sqrt(n);
	nk=n/block;
	if(n%block!=0) nk++;
	for(int i=1;i<=nk;i++) st[i]=(i-1)*block+1,ed[i]=st[i]+block-1,sum[i]=1;
	ed[nk]=n;
	for(int i=1;i<=n;i++) p[i]=(i-1)/block+1;
}
void hy(int k)
{
	for(int i=st[k];i<=ed[k];i++) a[i]=(a[i]*sum[k]+add[k])%MOD;
	sum[k]=1,add[k]=0;
}
void update_j(int l,int r,ll num)
{
	int kl=p[l],kr=p[r];
	if(kl==kr)
	{
		hy(kl);
		for(int i=l;i<=r;i++) a[i]+=num,a[i]%=MOD;
	}
	else
	{
		for(int i=kl+1;i<=kr-1;i++) add[i]+=num,add[i]%=MOD;
		hy(kl);
		for(int i=l;i<=ed[kl];i++) a[i]+=num,a[i]%=MOD;
		hy(kr);
		for(int i=st[kr];i<=r;i++) a[i]+=num,a[i]%=MOD;
	}
}
void update_c(int l,int r,ll num)
{
	int kl=p[l],kr=p[r];
	if(kl==kr)
	{
		hy(kl);
		for(int i=l;i<=r;i++) a[i]*=num,a[i]%=MOD;
	}
	else
	{
		for(int i=kl+1;i<=kr-1;i++) sum[i]*=num,sum[i]%=MOD,add[i]*=num,add[i]%=MOD;
		hy(kl);
		for(int i=l;i<=ed[kl];i++) a[i]*=num,a[i]%=MOD;
		hy(kr);
		for(int i=st[kr];i<=r;i++) a[i]*=num,a[i]%=MOD;
	}
}
ll query(int l)
{
	int k=p[l];
	ll ans=(a[l]*sum[k]+add[k])%MOD;
	return ans;
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%lld",&a[i]),a[i]%=MOD;
	build();
	for(int i=1;i<=n;i++)
	{
		scanf("%d%d%d%lld",&op,&x,&y,&c);
		c%=MOD;
		if(!op) update_j(x,y,c);
		else if(op==1) update_c(x,y,c);
		else printf("%lld\n",query(y));
	}
	return 0;
}

我的石山唐诗代码

#include<stdio.h>
#include<bits/stdc++.h>
#define N 100010
#define M 1010
#define MOD 10007
#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)))
#define lc (u<<1)
#define rc (u<<1|1)
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,op,x,y,kc,ks,L[M],R[M],whe[N];
LL a[N],mu[M],ad[M],c;
LL q(LL l)
{
	return (a[l]*mu[whe[l]]+ad[whe[l]])%MOD;
}
void baoli(LL k)
{
	rep(i,L[k],R[k],1)a[i]=(a[i]*mu[k]+ad[k])%MOD;
	mu[k]=1;
	ad[k]=0;
}
void bui()
{
	kc=sqrt(n);
	ks=n/kc;
	if(n%kc!=0)ks++;
	rep(i,1,ks,1)
	{
		L[i]=(i-1)*kc+1;
		R[i]=L[i]+kc-1;
		mu[i]=1;
	}
	R[ks]=n;
	rep(i,1,n,1)whe[i]=(i-1)/kc+1;
}
void add(LL l,LL r,LL num)
{
	LL kl=whe[l],kr=whe[r];
	if(kl==kr)
	{
		baoli(kr);
		rep(i,l,r,1)
		{
			a[i]=(a[i]+num)%MOD;
		}
	}
	else
	{
		baoli(kl);
		baoli(kr);
		repn(i,kl+1,kr,1)
		{
			ad[i]=(ad[i]+num)%MOD;
		}
		rep(i,l,R[kl],1)
		{
			a[i]=(a[i]+num)%MOD;
		}
		rep(i,L[kr],r,1)
		{
			a[i]=(a[i]+num)%MOD;
		}
	}
}
void tim(LL l,LL r,LL num)
{
	LL kl=whe[l],kr=whe[r];
	if(kl==kr)
	{
		baoli(kr);
		rep(i,l,r,1)
		{
			a[i]=(a[i]*num)%MOD;
		}
	}
	else
	{
		repn(i,kl+1,kr,1)
		{
			ad[i]=(ad[i]*num)%MOD;
			mu[i]=(mu[i]*num)%MOD;
		}
		baoli(kl);
		baoli(kr);
		rep(i,l,R[kl],1)
		{
			a[i]=(a[i]*num)%MOD;
		}
		rep(i,L[kr],r,1)
		{
			a[i]=(a[i]*num)%MOD;
		}
	}
}
int main()
{
	cin>>n;
	rep(i,1,n,1)
	{
	    read(a[i]);
		a[i]%=MOD;
	}
	bui();
	rep(i,1,n,1)
	{
	    read(op);
	    read(x);
	    read(y);
	    read(c);
		c%=MOD;
		if(op==0)add(x,y,c);
		if(op==1)tim(x,y,c);
		if(op==2)
		{
    		write(q(y));
    		puts("");
		}
	}
	return 0;
}

T8

最开始的做法会 T 飞起，寄。

以下是另一位同学教我的思路：

修改的话，对于整块，打标记。

散块就暴力改，如果有标记先把标记更新再改。

然后是查询。

对于整块，有标记的直接统计，没有就暴力。

我们的散块还是暴力。

标记就是当前这一块统一的值。

然后我这个飞舞在标记上调了很久…

本题代码

#include<stdio.h>
#include<bits/stdc++.h>
#define N 100010
#define M 1010
#define MOD 10007
#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)))
#define lc (u<<1)
#define rc (u<<1|1)
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,x,y,c,kc,ks,L[M],R[M],whe[N],a[N],tg[N];
void bui()
{
	kc=sqrt(n);
	ks=n/kc;
	if(n%kc!=0)ks++;
	rep(i,1,ks,1)
	{
		L[i]=(i-1)*kc+1;
		R[i]=L[i]+kc-1;
		tg[i]=INF;
	}
	R[ks]=n;
	rep(i,1,n,1)whe[i]=(i-1)/kc+1;
}
LL q(LL x,LL y,LL c)
{
	LL lin=whe[x],rin=whe[y],sum=0;
	rep(i,lin,rin,1)
	{
		if(tg[i]!=INF&&tg[i]==c)sum+=min(R[i],y)-max(x,L[i])+1;
		else
		{
			if(tg[i]!=INF)continue;
			rep(j,max(x,L[i]),min(R[i],y),1)
			{
				if(a[j]==c)sum++;
			}
		}
	}
	return sum;
}
void chg(LL x,LL y,LL c)
{
	LL lin=whe[x],rin=whe[y];
	if(tg[lin]!=INF)
	{
		rep(i,L[lin],R[lin],1)
		{
			a[i]=tg[lin];
		}
		tg[lin]=INF;
	}
	if(tg[rin]!=INF)
	{
		rep(i,L[rin],R[rin],1)
		{
			a[i]=tg[rin];
		}
		tg[rin]=INF;
	}
	if(lin==rin)
	{
		rep(i,x,y,1)
		{
			a[i]=c;
		}
	}
	else
	{
		rep(i,x,R[lin],1)a[i]=c;
		rep(i,L[rin],y,1)a[i]=c;
		repn(i,lin+1,rin,1)tg[i]=c;
	}
}
void sol(LL x,LL y,LL c)
{
	cout<<q(x,y,c)<<endl;
	chg(x,y,c);
}
int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);
	cin>>n;
	rep(i,1,n,1)
	{
		cin>>a[i];
	}
	bui();
	rep(i,1,n,1)
	{
		cin>>x>>y>>c;
		sol(x,y,c);
	}
	return 0;
}

（加餐）T10

题外话：这玩意要不是 LCT 我觉得评绿。

\(0\sim n−1\)！！

挂分现场，没看清这【】【】数据范围。

然后我的思路是：一个数组 tp（题外话：这是我的世界里面传送的指令），这个位置表示跳出该块后的位置。

和一个数组 se（别想歪了），这个位置表示跳出这个块要几步。

跳的话就暴力跳。

水题，真的不配紫。

本题代码

#include<stdio.h>
#include<bits/stdc++.h>
#define N 200010
#define M 1010
#define MOD 10007
#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)))
#define lc (u<<1)
#define rc (u<<1|1)
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,m,op,x,y,kc,ks,L[M],R[M],whe[N];
LL a[N],se[N],tp[N];
void chg(LL l,LL r)
{
	rem(i,r,l,1)
	{
		if(i+a[i]>R[whe[i]])
		{
			tp[i]=i+a[i];
			se[i]=1;
		}
		else
		{
			tp[i]=tp[i+a[i]];
			se[i]=se[i+a[i]]+1;
		}
	}
}
LL baoli(LL x)
{
	LL sum=0;
	while(x<=n)
	{
		sum+=se[x];
		x=tp[x];
	}
	return sum;
}
void bui()
{
	kc=sqrt(n);
	ks=n/kc;
	if(n%kc!=0)ks++;
	rep(i,1,ks,1)
	{
		L[i]=(i-1)*kc+1;
		R[i]=L[i]+kc-1;
	}
	R[ks]=n;
	rep(i,1,n,1)whe[i]=(i-1)/kc+1;
	chg(1,n);
}
int main()
{
	cin>>n;
	rep(i,1,n,1)
	{
		cin>>a[i];
	}
	bui();
	cin>>m;
	while(m--)
	{
		cin>>op>>x;
		x++;
		if(op==1)
		{
			cout<<baoli(x)<<endl;
		}
		else
		{
			cin>>y;
			a[x]=y;
			chg(L[whe[x]],R[whe[x]]);
		}
	}
	return 0;
}

小结

题外话：你们可能看着这个博客很短，但我写到后面都卡了…

个人对分块的总结：简洁、复杂度较优、可以在线、后面可以用这个思想写莫队。

当然莫队肯定也要学的但是我发挥传统异能咕咕咕啦！

分块思考比较简单，代码也简单，但对我这个蒟蒻就 emm 了…

posted @ 2024-03-30 16:04 cppom 阅读(138) 评论(0) 编辑收藏举报

刷新页面返回顶部

登录后才能查看或发表评论，立即登录或者逛逛博客园首页

相关博文：

· CSP-S 2023 题解

· CSP-J 2023 题解

· 数列分块入门1-9

· 分块/莫队总结（未完结）

· 分块学习笔记

阅读排行：
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型，支持深度思考和联网搜索！
· 使用C#创建一个MCP客户端
· 基于 Docker 搭建 FRP 内网穿透开源项目（很简单哒）
· ollama系列1：轻松3步本地部署deepseek，普通电脑可用
· 按钮权限的设计及实现

公告

昵称： cppom
园龄： 1年
粉丝： 3
关注： 9

+加关注

2025年3月

日

一

二

三

四

五

六

飞舞 OIer cppom

世末 OIer

分块学习笔记

T1

T2

T3

T4

T5

T6

T9

T7

T8

（加餐）T10

小结

公告

搜索

常用链接

最新随笔

积分与排名

随笔分类 (13)

随笔档案 (13)

阅读排行榜

评论排行榜

推荐排行榜

最新评论