結果
問題 | No.1711 Divide LCM |
ユーザー |
|
提出日時 | 2021-08-08 00:14:20 |
言語 | PyPy3 (7.3.15) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,836 bytes |
コンパイル時間 | 154 ms |
コンパイル使用メモリ | 82,420 KB |
実行使用メモリ | 73,244 KB |
最終ジャッジ日時 | 2024-09-17 17:02:12 |
合計ジャッジ時間 | 12,146 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | WA * 1 RE * 2 |
other | AC * 5 WA * 30 RE * 7 |
ソースコード
def eratosthenes(n):res=[0 for i in range(n+1)]prime=set([])for i in range(2,n+1):if not res[i]:prime.add(i)for j in range(1,n//i+1):res[i*j]=1return primeprime = eratosthenes(32000)def factorization(n):res = []for p in prime:cnt = 0while n%p==0:cnt += 1n //= pif cnt:res.append((p,cnt))if n!=1:res.append((n,1))return resdef divisors_prime(prime):res = [1]for p,cnt in prime:newres = []for d in res:t = dfor j in range(cnt+1):newres.append(t)t *= pres = newresreturn resimport sys,random,bisectfrom collections import deque,defaultdictfrom heapq import heapify,heappop,heappushfrom itertools import permutationsfrom math import gcd,loginput = lambda :sys.stdin.readline().rstrip()mi = lambda :map(int,input().split())li = lambda :list(mi())M = 10**9N = int(input())A = li()all_prime = {}all_divisors = set()for a in A:a_prime = factorization(a)for p,cnt in a_prime:if p not in all_prime:all_prime[p] = 0all_prime[p] = max(all_prime[p],pow(p,cnt))div = divisors_prime(a_prime)for d in div:all_divisors.add(d)P = [all_prime[p] for p in all_prime]P.sort(reverse=True)n = len(P)selected_P = [P[0]]prod = P[0]for i in range(1,n):selected_P.append(P[i])prod *= P[i]if prod > M:breakelse:exit(print(-1))k = len(selected_P)S = [[] for i in range(k)]for i in range(N):for j in range(k):if A[i]%selected_P[j]:S[j].append(A[i])breakprint(k)for i in range(k):S[i] = [len(S[i])] + S[i]print(*S[i])