https://nanti.jisuanke.com/t/20689
在一个星光摧残的夜晚，蒜头君一颗一颗的数这天上的星星。
蒜头君给在天上巧妙的画了一个直角坐标系，让所有的星星都分布在第一象。天上有 nn 颗星星，他能知道每一颗星星的坐标和亮度。
现在，蒜头君问自己 qq 次，每次他问自己每个矩形区域的星星的亮度和是多少（包含边界上的星星）。
输入格式
第一行输入一个整数 n(1 le n le 50000)n(1≤n≤50000) 表示星星的数量。
接下里 nn 行，每行输入三个整数 x,y,w(0 le x, y, wle 2000)x,y,w(0≤x,y,w≤2000)，表示在坐标 (x,y)(x,y) 有一颗亮度为 ww 的星星。注意一个点可能有多个星星。
接下来一行输入一个整数 q(1 le q le 50000)q(1≤q≤50000)，表示查询的次数。
接下来 qq 行，每行输入四个整数 x_1, y_1, x_2, y_2x1,y1,x2,y2，其中 (x_1, y_1)(x1,y1) 表示查询的矩形的左下角的坐标，(x_2, y_2)(x2,y2) 表示查询的矩形的右上角的坐标，0 le x_1 le x_2 le 20000≤x1≤x2≤2000，0 le y_1 le y_2 le 20000≤y1≤y2≤2000。
输出格式
对于每一次查询，输出一行一个整数，表示查询的矩形区域内的星星的亮度总和。
样例输入
5
5 0 6
7 9 7
8 6 13
9 7 1
3 0 19
4
0 8 7 9
0 0 7 10
2 7 10 9
5 4 7 5
样例输出
7
32
8
0
暴力超时，还是dp吧：
#include <bits/stdc++.h>
using namespace std;
int mp[2005][2005];
int dp[2005][2005];
int visit[2005][2005];
int mydp(int i, int j){
// if(dp[i][j])
if(visit[i][j] != 0)
return dp[i][j];
visit[i][j] = 1;
if(i < 0 || j < 0)
return 0;
if(i == 0 && j == 0)
return dp[i][j] = mp[i][j];
if(i == 0)
return dp[i][j] = mp[i][j] + mydp(i, j - 1);
if(j == 0)
return dp[i][j] = mp[i][j] + mydp(i - 1, j);
return dp[i][j] = mp[i][j] + mydp(i - 1, j) + mydp(i, j - 1) - mydp(i - 1, j - 1);
}
int main(){
int n, x, y, w;
cin >> n;
for(int i = 0; i < n; i++){
cin >> x >> y >> w;
mp[x][y] += w;
}
int m, x1, x2, y1, y2;
cin >> m;
mydp(2000, 2000);
for(int i = 0; i < m; i++){
cin >> x1 >> y1 >> x2 >> y2;
// cout << "x2, y2: " << x2 << " " << y2 << "; " << dp[x2][y2] << endl;
cout << dp[x2][y2] - dp[x1 - 1][y2] - dp[x2][y1 - 1] + dp[x1 - 1][y1 - 1] << endl;
}
return 0;
}