結果

問題 No.843 Triple Primes
コンテスト
ユーザー kenta255
提出日時 2019-12-11 01:08:40
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,374 bytes
コンパイル時間 1,980 ms
コンパイル使用メモリ 205,524 KB
最終ジャッジ日時 2025-01-08 10:30:10
ジャッジサーバーID
（参考情報） 		judge1 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 TLE * 1
other AC * 16 TLE * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define P pair<ll,ll>
#define FOR(I,A,B) for(ll I = ll(A); I < ll(B); ++I)
#define FORR(I,A,B) for(ll I = ll((B)-1); I >= ll(A); --I)
#define TO(x,t,f) ((x)?(t):(f))
#define SORT(x) (sort(x.begin(),x.end())) // 0 2 2 3 4 5 8 9
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin()) //xi>=v  x is sorted
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin()) //xi>v  x is sorted
#define NUM(x,v) (POSU(x,v)-POSL(x,v))  //x is sorted
#define REV(x) (reverse(x.begin(),x.end())) //reverse
ll gcd(ll a,ll b){if(a%b==0)return b;return gcd(b,a%b);}
ll lcm(ll a,ll b){ll c=gcd(a,b);return ((a/c)*(b/c)*c);}
#define NEXTP(x) next_permutation(x.begin(),x.end())
const ll INF=ll(1e18)+ll(7);
const ll MOD=1000000007LL;
#define out(a) cout<<fixed<<setprecision((a))


vector<int> make_p(int x){
	vector<bool> isp(x+1,true);
	isp[0]=isp[1]=false;
	for(int i=2;i*i<=x;i++){
		if(isp[i]){
			for(int j=2*i;j<=x;j+=i){
				isp[j] = false;
			}
		}
	}
	vector<int> res;
	for(int i=2;i<=x;i++){
		if(isp[i])res.push_back(i);
	}
	return res;
}


int main(){
	int N;
	cin >> N;
	auto p = make_p(N);
	vector<ll> p2;
	set<ll> sep2;
	for(auto x:p){
		ll a = x;
		p2.push_back(a*a);
		sep2.insert(a*a);
	}
	int ans = 0;
	for(auto x:p){
		for(auto y:p){
			if(sep2.count(x+y))ans++;
		}
	}
	cout << ans << endl;
}
0