#476625 (C++17) No.843 Triple Primes

• 問題: No.843 Triple Primes
• ユーザー: fura
• 提出日時: 2020-05-05 16:37:29
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: WA
• 実行時間: -
• コード長: 796 bytes
• コンパイル時間: 2,698 ms
• コンパイル使用メモリ: 202,892 KB
• 実行使用メモリ: 4,384 KB
• 最終ジャッジ日時: 2023-09-09 03:55:47
• 合計ジャッジ時間: 25,771 ms
• ジャッジサーバーID （参考情報）: judge12 / judge13

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 4,380 KB
testcase_01	AC	1,147 ms; 4,380 KB
testcase_02	AC	2 ms; 4,376 KB
testcase_03	AC	2 ms; 4,380 KB
testcase_04	AC	2 ms; 4,376 KB
testcase_05	AC	2 ms; 4,380 KB
testcase_06	AC	3 ms; 4,380 KB
testcase_07	AC	770 ms; 4,376 KB
testcase_08	AC	905 ms; 4,380 KB
testcase_09	AC	1,133 ms; 4,376 KB
testcase_10	AC	822 ms; 4,380 KB
testcase_11	AC	979 ms; 4,380 KB
testcase_12	AC	1,084 ms; 4,380 KB
testcase_13	AC	1,034 ms; 4,376 KB
testcase_14	AC	1,039 ms; 4,376 KB
testcase_15	AC	792 ms; 4,380 KB
testcase_16	AC	813 ms; 4,380 KB
testcase_17	AC	2 ms; 4,380 KB
testcase_18	WA	-
testcase_19	AC	1 ms; 4,376 KB
testcase_20	AC	108 ms; 4,380 KB
testcase_21	AC	38 ms; 4,380 KB
testcase_22	AC	392 ms; 4,380 KB
testcase_23	AC	423 ms; 4,376 KB
testcase_24	AC	119 ms; 4,376 KB
testcase_25	AC	74 ms; 4,380 KB
testcase_26	AC	1,131 ms; 4,380 KB
testcase_27	AC	8 ms; 4,380 KB
testcase_28	AC	1,053 ms; 4,380 KB
testcase_29	AC	135 ms; 4,384 KB
testcase_30	AC	1,096 ms; 4,384 KB
testcase_31	AC	17 ms; 4,380 KB
testcase_32	AC	7 ms; 4,376 KB
testcase_33	AC	160 ms; 4,380 KB
testcase_34	AC	337 ms; 4,380 KB
testcase_35	AC	1,134 ms; 4,376 KB
testcase_36	AC	61 ms; 4,380 KB
testcase_37	AC	716 ms; 4,380 KB
testcase_38	AC	522 ms; 4,380 KB
testcase_39	AC	1,073 ms; 4,384 KB
testcase_40	WA	-
testcase_41	AC	2 ms; 4,380 KB
testcase_42	WA	-
testcase_43	WA	-

ソースコード

#include <bits/stdc++.h>

#define rep(i,n) for(int i=0;i<(n);i++)

using namespace std;

class Eratosthenes_sieve{
	vector<bool> er;
	vector<int> p;
public:
	Eratosthenes_sieve(int n):er(n+1,true){
		if(n>=0) er[0]=false;
		if(n>=1) er[1]=false;
		for(int i=2;i*i<=n;i++) if(er[i]) for(int j=i*i;j<=n;j+=i) er[j]=false;
		rep(i,n+1) if(er[i]) p.emplace_back(i);
	}

	vector<int> primes()const{ return p; }

	bool is_prime(int a)const{
		assert(a<=(int)er.size()-1);
		return a>=0 && er[a];
	}
};

int main(){
	int n; scanf("%d",&n);
	Eratosthenes_sieve E(n);

	int ans=0;
	auto p=E.primes();
	rep(i,p.size()){
		for(int j=i;j<p.size();j++) {
			int r2=p[i]+p[j];
			if(r2>n) break;
			int r=sqrt(r2)+0.5;
			if(r*r==r2 && E.is_prime(r)) ans+=(i==j?1:2);
		}
	}
	printf("%d\n",ans);

	return 0;
}