結果
問題 | No.2081 Make a Test Case of GCD Subset |
ユーザー | taiga0629kyopro |
提出日時 | 2022-09-06 11:32:52 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 55 ms / 2,000 ms |
コード長 | 2,241 bytes |
コンパイル時間 | 164 ms |
コンパイル使用メモリ | 82,240 KB |
実行使用メモリ | 71,336 KB |
最終ジャッジ日時 | 2024-05-03 11:12:13 |
合計ジャッジ時間 | 4,564 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 52 ms
70,928 KB |
testcase_01 | AC | 49 ms
70,808 KB |
testcase_02 | AC | 51 ms
71,036 KB |
testcase_03 | AC | 49 ms
69,540 KB |
testcase_04 | AC | 53 ms
71,120 KB |
testcase_05 | AC | 53 ms
70,452 KB |
testcase_06 | AC | 51 ms
69,792 KB |
testcase_07 | AC | 51 ms
70,284 KB |
testcase_08 | AC | 50 ms
69,748 KB |
testcase_09 | AC | 49 ms
70,360 KB |
testcase_10 | AC | 49 ms
70,564 KB |
testcase_11 | AC | 51 ms
70,344 KB |
testcase_12 | AC | 55 ms
70,120 KB |
testcase_13 | AC | 51 ms
70,464 KB |
testcase_14 | AC | 51 ms
70,976 KB |
testcase_15 | AC | 51 ms
70,108 KB |
testcase_16 | AC | 51 ms
71,108 KB |
testcase_17 | AC | 51 ms
70,196 KB |
testcase_18 | AC | 50 ms
70,956 KB |
testcase_19 | AC | 55 ms
70,468 KB |
testcase_20 | AC | 51 ms
70,656 KB |
testcase_21 | AC | 52 ms
70,552 KB |
testcase_22 | AC | 51 ms
70,728 KB |
testcase_23 | AC | 50 ms
70,620 KB |
testcase_24 | AC | 50 ms
71,336 KB |
testcase_25 | AC | 51 ms
70,788 KB |
testcase_26 | AC | 50 ms
70,592 KB |
testcase_27 | AC | 52 ms
70,704 KB |
ソースコード
class primes(): def __init__(self, n): self.prime_num = n self.min_prime = [-1] * (self.prime_num + 1) # 2以上の自然数に対して最小の素因数を表す self.min_prime[0] = 0 self.min_prime[1] = 1 i = 2 self.prime = [] self.memo_prifac = {} while i <= self.prime_num: if self.min_prime[i] == -1: self.min_prime[i] = i self.prime.append(i) for j in self.prime: if i * j > self.prime_num or j > self.min_prime[i]: break self.min_prime[j * i] = j i += 1 def prifac(self, n): # 素因数分解した結果を返す if n in self.memo_prifac: return self.memo_prifac[n] res = {} x = n while x > 1: p = self.min_prime[x] if p in res: res[p] += 1 else: res[p] = 1 x //= p # self.memo_prifac[n] = res #場合によってはこの行を消すと高速化 return res def divisors(self, n): # 約数列挙 メモした方がいいかも if n== 1: return [1] prf = self.prifac(n) keys = [key for key in prf] def divsearch(i): if i == len(keys) - 1: return [keys[i] ** j for j in range(prf[keys[i]] + 1)] else: res = [] subres = divsearch(i + 1) p = keys[i] for j in range(prf[p] + 1): for node in subres: res.append(node * p ** j) return res return divsearch(0) pri=primes(10**5+10) def make(m): if m==0: m=998244353 p=pri.prime ind=0 ind2=-1 while p[ind2]>10**5:ind2-=1 a=[] for i in range(29,-1,-1): if m-2**i>=0: for _ in range(i): a.append(p[ind]) ind+=1 a.append(p[ind2]) ind2-=1 m-=2**i n=len(a) for i in range(n-1,-1,-1): if a.count(a[i])>=2: a[i]*=p[ind] ind+=1 return n,a m=int(input()) n,a=make(m) print(n) print(*a)