結果
| 問題 | No.3776 Magic Square of Product |
| コンテスト | |
| ユーザー |
titia
|
| 提出日時 | 2026-10-24 03:13:11 |
| 言語 | PyPy3 (8.0.0 + ACL) |
| 結果 |
AC
不安定
|
| 実行時間 | 141 ms / 2,500 ms |
| + 53µs | |
| コード長 | 1,211 bytes |
| 記録 | |
| コンパイル時間 | 92 ms |
| コンパイル使用メモリ | 64,356 KB |
| 実行使用メモリ | 72,988 KB |
| 最終ジャッジ日時 | 2026-10-24 03:13:18 |
| 合計ジャッジ時間 | 6,023 ms |
|
ジャッジサーバーID (参考情報) |
judge4_0 / judge3_0 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 25 |
ソースコード
import sys
input = sys.stdin.readline
X=[1, 6, 21, 55, 120, 231, 406, 666, 1035, 1540, 2211, 3081, 4186, 5565, 7260, 9316, 11781, 14706, 18145, 22155, 26796, 32131, 38226, 45150, 52975, 61776, 71631, 82621, 94830, 108345]
# エラトステネスの篩を用いた素因数分解・約数列挙
MAX=10**6+10 # 使いたい最大値を指定
# Sieve[i]で、iの最も小さい1でない約数を返す。
Sieve=[i for i in range(MAX)]
for i in range(2,MAX):
if Sieve[i]!=i:
continue
for j in range(i,MAX,i):
if Sieve[j]==j:
Sieve[j]=i
# 素因数分解
def fact(x):
D=dict()
while x!=1:
k=Sieve[x]
if k in D:
D[k]+=1
else:
D[k]=1
x//=k
return D
# 約数列挙
def faclist(x):
LIST=[1]
while x!=1:
k=Sieve[x]
count=0
while x%k==0:
count+=1
x//=k
LIST2=[]
for l in LIST:
for i in range(1,count+1):
LIST2.append(l*k**i)
LIST+=LIST2
return LIST
T=int(input())
for tests in range(T):
N=int(input())
F=fact(N)
ANS=1
for f in F:
ANS=ANS*X[F[f]]
print(ANS)
titia