#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 #include #include using namespace std; using ll = long long; using P = pair; constexpr int INF = 1001001001; // constexpr int mod = 1000000007; constexpr int mod = 998244353; template inline bool chmax(T& x, T y){ if(x < y){ x = y; return true; } return false; } template inline bool chmin(T& x, T y){ if(x > y){ x = y; return true; } return false; } struct ExtendedSieve{ vector table; const int n; ExtendedSieve(int n) : n(n) { table.assign(n + 1, 0); table[0] = table[1] = 1; for(int i = 2; i <= n; ++i){ if(table[i] > 0) continue; table[i] = i; for(int j = i + i; j <= n; j += i){ if(table[j] == 0) table[j] = i; } } } int operator[](const int& k){ return table[k]; } bool is_prime(int x){ if(x < 2) return false; return x == table[x]; } map prime_factorization(int x){ map res; while(x > 1){ ++res[table[x]]; x /= table[x]; } return res; } }; int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int N; cin >> N; ExtendedSieve sieve(N); vector prime, prime2; for(int i = 2; i <= N; ++i){ if(!sieve.is_prime(i)) continue; prime.emplace_back(i); if((ll)i * i <= N * 2) prime2.emplace_back(i * i); } ll ans = 0; for(auto p : prime){ for(auto s : prime2){ ans += upper_bound(begin(prime), end(prime), s - p) - lower_bound(begin(prime), end(prime), s - p); } } cout << ans << endl; return 0; }