ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

二维树状数组基本操作

2021-06-11 20:01:45  阅读:163  来源: 互联网

标签:树状 sum 二维 数组 基本操作 ll getchar


单点修改区间查询

LOJ #133. 二维树状数组 1:单点修改,区间查询

根据二维前缀和的思想对普通树状数组优化:

const int N = 4106;

inline ll Read()
{
	ll x = 0, f = 1;
	char c = getchar();
	while (c != '-' && (c < '0' || c > '9')) c = getchar();
	if (c == '-') f = -f, c = getchar();
	while (c >= '0' && c <= '9') x = (x << 3) + (x << 1) + c - '0', c = getchar();
	return x * f;
}

int n, m;
ll t[N][N];

void modify(int x, int y, ll val)
{
	for (int i = x; i <= n; i += i & -i)
		for (int j = y; j <= m; j += j & -j)
			t[i][j] += val;
}

ll query(int x, int y)
{
	ll ans = 0;
	for (int i = x; i; i -= i & -i)
		for (int j = y; j; j -= j & -j)
			ans += t[i][j];
	return ans;
}

int main()
{
	n = Read(), m = Read(); 
	for (int op, a, b, c, d; scanf ("%d", &op) != EOF; )
	{
		if(op == 2) a = Read(), b = Read(), c = Read(), d = Read(), 
			printf ("%lld\n", query(c, d) - query(a - 1, d) - query(c, b - 1) + query(a - 1, b - 1));
		else a = Read(), b = Read(), c = Read(), modify(a, b, c);
	}
	return 0;
}

单点修改区间查询

LOJ #135. 二维树状数组 3:区间修改,区间查询

一般区间修改的树状数组维护的都是差分数组,那么二维的也应该维护二维差分数组,接着是求和:

\[\begin{aligned}&\sum_{x=1}^{a}\sum_{y=1}^{b}\sum_{i=1}^{x}\sum_{j=1}^{y}t_{i,j}\\ =&\sum_{i=1}^{a}\sum_{j=1}^{b}(a-i+1)(b-j+1)t_{i,j}\\ =&(a+1)(b+1)\sum_{i=1}^{a}\sum_{j=1}^{b}t_{i,j}-\\ &(b+1)\sum_{i=1}^{a}\sum_{j=1}^{b}t_{i,j}\cdot i-\\ &(a+1)\sum_{i=1}^{a}\sum_{j=1}^{b}t_{i,j}\cdot j+\\ &\sum_{i=1}^{a}\sum_{j=1}^{b}t_{i,j}\cdot ij\end{aligned}\]

const int N = 4106;

inline ll Read()
{
	ll x = 0, f = 1;
	char c = getchar();
	while (c != '-' && (c < '0' || c > '9')) c = getchar();
	if (c == '-') f = -f, c = getchar();
	while (c >= '0' && c <= '9') x = (x << 3) + (x << 1) + c - '0', c = getchar();
	return x * f;
}

int n, m;
ll t[4][N][N];

void modify(int x, int y, ll val)
{
	for (int i = x; i <= n; i += i & -i)
		for (int j = y; j <= m; j += j & -j)
			t[0][i][j] += val,
			t[1][i][j] += val * x,
			t[2][i][j] += val * y,
			t[3][i][j] += val * x * y;
}

ll query(int x, int y)
{
	ll ans = 0;
	for (int i = x; i; i -= i & -i)
		for (int j = y; j; j -= j & -j)
			ans += (x + 1) * (y + 1) * t[0][i][j] - 
			       (y + 1) * t[1][i][j] -
				   (x + 1) * t[2][i][j] + 
				   t[3][i][j];
	return ans;
}

int main()
{
	n = Read(), m = Read(); 
	for (int op, a, b, c, d, k; scanf ("%d", &op) != EOF; )
	{
		if(op == 2) a = Read(), b = Read(), c = Read(), d = Read(), 
			printf ("%lld\n", query(c, d) - query(a - 1, d) - query(c, b - 1) + query(a - 1, b - 1));
		else a = Read(), b = Read(), c = Read(), d = Read(), k = Read(), 
			modify(a, b, k), modify(a, d + 1, -k), modify(c + 1, b, -k), modify(c + 1, d + 1, k);
	}
	return 0;
}

标签:树状,sum,二维,数组,基本操作,ll,getchar
来源: https://www.cnblogs.com/GJY-JURUO/p/14876365.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有