Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 80671 Accepted Submission(s): 34073
题目链接：HDU 1166
线段树入门题，但是由于最近似乎碰到了很多涉及到一些算法用分块优化，否则会T的题，只能再膜几下分块算法了，可惜网上的代码和讲解不是很多，找的基本都是只有理论的嘴强王者或者基本没有预处理写的很丑的代码（也许这份代码也很丑） ，搞来搞去就先拿这水题练练手好了……
$block[]$代表每一个块所储存的信息（这题当然是代表块内元素值之和sum，也可以是最大值max、min什么的）；
$unit$代表一个块的大小，$bcnt$代表块的个数；
$L_i$代表一个块$block_i$的左闭区间，$R_i$同理为右闭区间；
$belong[]$代表一个点所属于的块；
代码：
#include <stdio.h>
#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define CLR(arr,val) memset(arr,val,sizeof(arr))
#define LC(x) (x<<1)
#define RC(x) ((x<<1)+1)
#define MID(x,y) ((x+y)>>1)
typedef pair<int,int> pii;
typedef long long LL;
const double PI=acos(-1.0);
const int N=5e4+7;
int L[N],R[N],belong[N],unit,bcnt;
int n,m;
int arr[N],block[N];
void init()
{
CLR(L,0);
CLR(R,0);
CLR(belong,0);
CLR(arr,0);
CLR(block,0);
}
void buildblock()//预处理各个数组
{
unit=sqrt(n);
bcnt=n/unit;
if(n%unit)
++unit;
for (int i=1; i<=bcnt; ++i)
{
L[i]=(i-1)*unit+1;
R[i]=i*unit;
}
R[bcnt]=n;
for (int i=1; i<=n; ++i)
belong[i]=(i-1)/unit+1;
}
void update(int x,int val)//块和原来的元素都要更新
{
arr[x]+=val;
block[belong[x]]+=val;
}
int query(int l,int r)
{
int ret=0;
int i;
if(belong[l]==belong[r])//特判同一个块内
{
for (i=l; i<=r; ++i)
ret+=arr[i];
}
else//不同块：左右边暴力+中间分块
{
for (i=l; i<=R[belong[l]]; ++i)
ret+=arr[i];
for (i=belong[l]+1; i<belong[r]; ++i)
ret+=block[i];
for (i=L[belong[r]]; i<=r; ++i)
ret+=arr[i];
}
return ret;
}
int main(void)
{
int tcase,i,j,l,r,x,v;
char ops[10];
scanf("%d",&tcase);
for (int q=1; q<=tcase; ++q)
{
init();
scanf("%d",&n);
for (i=1; i<=n; ++i)
scanf("%d",&arr[i]);
buildblock();
for (i=1; i<=bcnt; ++i)
{
for (j=L[i]; j<=R[i]; ++j)
block[i]+=arr[j];
}
printf("Case %d:\n",q);
while (~scanf("%s",ops)&&strcmp(ops,"End"))
{
if(ops[0]==‘A‘)
{
scanf("%d%d",&x,&v);
update(x,v);
}
else if(ops[0]==‘S‘)
{
scanf("%d%d",&x,&v);
update(x,-v);
}
else
{
scanf("%d%d",&l,&r);
printf("%d\n",query(l,r));
}
}
}
return 0;
}
原文：http://www.cnblogs.com/Blackops/p/6124195.html