結果

問題 No.732 3PrimeCounting
ユーザー fura
提出日時 2020-06-10 21:13:20
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 349 ms / 3,000 ms
コード長 817 bytes
コンパイル時間 2,694 ms
コンパイル使用メモリ 197,144 KB
最終ジャッジ日時 2025-01-11 00:54:43
ジャッジサーバーID
(参考情報)
judge5 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 89
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:28:21: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   28 |         int n; scanf("%d",&n);
      |                ~~~~~^~~~~~~~~

ソースコード

diff #

#include <bits/stdc++.h>

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

using namespace std;
using lint=long long;

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 ES(3*n);
	auto P=ES.primes();

	lint ans=0;
	vector<int> cnt(3*n+1);
	rep(k,P.size()){
		if(P[k]>n) break;
		for(int l=k+1;l<P.size();l++) ans+=cnt[P[l]-P[k]];

		int j=k;
		rep(i,j) cnt[P[i]+P[j]]++;
	}
	printf("%lld\n",ans);

	return 0;
}
0