Lưu ý: Các code mẫu dưới đây chỉ mang tính tham khảo và có thể không AC được bài tập này

Code mẫu của RR

#include <iostream>
#include <algorithm>
#define FOR(i,a,b) for(int i=a; i<=b; i++)
#define BASE 2009
#define ll long long
using namespace std;

ll X,Y,res;

void solve() {
    //x=0, y!=0
    FOR(y,-min(X,Y),min(X,Y))
    if (y!=0) {
        int u=y; if (y<0) u=-u;
        res=(res+(X-u+1)*(Y-u+1))%BASE;
    }

    //x!=0, y=0
    FOR(x,-min(X,Y),min(X,Y))
    if (x!=0)
        res=(res+(X-abs(x)+1)*(Y-abs(x)+1))%BASE;
}

void cal() {
    //x>0, y>0
    //x>y
    FOR(x,2,min(X,Y)) {
        int y=min(Y-x,(ll)x-1);
        if (y>0) {
            res+=2*(X-x+1)*(y*Y-y*x-(((y+1)*y)>>1)+y);
            res%=BASE;
        }
    }

    //x<y
    FOR(y,2,min(X,Y)) {
        int x=min(Y-y,(ll)y-1);
        if (x>0) {
            res+=2*(X-y+1)*(x*Y-x*y-(((x+1)*x)>>1)+x);
            res%=BASE;
        }
    }

    //x=y
    FOR(x,1,min(X,Y))
    if (Y>=x+x) {
        res+=2*(X-x+1)*(Y-x-x+1);
        res%=BASE;
    }
}

int main() {
    int test; scanf("%d",&test);
    while (test--) {
        res=0;
        scanf("%d %d",&X,&Y);
        cal();
        swap(X,Y);
        cal();
        solve();
        cout<<res<<endl;
    }
    return 0;
}

Code mẫu của hieult

#include <cstdio>
//#include <conio.h>
#define min(a,b) (a<b?a:b)
#define max(a,b) (a>b?a:b)

const long long du = 2009;
long long test,x,y;

long long chay(long long a,long long b,long long c){
    return (c*(c+1)*(2*c+1)/6+a*b*(c+1)-(c+1)*c/2*(a+b))%du;
}

long long tinh(long long x,long long y)
{
     long long kq = 0;
     for(int i = 1;i<=x&&i<=y;++i){
          if(2*i<=y){
                kq=(kq+(y-i+1+y-i-i+1)*(1+i)/2%du*(x-i+1))%du;
                int k=min(y-i+1,x);
                kq=(kq+chay(x+1,y-i+1,k)-chay(x+1,y-i+1,i))%du;
            }else kq=(kq+(y-i+1)*(y-i+1+1)/2%du*(x-i+1))%du;
     }
     return kq;
}

int main()
{
    scanf("%lld",&test);
    while(test--){
         scanf("%lld %lld",&x,&y);
         long long KQ = (tinh(x,y)+tinh(y,x))%du;
         printf("%lld\n",(KQ*2)%du);
    }
   // getch();
}