結果

問題 No.1661 Sum is Prime (Hard Version)
ユーザー ace_amuroace_amuro
提出日時 2021-08-31 17:54:10
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 890 ms / 3,000 ms
コード長 1,456 bytes
コンパイル時間 834 ms
コンパイル使用メモリ 82,320 KB
実行使用メモリ 7,168 KB
最終ジャッジ日時 2024-05-04 14:28:38
合計ジャッジ時間 9,319 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
5,248 KB
testcase_01 AC 3 ms
5,376 KB
testcase_02 AC 723 ms
7,168 KB
testcase_03 AC 3 ms
5,376 KB
testcase_04 AC 3 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 3 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 3 ms
5,376 KB
testcase_09 AC 3 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 3 ms
5,376 KB
testcase_12 AC 600 ms
6,656 KB
testcase_13 AC 581 ms
6,528 KB
testcase_14 AC 612 ms
6,656 KB
testcase_15 AC 817 ms
7,040 KB
testcase_16 AC 723 ms
6,912 KB
testcase_17 AC 675 ms
6,912 KB
testcase_18 AC 284 ms
6,016 KB
testcase_19 AC 460 ms
6,272 KB
testcase_20 AC 314 ms
6,016 KB
testcase_21 AC 541 ms
7,040 KB
testcase_22 AC 408 ms
7,040 KB
testcase_23 AC 890 ms
7,040 KB
権限があれば一括ダウンロードができます

ソースコード

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