結果

問題 No.843 Triple Primes
ユーザー ahe100
提出日時 2019-06-28 22:07:05
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 914 bytes
コンパイル時間 1,683 ms
コンパイル使用メモリ 168,976 KB
実行使用メモリ 14,120 KB
最終ジャッジ日時 2024-07-02 04:48:14
合計ジャッジ時間 38,605 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 TLE * 1
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for(int i = 0;i<((int)(n));i++)
#define reg(i,a,b) for(int i = ((int)(a));i<=((int)(b));i++)
#define irep(i,n) for(int i = ((int)(n)-1);i>=0;i--)
#define ireg(i,a,b) for(int i = ((int)(b));i>=((int)(a));i--)
typedef long long ll;

/*
*/

ll n,a[1000010]={},ans=0;
vector<ll> pl;

void prime_list(ll n,vector<ll>& p){
	ll visited[500010]={};
	reg(i,2,n){
		if(visited[i]==0){
			p.push_back(i);
			reg(j,2,n){
				if(j*i>n)break;
				visited[i*j]=1;
			}
		}
	}
}

void init(){
	cin>>n;
	prime_list(n,pl);
	for(ll r:pl){
		if(r*r>1000000)break;
		// cerr<<r<<endl;
		a[r*r]=r;
	}
}

int main(void){
	init();
	ll m=pl.size();
	// cerr<<m<<endl;
	rep(i,m){
		ll p=pl[i];
		rep(j,i){
			ll q=pl[j];
			if(a[p+q]!=0){
				ans++;
				// cerr<<p<<" "<<q<<" "<<a[p+q]<<endl;
			}
		}
	}
	ans*=2;
	if(n>=2)ans++;//r=2でいける
	cout<<ans<<endl;
	return 0;
}
0