結果

問題 No.843 Triple Primes
ユーザー claw88claw88
提出日時 2019-06-28 21:37:07
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,261 ms / 2,000 ms
コード長 1,651 bytes
コンパイル時間 1,473 ms
コンパイル使用メモリ 168,288 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-26 07:26:29
合計ジャッジ時間 25,711 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define mod 1000000007
#define dom 998244353
#define all(c) begin(c),end(c)
#define rep(i,n) repi(i,0,n)
#define repi(i,a,b) for(int i=(a);i<(b);i++)
#define rrep(i,a,b) for(int i=(b)-1;i>=(a);i--)
template <typename T> bool chmax(T &a, const T &b) { if (a < b) { a = b; return true; } return false; }
template <typename T> bool chmin(T &a, const T &b) { if (b < a) { a = b; return true; } return false; }
#define long long long
#define List vector
#define var auto
#define Count size()
#define Length size()
int dd[] = { 0, 1, 0, -1, 0 }; //→↓←↑

#define MAX 500000
vector<int> prime;
bool is_prime[MAX + 1];
void eratos(int n)
{
	rep(i, n + 1) is_prime[i] = true;
	is_prime[0] = is_prime[1] = false;
	rep(i, n + 1)
	{
		if (is_prime[i])
		{
			prime.emplace_back(i);
			for (int j = 2 * i; j <= n; j += i)
				is_prime[j] = false;
		}
	}
}

// オイラーのファイ関数
// 自然数nに対し、1からnまでの自然数の中で
// nと互いに素なものの数 
//int P[MAX + 1];
//void phi(int n)
//{
//	for (int i = 0; i <= n; i++)
//		P[i] = i;
//	for (int i = 2; i <= n; i++)
//	{
//		if (is_prime[i])
//		{
//			P[i] -= P[i] / i;
//			for (int j = 2 * i; j <= n; j += i)
//				P[j] -= P[j] / i;
//		}
//	}
//}
void solve()
{
	int N; cin >> N;
	eratos(N);
	int cnt = 0;
	for (auto r : prime)
	{
		long L = (long)r * r;
		for (auto p : prime)
		{
			if (p > L / 2)
			{
				break;
			}
			if (L - p <= N && is_prime[L - p])
			{
				cnt++;
				if (L - p != p)cnt++;
			}
		}
	}
	cout << cnt << endl;
}

int main()
{
	cin.tie(0);
	ios::sync_with_stdio(false);
	solve();
	return 0;
}
0