结果填空：方阵

 

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
int main()
{
    int n=1000;
    int k=500;
    ll sum=2*(n-1)*n%mod;
    for(int i=1; i<n; i++)
    {
        for(int j=1; j<n; j++)
        {
            if(__gcd(i,j)==1&&i*i+j*j<=k*k)
            {
                sum+=2*(n-i)*(n-j);
                sum%=mod;
            }
        }
    }
    printf("%lld\n",sum);
    return 0;
}