方阵

广场上的小朋友们排成了整齐的方阵。具体来说，我们可以把每个小朋友看做是一个点，那么小朋友们就形成了 n\times nn×n 的点阵。方阵中，小朋友 AA 和小朋友 BB 互相可以看见，当且仅当二人之间的连线不经过别的小朋友，且他们之间的距离不超过 kk （因为太远就看不见了）。我们想知道有多少对小朋友互相可以看见。(A,B)(A,B) 与 (B,A)(B,A) 算同一对。

例如， n=2,k=1n=2,k=1 时答案为 44， n=2,k=2n=2,k=2 时答案为 66（距离为 11 的有 44 对，距离为 \sqrt 22​ 的有 22 对）， n=3,k=2n=3,k=2 时答案为 2020 。

现在我们想要知道，当 n=1000,k=500n=1000,k=500 时的答案是多少。由于答案过大，请回答对 10^9+7109+7 取模后的结果。

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

 

本来感觉哇！！！好神奇为什么这么写，结果跟着他运行一遍恍然大悟 

• 判断条件为gcd(x, y) == 1和x*x+y*y<=k*k，其中x代表两点之间横坐标之差，y表示纵坐标之差。
• gcd(x, y) == 1  ，当且仅当这个时候小孩子才能互相看见，x*x+y*y<=k*k  利用勾股定理，只有在给出的条件下才能互相看见