codeforces-835C. Star sky(维护二维前缀和，dp)

 传送门

在不大于100*100的地图上有星星，每个星星的亮度随时间变化，它们初始亮度不同，但能达到的最高亮度相同，为C。初始亮度为p，t时刻亮度为（p+t)%(C+1)

输入

n q c（为别代表星星数目，查询次数，最高亮度）

n行

x y s(星星的横纵坐标和初始亮度)

q行

t x1 y1 x2 y2(查询t时刻矩阵内及边界上星星的总亮度，其中(x1,y1)是矩阵左上角, (x2,y2)是矩阵右下角）

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define INF 0x3f3f3f3f
#define MOD 1000000007
using namespace std;
typedef long long LL;

const int maxc = 15;
const int maxn = 110;
int cnt[maxc][maxn][maxn];

int N, Q, C;

void init() {
    for (int k = 0; k <= C; k++) {
        for (int i = 1; i <= 100; i++) {
            for (int j = 1; j <= 100; j++) {
                cnt[k][i][j] += cnt[k][i - 1][j] + cnt[k][i][j - 1] - cnt[k][i - 1][j - 1];
            }
        }
    }
}

int solve(int t, int x1, int y1, int x2, int y2) {
    int res = 0;
    for (int k = 0; k <= C; k++) {
        int tmp = cnt[k][x2][y2] + cnt[k][x1 - 1][y1 - 1] - cnt[k][x1 - 1][y2] - cnt[k][x2][y1 - 1];
        res = res + ((k + t) % (C + 1)) * tmp;
    }
    return res;
}

int main() {
    scanf("%d%d%d", &N, &Q, &C);
    int x, y, s;
    for (int i = 1; i <= N; i++) {
        scanf("%d%d%d", &x, &y, &s);
        cnt[s][x][y]++;
    }
    init();
    int t, x1, x2, y1, y2;
    for (int i = 1; i <= Q; i++) {
        scanf("%d%d%d%d%d", &t, &x1, &y1, &x2, &y2);
        printf("%d
", solve(t, x1, y1, x2, y2));
    }
    return 0;
}