結果
問題 | No.2081 Make a Test Case of GCD Subset |
ユーザー | taiga0629kyopro |
提出日時 | 2022-09-14 17:32:15 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 76 ms / 2,000 ms |
コード長 | 2,465 bytes |
コンパイル時間 | 238 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 77,824 KB |
最終ジャッジ日時 | 2024-05-08 21:43:43 |
合計ジャッジ時間 | 5,743 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 70 ms
77,312 KB |
testcase_01 | AC | 71 ms
77,312 KB |
testcase_02 | AC | 74 ms
77,312 KB |
testcase_03 | AC | 72 ms
77,696 KB |
testcase_04 | AC | 72 ms
77,312 KB |
testcase_05 | AC | 71 ms
77,312 KB |
testcase_06 | AC | 75 ms
77,696 KB |
testcase_07 | AC | 72 ms
77,184 KB |
testcase_08 | AC | 75 ms
77,056 KB |
testcase_09 | AC | 72 ms
77,312 KB |
testcase_10 | AC | 72 ms
77,312 KB |
testcase_11 | AC | 73 ms
77,056 KB |
testcase_12 | AC | 72 ms
77,184 KB |
testcase_13 | AC | 73 ms
77,696 KB |
testcase_14 | AC | 76 ms
77,568 KB |
testcase_15 | AC | 71 ms
77,824 KB |
testcase_16 | AC | 69 ms
77,312 KB |
testcase_17 | AC | 71 ms
77,312 KB |
testcase_18 | AC | 69 ms
77,184 KB |
testcase_19 | AC | 74 ms
77,184 KB |
testcase_20 | AC | 73 ms
77,312 KB |
testcase_21 | AC | 69 ms
76,672 KB |
testcase_22 | AC | 72 ms
77,184 KB |
testcase_23 | AC | 70 ms
77,568 KB |
testcase_24 | AC | 68 ms
77,056 KB |
testcase_25 | AC | 71 ms
77,568 KB |
testcase_26 | AC | 69 ms
77,568 KB |
testcase_27 | AC | 71 ms
77,696 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) prime=primes(2*10**5).prime m=int(input()) if m==0: print(1) print(1) exit() while prime[-1]>10**5:prime.pop() ind=0 ind2=-1 cost=10**9 res=[] for b in range(31): x=[] c=0 nokori=m if nokori-(2**b-1)>=0: x.append((b,0)) nokori-=2**b-1 c+=b for i in range(29,-1,-1): if nokori-2**i>=0: x.append((i,1)) nokori-=2**i c+=i if c<cost: cost=c res=x[:] ans=[] for num,f in res: for iii in range(num):ans.append(prime[ind]) ind+=1 if f: ans.append(prime[ind2]) ind2-=1 ans.sort(reverse=True) for i in range(len(ans)): if ans.count(ans[i])>1: ans[i]*=prime[ind] ind+=1 print(len(ans)) print(*ans)