結果

問題 No.1661 Sum is Prime (Hard Version)
コンテスト
ユーザー ace_amuro
提出日時 2021-08-31 17:54:10
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,088 ms / 3,000 ms
コード長 1,456 bytes
コンパイル時間 895 ms
コンパイル使用メモリ 82,492 KB
実行使用メモリ 7,168 KB
最終ジャッジ日時 2024-12-30 16:59:08
合計ジャッジ時間 11,944 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
using namespace std;
typedef long long LL;
const int M=1e6+10;
const LL MOD=1e9+7;

bool is_p[M];
int prime[M];
LL sqrtn;

int euler(int N) {
    int len = 0;
    memset(is_p, true, sizeof(is_p));
    for(int i = 2; i < N; i++) {
        if(is_p[i]) {
            len++;
            prime[len] = i;
        }
        for(int j = 1; j <= len && prime[j] <= i; j++) {
            if(i * prime[j] >= N) {
                break;
            }
            is_p[i * prime[j]] = false;
            if(i % prime[j] == 0) {
               
                break;
            } else {

            }
        }
    }
    return len;
}


LL L[M],R[M];
LL primepi(LL n){
    if(n<2) return 0;
    for(LL i=1;i<=sqrtn;++i)   R[i]=n/i-1;
    for(LL i=1;i<=sqrtn;++i)   L[i]=i-1;
    for(int s=1;prime[s]<=sqrtn;s++){
        LL ps=prime[s];
        for(LL i=1,tn=min(n/(ps*ps),sqrtn);i<=tn;++i){
            R[i] -= (i*ps<=sqrtn?R[i*ps]:L[n/(i*ps)])-L[ps-1];
        }
        for(LL i=sqrtn;i>=ps*ps;--i){
            L[i] -= L[i/ps]-L[ps-1];
        }
    }
    return R[1];
}


int main(){
    LL l,r;
    scanf("%lld%lld",&l,&r);
    sqrtn=ceil(sqrt(2*r));
    euler(sqrtn*2);
    LL ans1 = primepi(r)-primepi(l-1);
    LL ans2 = l==r?0:primepi(2*r-1)-primepi(2*l);
    printf("%lld\n",ans1+ans2);
    return 0;
}
0