704 二维差分

// 704 二维差分.cpp : 此文件包含 "main" 函数。程序执行将在此处开始并结束。
//

/*
http://oj.daimayuan.top/course/22/problem/895


给一个 n×m 的矩阵 a1,1,a1,2,…,a1,m,…,an,m和 q个修改操作。
一开始所有位置都是 0，每次修改给出五个数 x1,y1,x2,y2,d，令所有 ai,j(x1≤i≤x2,y1≤j≤y2) 加上 d。

求出所有操作结束之后，整个矩阵的值。

输入格式
第一行六个整数 n,m,q,A,B,C(1≤n,m≤2×103,1≤q≤106,0≤A,B,C<232)。

为了防止输入过大，询问用下面代码生成。
输出格式
为了防止输出过大，输出操作完之后矩阵中所有数的异或和。

样例输入
5 5 10 1 2 3
样例输出
3809069637
*/
#define  _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std;


const int N = 2010;
int n, m, q;
unsigned int A, B, C;
long long  s[N][N];

inline unsigned int rng61() {
    A ^= A << 16;
    A ^= A >> 5;
    A ^= A << 1;
    unsigned int t = A;
    A = B;
    B = C;
    C ^= t ^ A;
    return C;
}



int main() {
    scanf("%d%d%d%u%u%u", &n, &m, &q, &A, &B, &C);
    for (int i = 1; i <= q; i++) {
        int x1 = rng61() % n + 1, x2 = rng61() % n + 1;
        int y1 = rng61() % m + 1, y2 = rng61() % m + 1;
        unsigned int d = rng61();
        if (x1 > x2) swap(x1, x2);
        if (y1 > y2) swap(y1, y2);
        // 执行修改操作
        s[x1][y1] += d;
        s[x2 + 1][y2 + 1] += d;
        s[x1][y2 + 1] -= d;
        s[x2 + 1][y1] -= d;
    }
    long long ans = 0;
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= m; j++) {
            s[i][j] += s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1];
            ans ^= s[i][j];
        }

    cout << ans << endl;


}