#include <math.h>
#include <windows.h>
#include <iostream>
typedef unsigned int		u32;
typedef unsigned __int64	u64;

#define N 10000
#define MIN(a,b) ((a)<(b)?(a):(b))
#define MAX(a,b) ((a)>(b)?(a):(b))
using namespace std;

u64 GCD(u64 x, u64 y) {
	return x ? GCD(y % x, x) : y;
}

u64 Sqrt(u64 n) {
	return (u64)( sqrt((double)n + 0.000000000001) );
}

int main(int argc, char* argv[])
{
	// スレッド優先度を下げる.
	SetThreadPriority(GetCurrentThread(), THREAD_PRIORITY_LOWEST);

	// 面積は r^2 * (b^2-a^2) / (2*(b^2+a^2)).
	// この分母をキャンセルする r をここで先に求めておく.
	// q は k=(a^2+b^2-1)/4 をインデックスにした配列で,
	// r が q[k] の倍数のとき分母がキャンセルされる.
	u32 * q = new u32 [N*N/4];
	for (u32 i=0; i<N*N/4; i++) q[i] = 1;
	for (u64 i=5; i<N*N; i+=4) {
		if (q[(u32)((i-1)/4)] != 1) continue;
		for (u64 j = i; j<N*N; j *= i*i) {
			for (u64 k = j; k<N*N; k += j) {
				if (k%4 != 1) continue;
				q[(u32)((k-1)/4)] *= (u32)i;
			}
		}
	}

	u64 s = 0;
	// gcd(a,b)=1 となる a,b を総当り.
	for (u64 b=1; b<N; b++) {
		for (u64 a=1; a<b && a*a+b*b<N*N; a++) {
			if (GCD(a, b) != 1) continue;
			u64 n = b*b - a*a;
			u64 d = b*b + a*a;
			u64 g = (a%2==1&&b%2==1) ? 2 : 1;
			n /= g;
			d /= g;
			d = q[(d-1)/4];
			if (n % 2 == 0) n /= 2;
			else d *= 2;
			u64 h = a*a+b*b;
			for (u64 i=1; h*i*i<N*N; i++) {
				u64 r = Sqrt(i*i*h-1);
				u64 s0 = r / d;
				s += s0;
			}
		}
	}
	cout << s << endl;

	return 0;
}