一种新的二维差分

#include<bits/stdc++.h>
#define int long long
typedef long long ll;
using namespace std;
const int N = 2e3 + 10;
const int mod  = 1e9 + 7;
void solve(){
    int n,m;
    cin>>n>>m;
    int a[N][N];
    memset(a,0,sizeof a);
    while (m--){
        int x1,y1,x2,y2;
        cin>>x1>>y1>>x2>>y2;
        a[x1][y1] ++;
        a[x1][y2+1] --;
        a[x2+1][y1] --;
        a[x2+1][y2+1]++;
    }
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= n; ++j) {
            a[i][j] += a[i-1][j] + a[i][j-1] - a[i-1][j-1];
        }
    }
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= n; ++j) {
            if(a[i][j] % 2 == 0){
                cout<<"0";
            } else{
                cout<<"1";
            }
        }
        cout<<endl;
    }
}
signed main(){
    ios::sync_with_stdio(false);
    cin.tie(NULL),cout.tie(NULL);
    int t = 1;
    //cin>>t;
    while(t--){
        solve();
    }
    return 0 ;
}