#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; const ll PMAX = 300000; //unordered_setを使って高速化できるのだろうか //O(π(N)×μ(N))≒10^8 //(μ(N)は1~Nまでの素数2つの和として表せる数の個数とする) int main(void) { int n; vectorprime; vectorflag(PMAX); cin >> n; //primeにn以下のすべての素数が入る //flag[i-2]はiが素数の時false for (int x = 2; x <= PMAX; ++x) { if (flag[x - 2])continue; if (x <= n)prime.push_back(x); for (int y = x + x; y <= PMAX; y += x) { flag[y - 2] = true; } } vectorcrush(n + n + 1); unordered_setsum; for (int a : prime) { for (int b : prime) { sum.insert(a + b); } } //crush[x]:=N以下の素数pのうち、x+pも素数である個数。 //xが2つの素数の和で表せる場合のみを算出 for (int x : sum) { for (int p : prime) { if (!flag[p + x - 2]) { crush[x]++; } } } ll ans = 0; //a,bを総当たり for (int a : prime) { for (int b : prime) { //cとして考えられる数を総当たり ans += crush[a + b]; } } //2つの素数を含む組みを消す for (int a : prime) { ans -= crush[a + a] * 3; } cout << ans / 6 << endl; return 0; }