結果
問題 | No.917 Make One With GCD |
ユーザー |
![]() |
提出日時 | 2019-10-25 22:32:55 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 184 ms / 2,000 ms |
コード長 | 1,435 bytes |
コンパイル時間 | 180 ms |
コンパイル使用メモリ | 82,368 KB |
実行使用メモリ | 70,976 KB |
最終ジャッジ日時 | 2024-06-24 11:04:45 |
合計ジャッジ時間 | 3,411 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 32 |
ソースコード
# coding: utf-8# Your code here!import sys#sys.setrecursionlimit(10**6)readline = sys.stdin.readlinedef prime_factorize(N): #素因数分解exponent = 0while N%2 == 0:exponent += 1N //= 2if exponent: factorization = [[2,exponent]]else: factorization = []i=1while i*i <=N:i += 2if N%i: continueexponent = 0while N%i == 0:exponent += 1N //= ifactorization.append([i,exponent])if N!= 1: factorization.append([N,1])assert N != 0, "zero"return factorizationdef divisor_list(N): #約数のリストif N == 1: return [1]res = []for i in range(1,N):if i*i >= N: breakif N%i == 0:res.append(i)res.append(N//i)if i*i == N: res.append(i)return resn = int(input())a = [int(i) for i in readline().split()]res = set()pres = set()for i in a:res|=set(divisor_list(i))P = prime_factorize(i)P = [pi[0] for pi in P]pres|=set(P)ans = {i:0 for i in res}for i in a: ans[i] += 1res = sorted(list(res),reverse=True)for p in pres:for c in res:if c%p==0: ans[c//p] += ans[c]p2 = [1]for i in range(100): p2.append(p2[-1]*2)for i in ans:ans[i] = p2[ans[i]]-1res = list(reversed(res))for p in pres:for c in res:if c*p in res: ans[c] -= ans[c*p]print(ans[1])