結果
問題 |
No.371 ぼく悪いプライムじゃないよ
|
ユーザー |
![]() |
提出日時 | 2016-05-14 16:38:14 |
言語 | Python2 (2.7.18) |
結果 |
WA
|
実行時間 | - |
コード長 | 961 bytes |
コンパイル時間 | 394 ms |
コンパイル使用メモリ | 6,912 KB |
実行使用メモリ | 13,856 KB |
最終ジャッジ日時 | 2024-10-06 02:39:32 |
合計ジャッジ時間 | 6,055 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 20 WA * 5 TLE * 2 -- * 15 |
ソースコード
import math L,H=map(int,raw_input().split()) upperBound = 100001 isPrime=[True for i in range((upperBound-3)/2+1)] for i in range((int(math.sqrt(upperBound-3))-3)/2+1): # k*k = (2i+3)^2 <= upperBound-3 if isPrime[i]: k = i+i+3 # this is prime number # false out non_prime number from k*k # k*k = 4i^2+12i+9 = 2*(2i^2+6i+3)+3 # k*(k+1) = 2*(2i^2+6i+3)+3 + 2i+3 ... even number. skip # k*(k+2) = 2*(2i^2+6i+3+k)+3 non_prime = k*(i+1)+i for j in range(non_prime,(upperBound-3)/2+1,k): isPrime[j]=False primes = [2] + [i+i+3 for i in range((upperBound-3)/2+1) if isPrime[i]==True] primes.reverse() for i,p in enumerate(primes): if H/p <= 1: continue for c in range((H/p)*p, L-1, -p): ok = True for j in range(i+1,len(primes)): q=primes[j] if c%q==0: ok=False break if ok: print c exit()