結果
問題 | No.843 Triple Primes |
ユーザー | LayCurse |
提出日時 | 2019-06-28 22:00:41 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 2,309 bytes |
コンパイル時間 | 2,168 ms |
コンパイル使用メモリ | 201,116 KB |
実行使用メモリ | 5,632 KB |
最終ジャッジ日時 | 2024-09-19 14:01:34 |
合計ジャッジ時間 | 3,471 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 5 ms
5,504 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 3 ms
5,376 KB |
testcase_07 | AC | 4 ms
5,376 KB |
testcase_08 | AC | 4 ms
5,376 KB |
testcase_09 | AC | 5 ms
5,504 KB |
testcase_10 | AC | 4 ms
5,376 KB |
testcase_11 | AC | 4 ms
5,504 KB |
testcase_12 | AC | 4 ms
5,504 KB |
testcase_13 | AC | 5 ms
5,504 KB |
testcase_14 | AC | 4 ms
5,376 KB |
testcase_15 | AC | 4 ms
5,376 KB |
testcase_16 | AC | 4 ms
5,376 KB |
testcase_17 | AC | 2 ms
5,376 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 2 ms
5,376 KB |
testcase_20 | AC | 2 ms
5,376 KB |
testcase_21 | AC | 2 ms
5,376 KB |
testcase_22 | AC | 3 ms
5,376 KB |
testcase_23 | AC | 3 ms
5,376 KB |
testcase_24 | AC | 3 ms
5,376 KB |
testcase_25 | AC | 2 ms
5,376 KB |
testcase_26 | AC | 5 ms
5,504 KB |
testcase_27 | AC | 3 ms
5,376 KB |
testcase_28 | AC | 5 ms
5,504 KB |
testcase_29 | AC | 3 ms
5,376 KB |
testcase_30 | AC | 5 ms
5,632 KB |
testcase_31 | AC | 2 ms
5,376 KB |
testcase_32 | AC | 3 ms
5,376 KB |
testcase_33 | AC | 3 ms
5,376 KB |
testcase_34 | AC | 4 ms
5,376 KB |
testcase_35 | AC | 5 ms
5,504 KB |
testcase_36 | AC | 3 ms
5,376 KB |
testcase_37 | AC | 4 ms
5,376 KB |
testcase_38 | AC | 3 ms
5,376 KB |
testcase_39 | AC | 5 ms
5,504 KB |
testcase_40 | AC | 2 ms
5,376 KB |
testcase_41 | AC | 2 ms
5,376 KB |
testcase_42 | WA | - |
testcase_43 | AC | 5 ms
5,504 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; void *wmem; template<class S, class T> inline S min_L(S a,T b){ return a<=b?a:b; } inline void rd(int &x){ int k, m=0; x=0; for(;;){ k = getchar_unlocked(); if(k=='-'){ m=1; break; } if('0'<=k&&k<='9'){ x=k-'0'; break; } } for(;;){ k = getchar_unlocked(); if(k<'0'||k>'9'){ break; } x=x*10+k-'0'; } if(m){ x=-x; } } inline void wt_L(char a){ putchar_unlocked(a); } inline void wt_L(int x){ char f[10]; int m=0, s=0; if(x<0){ m=1; x=-x; } while(x){ f[s++]=x%10; x/=10; } if(!s){ f[s++]=0; } if(m){ putchar_unlocked('-'); } while(s--){ putchar_unlocked(f[s]+'0'); } } int Prime_L(int N, int res[], void *mem=wmem){ bool *isprime; const int r=23000; int a, b, i, *sf, ss=1, sz=1; isprime = (bool*)mem; sf = (int*)(isprime + r); N /= 2; res[0] = 2; b =min_L(r, N); for(i=1;i<b;i++){ isprime[i] = 1; } for(i=1;i<b;i++){ if(isprime[i]){ res[sz++] = 2*i+1; sf[ss] = 2*i*(i+1); if(sf[ss] < N){ while(sf[ss] < r){ isprime[sf[ss]] = 0; sf[ss] += res[ss]; } ss++; } } } for(a=r; a<N; a+=r){ b =min_L(a + r, N); isprime -= r; for(i=a;i<b;i++){ isprime[i] = 1; } for(i=1;i<ss;i++){ while(sf[i] < b){ isprime[sf[i]] = 0; sf[i] += res[i]; } } for(i=a;i<b;i++){ if(isprime[i]){ res[sz++] = 2*i+1; } } } return sz; } char memarr[96000000]; int N; int ps; int p[42000]; int isp[500001]; int main(){ int i, j, k, m, res=0; wmem = memarr; rd(N); ps =Prime_L(N, p); for(i=0;i<ps;i++){ isp[p[i]] = 1; } if(N >= 2){ res++; } for(i=1;i<ps;i++){ j = p[i]*p[i]; if(j-2 > N){ break; } res += 2*isp[j-2]; } wt_L(res); wt_L('\n'); return 0; } // cLay varsion 20190626-1 // --- original code --- // int N; // int ps, p[42000]; // int isp[500001]; // { // int i, j, k, m, res = 0; // rd(N); // ps = Prime(N, p); // // rep(i,ps) isp[p[i]] = 1; // // if(N >= 2) res++; // rep(i,1,ps){ // j = p[i]*p[i]; // if(j-2 > N) break; // res += 2isp[j-2]; // } // // wt(res); // }