1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
#include<bits/stdc++.h>
using namespace std;
const int N = 5001;
int n,DP[N*N],sq[N*N];
bitset<N*N>p2;
int main(){
    ios_base::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);

    cin >> n;
    long long res = 0;
    for(int i = 1;i <= n;i++)p2[i*i] = 1,sq[i*i] = i;
    for(int i = 1;i <= n;i++)
        for(int j = i;i*i + j*j <= n*n;j++)
           DP[i*i + j*j]++;
    for(int i = 1;i <= n;i++){
        for(int j = i-1;j > 0;j--){
            res += DP[i*i - j*j];
        }
    }
    cout << res << "\n";
}