結果

問題 No.1746 Sqrt Integer Segments
ユーザー chineristAC
提出日時 2021-11-18 01:20:43
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,264 bytes
コンパイル時間 1,763 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 130,548 KB
最終ジャッジ日時 2024-12-24 11:09:23
合計ジャッジ時間 9,902 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 14 WA * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys,random,bisect
from collections import deque,defaultdict
from heapq import heapify,heappop,heappush
from itertools import permutations
from math import log,gcd

input = lambda :sys.stdin.readline()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

M = 10**6
spf = [i for i in range(M+1)]
for i in range(2,M+1):
    if spf[i]==i:
        for j in range(i,M+1,i):
            spf[j] = i


def factorize(n):
    res = []
    while n!=1:
        p = spf[n]
        cnt = 0
        while n%p==0:
            cnt += 1
            n //= p
        if cnt&1:
            res.append(p)
    return res
    
N = int(input())
A = [1] + li()

mod = 998244353
r = random.randint(2,mod-2)
pow_r = [1] * (M+1)
for i in range(1,M+1):
    pow_r[i] = pow_r[i-1] * r % mod

dic = {}

PS = [0]
for i in range(N):
    tmp_p = factorize(A[i+1])
    tmp_PS = PS[-1]
    for p in tmp_p:
        if p not in dic:
            dic[p] = 0
        tmp_PS += pow_r[p] * (1-2*dic[p]) % mod
        tmp_PS %= mod
        dic[p] = 1 - dic[p]
    PS.append(tmp_PS)

res = 0
PS.sort()
val = PS[0]
tmp = 1
for i in range(1,N+1):
    if val==PS[i]:
        tmp += 1
    else:
        res += tmp * (tmp-1)//2
        val = PS[i]
        tmp = 1
res += tmp * (tmp-1)//2
print(res)
0