結果

問題 No.843 Triple Primes
コンテスト
ユーザー Naoto Fujiwara
提出日時 2023-05-17 12:20:12
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 7 ms / 2,000 ms
コード長 1,222 bytes
コンパイル時間 952 ms
コンパイル使用メモリ 87,960 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-12-15 07:07:59
合計ジャッジ時間 2,684 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<set>
#include<algorithm>
#include<vector>
#include<string>
#include<set>
#include<map>
#include<numeric>
#include<queue>
#include<cmath>
using namespace std;
typedef long long ll;
const ll INF=1LL<<60;
typedef pair<int,int> P;
typedef pair<int,P> PP;
const ll MOD=998244353;

vector<ll> prime(ll n){
    vector<bool> isprime(n+1,true);
    isprime[0]=isprime[1]=false;
    vector<ll> res;
    for(int i=2;i<=n;i++){
        if(isprime[i]){
            res.push_back(i);
            for(ll p=2*i;p<=n;p+=i){
                isprime[p]=false;
            }
        }
    }

    return res;
}

int main(){
    ll N;
    cin>>N;

    if(N==1){
        cout<<0<<endl;
        return 0;
    }

    //cout<<"sz="<<prime(500000).size()<<endl;//素数の数41538

    vector<ll> pr=prime(N);

    vector<bool> isprime(N+1,false);
    for(ll pnum: pr){
        isprime[pnum]=true;
    }

    ll cnt=0;
    //p=2
    for(ll r:pr){
        if(r*r-2>N) continue;
        if(r*r-2>=0 && isprime[r*r-2])cnt++;
    }

    //q=2
    for(ll r:pr){
        if(r*r-2>N) continue;
        if(r*r-2>=0 && isprime[r*r-2])cnt++;
    }

    cnt--;//p=2,q=2を2回カウントしているので
    cout<<cnt<<endl;



}
0