1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
#include <bits/stdc++.h>
#define all(x) (x).begin(),(x).end()
using namespace std;

using ll = long long;
using ld = long double;

//#define int ll
#define sz(x) ((int)(x).size())

using pii = pair<int,int>;
using tii = tuple<int,int,int>;

const int nmax = 5e3 + 5;
int freq[nmax * nmax];

signed main() {
   cin.tie(0) -> sync_with_stdio(0);
   int n;
   cin >> n;
   
   for(int a = 1; a <= n; a++) {
      for(int b = a; a * a + b * b <= n * n; b++) {
         freq[a * a + b * b]++;
      }
   }
   ll rez = 0;
   for(int d = 1; d <= n; d++) {
      for(int c = 1; d * d - c * c > 0; c++) {
         rez += freq[d * d - c * c];
      }
   }
   cout << rez << '\n';
   
}

/**
      Binecuvanteaza Doamne Ukraina.
--
*/