結果
問題 | No.12 限定された素数 |
ユーザー | なお |
提出日時 | 2014-10-06 23:56:09 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 182 ms / 5,000 ms |
コード長 | 2,673 bytes |
コンパイル時間 | 626 ms |
コンパイル使用メモリ | 72,084 KB |
実行使用メモリ | 8,356 KB |
最終ジャッジ日時 | 2024-05-03 08:51:55 |
合計ジャッジ時間 | 2,931 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 19 ms
8,172 KB |
testcase_01 | AC | 36 ms
8,204 KB |
testcase_02 | AC | 24 ms
8,200 KB |
testcase_03 | AC | 182 ms
8,192 KB |
testcase_04 | AC | 27 ms
8,276 KB |
testcase_05 | AC | 53 ms
8,180 KB |
testcase_06 | AC | 106 ms
8,276 KB |
testcase_07 | AC | 119 ms
8,240 KB |
testcase_08 | AC | 43 ms
8,268 KB |
testcase_09 | AC | 42 ms
8,228 KB |
testcase_10 | AC | 36 ms
8,312 KB |
testcase_11 | AC | 155 ms
8,192 KB |
testcase_12 | AC | 116 ms
8,232 KB |
testcase_13 | AC | 85 ms
8,280 KB |
testcase_14 | AC | 44 ms
8,240 KB |
testcase_15 | AC | 86 ms
8,280 KB |
testcase_16 | AC | 132 ms
8,316 KB |
testcase_17 | AC | 22 ms
8,256 KB |
testcase_18 | AC | 18 ms
8,248 KB |
testcase_19 | AC | 19 ms
8,272 KB |
testcase_20 | AC | 24 ms
8,348 KB |
testcase_21 | AC | 34 ms
8,300 KB |
testcase_22 | AC | 22 ms
8,196 KB |
testcase_23 | AC | 21 ms
8,356 KB |
testcase_24 | AC | 23 ms
8,172 KB |
testcase_25 | AC | 53 ms
8,244 KB |
ソースコード
#include <iostream> #include <cmath> #include <set> using namespace std; #define REP(i, n) for(int(i)=0;(i)<(n);++(i)) const int LIM = 5000000; bool isprime[LIM+1]; void sieve_of_atkin(int N) { int n; int sqrtN = (int)sqrt((double)N); for(int z = 1; z <= 5; z += 4){ for(int y = z; y <= sqrtN; y += 6){ for(int x = 1; x <= sqrtN && (n = 4*x*x+y*y) <= N; ++x) isprime[n] = !isprime[n]; for(int x = y+1; x <= sqrtN && (n = 3*x*x-y*y) <= N; x += 2) isprime[n] = !isprime[n]; } } for(int z = 2; z <= 4; z += 2){ for(int y = z; y <= sqrtN; y += 6){ for(int x = 1; x <= sqrtN && (n = 3*x*x+y*y) <= N; x += 2) isprime[n] = !isprime[n]; for(int x = y+1; x <= sqrtN && (n = 3*x*x-y*y) <= N; x += 2) isprime[n] = !isprime[n]; } } for(int y = 3; y <= sqrtN; y += 6){ for(int z = 1; z <= 2; ++z){ for(int x = z; x <= sqrtN && (n = 4*x*x+y*y) <= N; x += 3) isprime[n] = !isprime[n]; } } for(int n = 5; n <= sqrtN; ++n) if(isprime[n]) for(int k = n*n; k <= N; k+=n*n) isprime[k] = false; isprime[2] = isprime[3] = true; } int main(){ int N, n; set<int> setN; cin >> N; REP(i,N){ cin >> n; setN.insert(n); } sieve_of_atkin(LIM); set<int> setK; int low = 1; int maxv = -1; for(int i = 1; i <= LIM; i++){ if(!isprime[i]) continue; int j = i; bool flag = true; set<int> add; while(j){ if(setN.count(j%10)){ add.insert(j%10); j /= 10; } else { flag = false; break; } } if(!flag){ if(low && setK.size() == N){ maxv = max(maxv, (i-1)-low); #ifdef _DEBUG if(maxv == (i-1)-low){ cout << low << ":" << i-1 << ":" << ((i-1)-low) << endl; for(int j = low; j <= i-1; j++){ if(isprime[j]) cout << j << " "; } cout << endl; } #endif } low = i+1; setK.clear(); } else { setK.insert(add.begin(), add.end()); } } if(setK.size() == N){ maxv = max(maxv, LIM-low); #ifdef _DEBUG if(maxv == (LIM-low)){ cout << low << ":" << LIM << ":" << (LIM-low) << endl; for(int j = low; j <= LIM; j++){ if(isprime[j]) cout << j << " "; } cout << endl; } #endif } cout << maxv << endl; return 0; }