Задача "Суммы на прямоугольниках"

#include <iostream>
#include <vector>
#include <algorithm>
 
using namespace std;
 
typedef long long ll;
 
ll AT(vector<vector<ll>> p, ll i, ll j) {
    if (i < 0 || j < 0) {
        return 0ll;
    }
    return p[i][j];
}
 
int main() {
    ll n, m, k;
    cin >> n >> m >> k;
    vector<vector<ll>> a(n, vector<ll>(m)), pref(n, vector<ll>(m));
    for (ll i = 0; i < n; i++) {
        for (ll j = 0; j < m; j++) {
            cin >> a[i][j];
        }
    }
    pref[0][0] = a[0][0];
    for (ll i = 0; i < n; i++) {
        for (ll j = 0; j < m; j++) {
            if (i >= 1 && j >= 1) {
                pref[i][j] += pref[i - 1][j] + pref[i][j - 1] + a[i][j] - pref[i-1][j-1];
            }
            if (i >= 1 && j < 1) {
                pref[i][j] += pref[i - 1][j] + a[i][j];
            }
            if (i < 1 && j >= 1) {
                pref[i][j] += pref[i][j - 1] + a[i][j];
            }
        }
    }
    for (ll i = 0; i < k; i++) {
        ll x1, x2, y1, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        x1--;
        y1--;
        x2--;
        y2--;
        cout << AT(pref, x2, y2) - AT(pref, x1 - 1, y2) - AT(pref, x2, y1 - 1) + AT(pref, x1 - 1, y1 - 1) << endl;
    }
}

#include <iostream>
#include <vector>
 
int main()
{
    std::ios_base::sync_with_stdio(false);
    std::cin.tie(0);
 
    size_t n {}, m {}, k {};
    std::cin >> n >> m >> k;
    ++n;
    ++m;
 
    std::vector<std::vector<uint64_t>> mx(n, std::vector<uint64_t>(m, 0));
 
    for(size_t y { 1 }; y < n; ++y)
        for(size_t x { 1 }; x < m; ++x)
        {
            uint64_t val {};
            std::cin >> val;
            mx[y][x] = val + mx[y - 1][x] + mx[y][x - 1] - mx[y - 1][x - 1];
        }
 
    size_t a {}, b {}, c {}, d {};
    for(size_t i { 0 }; i < k; ++i)
    {
        std::cin >> b >> a >> d >> c;
        std::cout << mx[d][c] + mx[b - 1][a - 1] - mx[d][a - 1] - mx[b - 1][c] << "\n";
    }
 
    return 0;
}