結果

問題 No.1036 Make One With GCD 2
ユーザー 👑 SPD_9X2SPD_9X2
提出日時 2020-04-24 22:18:56
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
RE  
実行時間 -
コード長 1,066 bytes
コンパイル時間 165 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 74,224 KB
最終ジャッジ日時 2024-04-24 17:47:58
合計ジャッジ時間 8,939 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
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 -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
testcase_38 RE -
testcase_39 RE -
testcase_40 RE -
testcase_41 RE -
testcase_42 RE -
testcase_43 RE -
testcase_44 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 


#RMQ用のSparseTableを作成
import math
import fractions
def make_SparseTable(lis):

    spt_len = 0
    while spt_len ** 2 < len(lis):
        spt_len += 1

    N = len(lis)
    ret = [[float("inf")] * spt_len for i in range(N)]

    for i in range(N):
        ret[i][0] = lis[i]
    
    for j in range(spt_len-1):
        j += 1

        for i in range(N):

            if i + 2**(j-1) < N:
                ret[i][j] = fractions.gcd(ret[i][j-1] , ret[i+2**(j-1)][j-1])
            else:
                ret[i][j] = ret[i][j-1]

    return ret

#[a,b]の最小値を求める
def RMQ_SPT(a,b,spt):

    if a == 0 and b == N-1:
        return al

    qlen = b-a+1
    ind = int(math.log(qlen,2))

    return fractions.gcd(spt[a][ind] , spt[b-2**ind+1][ind])

N = int(input())

A = list(map(int,input().split()))

al = A[0]
for i in A:
    al = fractions.gcd(i,al)
    

spt = make_SparseTable(A)

l = 0

ans = N*(N-1)//2 + N

for r in range(N):

    if A[r] == 1:
        continue

    while RMQ_SPT(l,r,spt) == 1:
        l += 1

    ans -= r-l+1

print (ans)
0