結果

問題 No.843 Triple Primes
ユーザー furafura
提出日時 2020-05-05 16:37:29
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 796 bytes
コンパイル時間 2,937 ms
コンパイル使用メモリ 198,036 KB
最終ジャッジ日時 2025-01-10 06:50:08
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 38 WA * 4
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:27:21: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   27 |         int n; scanf("%d",&n);
      |                ~~~~~^~~~~~~~~

ソースコード

diff # 

#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;
}
0