結果
| 問題 |
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])