#define _CRT_SECURE_NO_WARNINGS #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define aLL(g) (g).begin(),(g).end() #define REP(i, x, n) for(int i = x; i < n; i++) #define rep(i,n) REP(i,0,n) #define F(i,j,k) fill(i[0],i[0]+j*j,k) #define P(p) cout<<(p)< vi; typedef vector vl; typedef vector vd; typedef pair pii; typedef pair pll; typedef long long ll; int dx[] = { 0, 1, 0, -1 }; int dy[] = { -1, 0, 1, 0 }; int sttoi(std::string str) { int ret; std::stringstream ss; ss << str; ss >> ret; return ret; } bool sort_greater(const pair &a, const pair &b) { return a.second > b.second; } int isPrime(int a, int b){ if (b == 0)return a; return isPrime(b,a%b); } void solve() { int L; cin >> L; int l = L / 4; set deg; int count = 0; for (int i = 0; i*i < l; i++){ for (int j = 1; j < i; j++){ if (isPrime(i, j) != 1)continue; if ((i - j) % 2 == 0)continue; int a = i*i-j*j;int b = 2*i*j;int c = i*i+j*j; if (a + b + c <= l){ deg.insert((double)a / (double)c); } } } count = deg.size(); count %= 1000003; P(count); } int main() { solve(); return 0; }