結果
問題 | 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]=1 return prime prime = eratosthenes(32000) def factorization(n): res = [] for p in prime: cnt = 0 while n%p==0: cnt += 1 n //= p if cnt: res.append((p,cnt)) if n!=1: res.append((n,1)) return res def divisors_prime(prime): res = [1] for p,cnt in prime: newres = [] for d in res: t = d for j in range(cnt+1): newres.append(t) t *= p res = newres return res import sys,random,bisect from collections import deque,defaultdict from heapq import heapify,heappop,heappush from itertools import permutations from math import gcd,log input = lambda :sys.stdin.readline().rstrip() mi = lambda :map(int,input().split()) li = lambda :list(mi()) M = 10**9 N = 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] = 0 all_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: break else: 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]) break print(k) for i in range(k): S[i] = [len(S[i])] + S[i] print(*S[i])