Star sky CodeForces - 835C

用一个三维数组cnt[x][y][k]表示从(1, 1)到(x, y)亮度为k的个数，然后查询的时候就是对于每一个亮度，计算出这个亮度t秒后的亮度和当前这个亮度的个数，答案就是他们的乘积， 然后遍历每种亮度就可以算出答案

#include<map>
#include<set>
#include<ctime>
#include<cmath>
#include<stack>
#include<queue>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define lowbit(x) (x & (-x))

typedef unsigned long long int ull;
typedef long long int ll;
const double pi = 4.0*atan(1.0);
const int inf = 0x3f3f3f3f;
const int maxn = ;
const int maxm = ;
const int mod = 1e9+;
using namespace std;

int n, m;
int T, tol;
int cnt[maxn][maxn][maxn/];

void init() {
    memset(cnt, , sizeof cnt);
}

int main() {
    init();
    int q, c;
    scanf("%d%d%d", &n, &q, &c);
    int mx=, my=;
    int x, y, s;
    for(int i=; i<=n; i++) {
        scanf("%d%d%d", &x, &y, &s);
        cnt[x][y][s]++;
        mx = max(mx, x);
        my = max(my, y);
    }
    for(int i=; i<=; i++) {
        for(int j=; j<=; j++) {
            for(int k=; k<=c; k++) {
                cnt[i][j][k] = cnt[i][j][k] + cnt[i-][j][k] + cnt[i][j-][k] - cnt[i-][j-][k];
            }
        }
    }
    int t, x1, y1, x2, y2;
    while(q--) {
        scanf("%d%d%d%d%d", &t, &x1, &y1, &x2, &y2);
        int ans = ;
        for(int k=; k<=c; k++) {
            int tmp = (k+t) % (c+);
            int cx = cnt[x2][y2][k] - cnt[x1-][y2][k] - cnt[x2][y1-][k] + cnt[x1-][y1-][k];
            ans += tmp * cx;
        }
        printf("%d\n", ans);
    }
    return ;
}