Codeforces 835C - Star sky 【预处理二维前缀和】

本文介绍了一个星空亮度模拟查询的问题及解决思路。在一个100*100的星空上,有n颗星星亮度随时间周期性变化。文章通过预处理前缀和的方式,实现了快速查询在特定时间点内,任意给定矩形区域内的星星总亮度。
C. Star sky
time limit per test
2 seconds
memory limit per test
256 megabytes
input
standard input
output
standard output

The Cartesian coordinate system is set in the sky. There you can see n stars, the i-th has coordinates (xiyi), a maximum brightness c, equal for all stars, and an initial brightness si (0 ≤ si ≤ c).

Over time the stars twinkle. At moment 0 the i-th star has brightness si. Let at moment t some star has brightness x. Then at moment (t + 1) this star will have brightness x + 1, if x + 1 ≤ c, and 0, otherwise.

You want to look at the sky q times. In the i-th time you will look at the moment ti and you will see a rectangle with sides parallel to the coordinate axes, the lower left corner has coordinates (x1iy1i) and the upper right — (x2iy2i). For each view, you want to know the total brightness of the stars lying in the viewed rectangle.

A star lies in a rectangle if it lies on its border or lies strictly inside it.

Input

The first line contains three integers nqc (1 ≤ n, q ≤ 1051 ≤ c ≤ 10) — the number of the stars, the number of the views and the maximum brightness of the stars.

The next n lines contain the stars description. The i-th from these lines contains three integers xiyisi (1 ≤ xi, yi ≤ 1000 ≤ si ≤ c ≤ 10) — the coordinates of i-th star and its initial brightness.

The next q lines contain the views description. The i-th from these lines contains five integers tix1iy1ix2iy2i (0 ≤ ti ≤ 1091 ≤ x1i < x2i ≤ 1001 ≤ y1i < y2i ≤ 100) — the moment of the i-th view and the coordinates of the viewed rectangle.

Output

For each view print the total brightness of the viewed stars.

Examples
input
Copy
2 3 3
1 1 1
3 2 0
2 1 1 2 2
0 2 1 4 5
5 1 1 5 5
output
Copy
3
0
3
input
Copy
3 4 5
1 1 2
2 3 0
3 3 1
0 1 1 100 100
1 2 2 4 4
2 2 1 4 7
1 50 50 51 51
output
Copy
3
3
5
0
Note

Let's consider the first example.

At the first view, you can see only the first star. At moment 2 its brightness is 3, so the answer is 3.

At the second view, you can see only the second star. At moment 0 its brightness is 0, so the answer is 0.

At the third view, you can see both stars. At moment 5 brightness of the first is 2, and brightness of the second is 1, so the answer is 3.

【题意】

有一片100*100的星空,上面有n颗星星,每个星星有一个亮度,且在0~C范围内周期性变化,现在给出q个查询,每个查询给出时间和一个矩形,求在该时间时矩形内星星的亮度和。


【思路】:

  学树状数组的时候遇到的,一开始以为是个二维树状数组题,结果发现不需要更新数组,直接预处理前缀和就好了。

我们用sum[i][j][k]来表示左下角(0,0),右上角为(i,j),亮度为k的矩形星星总数,迭代维护sum[i][j][k]的数量,有点Dp的意思,状态转移方程

sum[i][j][k] += sum[i-1][j][k] + sum[i][j-1][k] - sum[i-1][j-1][k],方程中用到的容斥的思想,画图就很好理解了。

对于q个询问,求出每个矩形星星的亮度总和,再用系数(t+i)%(c+1)将时间转化到t时刻,累和得出结果


#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
const int Maxn = 1e5 + 6;
int n,q,c;
int sum[109][109][19];
using namespace std;
void init(){
	for (int i = 1; i <= 100; i++)
		for (int j = 1; j <=  100; j++)
			for (int k = 0; k <= 10; k++){
				sum[i][j][k] += sum[i-1][j][k] + sum[i][j-1][k] - sum[i-1][j-1][k];
			}
}
int main ()
{
	scanf("%d%d%d",&n,&q,&c);
	for (int i = 1; i <= n; i++) {
		int x,y,t;
		scanf("%d%d%d",&x,&y,&t);
		sum[x][y][t]++;
	}
	init();
	for (int i = 1; i <= q; i++){
		int cnt = 0;
		int t,x1,y1,x2,y2;
		scanf("%d%d%d%d%d",&t,&x1,&y1,&x2,&y2);
		for (int i = 0; i <= c; i++){
			int k = (t+i)%(c+1);
			int p = sum[x2][y2][i]+sum[x1-1][y1-1][i]-sum[x1-1][y2][i]-sum[x2][y1-1][i];
			cnt += k*p;
		}
		printf("%d\n",cnt);
	}
}





考虑柔性负荷的综合能源系统低碳经济优化调度【考虑碳交易机制】(Matlab代码实现)内容概要:本文围绕“考虑柔性负荷的综合能源系统低碳经济优化调度”展开,重点研究在碳交易机制下如何实现综合能源系统的低碳化与经济性协同优化。通过构建包含风电、光伏、储能、柔性负荷等多种能源形式的系统模型,结合碳交易成本与能源调度成本,提出优化调度策略,以降低碳排放并提升系统运行经济性。文中采用Matlab进行仿真代码实现,验证了所提模型在平衡能源供需、平抑可再生能源波动、引导柔性负荷参与调度等方面的有效性,为低碳能源系统的设计与运行提供了技术支撑。; 适合人群:具备一定电力系统、能源系统背景,熟悉Matlab编程,从事能源优化、低碳调度、综合能源系统等相关领域研究的研究生、科研人员及工程技术人员。; 使用场景及目标:①研究碳交易机制对综合能源系统调度决策的影响;②实现柔性负荷在削峰填谷、促进可再生能源消纳中的作用;③掌握基于Matlab的能源系统建模与优化求解方法;④为实际综合能源项目提供低碳经济调度方案参考。; 阅读建议:建议读者结合Matlab代码深入理解模型构建与求解过程,重点关注目标函数设计、约束条件设置及碳交易成本的量化方式,可进一步扩展至多能互补、需求响应等场景进行二次开发与仿真验证。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值