結果

問題 No.2074 Product is Square ?
ユーザー shotoyooshotoyoo
提出日時 2022-09-16 22:44:28
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 464 ms / 2,000 ms
コード長 2,673 bytes
コンパイル時間 149 ms
コンパイル使用メモリ 81,820 KB
実行使用メモリ 70,696 KB
最終ジャッジ日時 2024-06-01 13:45:31
合計ジャッジ時間 8,313 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 43 ms
56,112 KB
testcase_01 AC 60 ms
65,468 KB
testcase_02 AC 60 ms
66,308 KB
testcase_03 AC 59 ms
66,092 KB
testcase_04 AC 64 ms
67,392 KB
testcase_05 AC 61 ms
66,768 KB
testcase_06 AC 60 ms
65,316 KB
testcase_07 AC 59 ms
65,656 KB
testcase_08 AC 59 ms
67,008 KB
testcase_09 AC 61 ms
66,348 KB
testcase_10 AC 68 ms
69,416 KB
testcase_11 AC 309 ms
66,204 KB
testcase_12 AC 285 ms
66,392 KB
testcase_13 AC 463 ms
70,184 KB
testcase_14 AC 309 ms
66,956 KB
testcase_15 AC 306 ms
66,596 KB
testcase_16 AC 287 ms
66,256 KB
testcase_17 AC 464 ms
70,696 KB
testcase_18 AC 316 ms
67,432 KB
testcase_19 AC 265 ms
66,312 KB
testcase_20 AC 240 ms
66,200 KB
testcase_21 AC 462 ms
69,228 KB
testcase_22 AC 285 ms
66,244 KB
testcase_23 AC 268 ms
66,340 KB
testcase_24 AC 282 ms
66,492 KB
testcase_25 AC 346 ms
68,276 KB
testcase_26 AC 232 ms
66,220 KB
testcase_27 AC 307 ms
67,160 KB
testcase_28 AC 287 ms
66,820 KB
testcase_29 AC 301 ms
67,072 KB
testcase_30 AC 241 ms
66,252 KB
testcase_31 AC 48 ms
61,140 KB
testcase_32 AC 43 ms
55,584 KB
testcase_33 AC 61 ms
68,292 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys, random
input = lambda : sys.stdin.readline().rstrip()


write = lambda x: sys.stdout.write(x+"\n"); writef = lambda x: print("{:.12f}".format(x))
debug = lambda x: sys.stderr.write(x+"\n")
YES="Yes"; NO="No"; pans = lambda v: print(YES if v else NO); INF=10**18
LI = lambda : list(map(int, input().split())); II=lambda : int(input())
def debug(_l_):
    for s in _l_.split():
        print(f"{s}={eval(s)}", end=" ")
    print()
def dlist(*l, fill=0):
    if len(l)==1:
        return [fill]*l[0]
    ll = l[1:]
    return [dlist(*ll, fill=fill) for _ in range(l[0])]

def isqrt(n):
    """x*x<=nなる最大のx
    """
    x,y = n, (n+1)//2
    while y<x:
        x,y = y, (y + n//y) // 2
    return x
from math import gcd
import random
def is_prime(n):
    """miller_rabinによる素数判定
    ※ 1は素数と扱う
    """
    l = [2,3,5,7,11,13,17,19,23,29,31,37]
    if n==1 or n in l:
        return True
    d = n-1
    s = 0
    while d%2==0:
        s += 1
        d //= 2
    for a in l:
        v = pow(a,d,n)
        if v==1 or v==n-1:
            continue
        for _ in range(s):
            v = v*v % n
            if v==n-1:
                break
        else:
            return False
    return True
def rho(n):
    """nを割り切る3以上の素数を返す(素数のときnを返す)
    """
    if is_prime(n):
        return n
    while True:
        x = y = random.randint(1,n-1)
        g = 1
        while g==1:
            x = (x*x - 3) % n
            y = (y*y - 3) % n
            y = (y*y - 3) % n
            g = gcd((x-y), n)
        if g>1:
            return rho(g)
def factor(n):
    """高速な素因数分解
    """
    n0 = n
    if n in _d:
        return _d[n0]
    if n==1:
        return {}
    f = is_prime(n)
    if f:
        _d[n0] = {n:1}
        return {n:1}
    ans = {}
    while n%2==0:
        ans.setdefault(2, 0)
        ans[2] += 1
        n //= 2
    v = rho(n)
    while v!=n and n>1:
        ans.setdefault(v, 0)
        while n%v==0:
            n //= v
            ans[v] += 1
        if n>3 and is_prime(n):
            ans.setdefault(n,0)
            ans[n] += 1
            _d[n0] = ans
            return ans
        v = rho(n)
    if n>1:
        ans.setdefault(n, 0)
        ans[n] += 1
    _d[n0] = ans
    return ans
t = II()
_d = {}
for _ in range(t):
    n = int(input())
    a = list(map(int, input().split()))
    from math import gcd
    for i in range(n):
        for j in range(i+1,n):
            g = gcd(a[i], a[j])
            a[i] //= g
            a[j] //= g
        if isqrt(a[i])**2!=a[i]:
            ans = 0
            break
    else:
        ans = 1
    pans(ans)
0