結果

問題 No.368 LCM of K-products
ユーザー yaoshimaxyaoshimax
提出日時 2016-05-07 18:38:08
言語 Python2
(2.7.18)
結果
AC  
実行時間 1,983 ms / 2,000 ms
コード長 1,146 bytes
コンパイル時間 80 ms
コンパイル使用メモリ 7,040 KB
実行使用メモリ 20,036 KB
最終ジャッジ日時 2024-09-22 14:48:11
合計ジャッジ時間 30,316 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,439 ms
19,676 KB
testcase_01 AC 1,890 ms
19,956 KB
testcase_02 AC 1,899 ms
19,700 KB
testcase_03 AC 1,886 ms
19,828 KB
testcase_04 AC 1,889 ms
19,836 KB
testcase_05 AC 1,983 ms
19,828 KB
testcase_06 AC 71 ms
19,740 KB
testcase_07 AC 74 ms
19,868 KB
testcase_08 AC 75 ms
19,744 KB
testcase_09 AC 76 ms
19,868 KB
testcase_10 AC 78 ms
19,616 KB
testcase_11 AC 75 ms
19,616 KB
testcase_12 AC 90 ms
19,740 KB
testcase_13 AC 1,045 ms
19,796 KB
testcase_14 AC 1,577 ms
19,812 KB
testcase_15 AC 1,746 ms
19,772 KB
testcase_16 AC 1,528 ms
19,632 KB
testcase_17 AC 1,291 ms
19,676 KB
testcase_18 AC 1,788 ms
19,824 KB
testcase_19 AC 507 ms
19,996 KB
testcase_20 AC 1,189 ms
19,800 KB
testcase_21 AC 338 ms
19,756 KB
testcase_22 AC 1,711 ms
19,820 KB
testcase_23 AC 76 ms
19,744 KB
testcase_24 AC 75 ms
19,616 KB
testcase_25 AC 84 ms
19,484 KB
testcase_26 AC 74 ms
19,744 KB
testcase_27 AC 88 ms
19,616 KB
testcase_28 AC 77 ms
19,752 KB
testcase_29 AC 86 ms
19,868 KB
testcase_30 AC 89 ms
19,744 KB
testcase_31 AC 90 ms
19,744 KB
testcase_32 AC 131 ms
19,616 KB
testcase_33 AC 730 ms
20,036 KB
testcase_34 AC 1,907 ms
19,692 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import math
N,K=map(int,raw_input().split())
A=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]
vec=[[] for i in range(101001)]
for a in A:
    for i,p in enumerate(primes):
        cnt=0
        while a%p==0:
            a/=p
            cnt+=1
        if cnt!=0:
            vec[i].append(cnt)
    if a!=1:
        vec[len(primes)].append(1)
        primes.append(a)
ans=1
mod=1000000007
for i,p in enumerate(primes):
    vec[i].sort()
    vec[i].reverse()
    for j,v in enumerate(vec[i]):
        if j>=K:
            break
        ans*=pow(p,v)
        ans%=mod
print ans
0