看了电子科大的分块算法视频学会的,看视频来学算法感觉特别棒!!
当线段树或树状数组对区间问题还是很难解决时,可以考虑分块。
分块的思想是把范围为n的数列,分为num块,每块长度为sqrt(n),长度为num(n%sqrt(n)==0时)/num+1块。
block为块的长度,然后定义数组L[i],R[i],表示块i的左边界和右边界。belong[i]意思是第i个数据属于哪一个块。
这样就build结束,模版:
void bulid(){
block=sqrt(n);
num=n/block;if(n%block!=0)num++; //block为块的长度,num为块的数量
for(int i;i<=num;i++){ //l,r为这个块的左边界与右边界
l[i]=(i-1)*block+1;
r[i]=i*block;
}
r[num]=n;
单点更新时在改变其点的值的同时改变这个块的最值。
int update(int x,int y){
a[x]+=y;
MAX[belong[x]]=max(MAX[belong[x]],a[x]);//更新这个数的同时更新这个块
}
询问线段[x,y]最值时
就可以暴力x到r[belong[x]]的最值,然后是块状之间的最值l[belong[x]]+1 to r[belong[y]]-1 之间的最值,
再从l[belong[y]] to y 便可查询到最值
long long ask(int x,int y){
long long res=0;
if(belong[x]==belong[y]){
for(int i=x;i<=y;i++)
res=max(res,a[i]);
return res;
}
for(int i=x;i<=r[belong[x]];i++)
res=max(res,a[i]);
for(int i=belong[x]+1;i<=belong[y]-1;i++)
res=max(res,MAX[i]);
for(int i=l[belong[y]];i<=y;i++)
res=max(res,a[i]);
return res;
}
卿学姐与公主
Time Limit: 2000/1000MS (Java/Others) Memory Limit: 65535/65535KB (Java/Others)
某日,百无聊赖的卿学姐打开了某11区的某魔幻游戏
在这个魔幻的游戏里,生活着一个美丽的公主,但现在公主被关押在了魔王的城堡中。
英勇的卿学姐拔出利刃冲向了拯救公主的道路。
走过了荒野,翻越了高山,跨过了大洋,卿学姐来到了魔王的第一道城关。
在这个城关面前的是魔王的精锐部队,这些士兵成一字排开。
卿学姐的武器每次只能攻击一个士兵,并造成一定伤害,卿学姐想知道某时刻从 L 到 R 这个区间内,从开始到现在累计受伤最严重的士兵受到的伤害。
最开始每个士兵的受到的伤害都是0
Input
第一行两个整数 N,Q 表示总共有 N 个士兵编号从 1 到 N ,和 Q 个操作。
接下来 Q 行,每行三个整数,首先输入一个 t ,如果 t 是 1 ,那么输入 p,x ,表示卿学姐攻击了 p 这个位置的士兵,并造成了 x 的伤害。如果 t 是 2 ,那么输入 L,R ,表示卿学姐想知道现在 [L,R] 闭区间内,受伤最严重的士兵受到的伤害。
1≤N≤100000
1≤Q≤100000
1≤p≤N
1≤x≤100000
1≤L≤R≤N
Output
对于每个询问,回答相应的值
Sample input and output
Sample Input | Sample Output |
---|---|
5 4 2 1 2 1 2 4 1 3 5 2 3 3 | 0 5 |
Hint
注意可能会爆int哦
#include<iostream>
#include <stdio.h>
#include <string.h>
#include<stdlib.h>
#include <map>
#include <set>
#include <math.h>
#include <queue>
#include <algorithm>
#include<vector>
using namespace std;
const int maxn= 1e5+7;
int belong[maxn],num,l[maxn],r[maxn];
long long a[maxn],MAX[maxn],q,n;
int block;
void bulid(){
block=sqrt(n);
num=n/block;if(n%block!=0)num++; //block为块的长度,num为块的数量
for(int i;i<=num;i++){ //l,r为这个块的左边界与右边界
l[i]=(i-1)*block+1;
r[i]=i*block;
}
r[num]=n;
for(int i=1;i<=n;i++) //每一个数,其属于哪一个块
belong[i]=(i-1)/block+1;
for(int i=1;i<=num;i++) //上面部分是模版,下面这行开始是根据题意而初始化的
for(int j=l[i];j<=r[i];j++)
MAX[i]=max(MAX[i],a[j]);
}
int update(int x,int y){
a[x]+=y;
MAX[belong[x]]=max(MAX[belong[x]],a[x]);//更新这个数的同时更新这个块
}
long long ask(int x,int y){
long long res=0;
if(belong[x]==belong[y]){
for(int i=x;i<=y;i++)
res=max(res,a[i]);
return res;
}
for(int i=x;i<=r[belong[x]];i++)
res=max(res,a[i]);
for(int i=belong[x]+1;i<=belong[y]-1;i++)
res=max(res,MAX[i]);
for(int i=l[belong[y]];i<=y;i++)
res=max(res,a[i]);
return res;
}
int main(){
//freopen("in.txt","r",stdin);
int i,j,k,l1,f1,f2,f3,f4,f5,t4,t1,m;
scanf("%d %d",&n,&q);
for(i=1;i<=q;i++){
scanf("%d %d %d",&f1,&f2,&f3);
if(f1==1)update(f2,f3);
else printf("%lld\n",ask(f2,f3));
}
return 0;
}