結果

問題 No.1746 Sqrt Integer Segments
ユーザー shotoyooshotoyoo
提出日時 2021-11-18 20:45:02
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 1,919 bytes
コンパイル時間 313 ms
コンパイル使用メモリ 86,792 KB
実行使用メモリ 133,528 KB
最終ジャッジ日時 2023-08-27 14:44:06
合計ジャッジ時間 11,302 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
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)
LI = lambda : list(map(int, input().split()))
# sys.setrecursionlimit(3*10**5+10)

def hurui(n):
    """線形篩
    pl: 素数のリスト
    mpf: iを割り切る最小の素因数
    """
    pl = []
    mpf = [None]*(n+1)
    for d in range(2,n+1):
        if mpf[d] is None:
            mpf[d] = d
            pl.append(d)
        for p in pl:
            if p*d>n or p>mpf[d]:
                break
            mpf[p*d] = p
    return pl, mpf
from collections import defaultdict
def factor(num):
    d = defaultdict(int)
    if num==1:
        d.update({1:1})
        return d
    while num>1:
        d[mpf[num]] += 1
        num //= mpf[num]
    return d
def fs(num):
    f = factor(num)
    ans = [1]
    for k,v in f.items():
        tmp = []
        for i in range(len(ans)):
            val = 1
            for _ in range(v):
                val *= k
                ans.append(ans[i]*val)
    return ans
pl, mpf = hurui(10**6+10)
n = int(input())
a = list(map(int, input().split()))
# M = 10**9+7
# M = (1<<63) - 1
# M = 72057594037927931
# M = 92709568269121
vs = {}
d = {}
d[0] = 1
val = 0
s = set()
ans = 0
import random
b = random.randint(123,M-2)
def sub(k):
    if k in vs:
        return vs[k]
    vs[k] = random.randint(1,M-1)
#     vs[k] = pow(b, k, M)
    return vs[k]
for i in range(n):
    f = factor(a[i])
    for k,v in f.items():
        if k==1:
            continue
        if v%2==1:
            tmp = sub(k)
            if k in s:
                s.remove(k)
                val ^= tmp
            else:
                s.add(k)
                val ^= tmp
    ans += d.get(val, 0)
    d.setdefault(val, 0)
    d[val] += 1
print(ans)
0