結果

問題 No.2075 GCD Subsequence
ユーザー LyricalMaestroLyricalMaestro
提出日時 2024-08-31 02:54:30
言語 PyPy3
(7.3.15)
結果
MLE  
実行時間 -
コード長 1,628 bytes
コンパイル時間 320 ms
コンパイル使用メモリ 82,424 KB
実行使用メモリ 522,280 KB
最終ジャッジ日時 2024-08-31 02:55:05
合計ジャッジ時間 32,777 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
52,524 KB
testcase_01 AC 45 ms
52,412 KB
testcase_02 AC 37 ms
52,532 KB
testcase_03 AC 37 ms
53,244 KB
testcase_04 AC 35 ms
52,636 KB
testcase_05 AC 39 ms
52,380 KB
testcase_06 AC 38 ms
52,464 KB
testcase_07 AC 37 ms
54,092 KB
testcase_08 AC 375 ms
137,880 KB
testcase_09 AC 523 ms
168,624 KB
testcase_10 AC 397 ms
139,236 KB
testcase_11 AC 453 ms
155,020 KB
testcase_12 AC 457 ms
147,540 KB
testcase_13 AC 364 ms
130,628 KB
testcase_14 AC 470 ms
155,324 KB
testcase_15 AC 360 ms
134,372 KB
testcase_16 AC 412 ms
138,484 KB
testcase_17 AC 520 ms
167,828 KB
testcase_18 AC 2,271 ms
499,456 KB
testcase_19 AC 2,314 ms
499,132 KB
testcase_20 AC 2,362 ms
499,224 KB
testcase_21 AC 2,357 ms
499,440 KB
testcase_22 AC 2,345 ms
499,304 KB
testcase_23 AC 2,296 ms
499,360 KB
testcase_24 AC 2,361 ms
499,456 KB
testcase_25 AC 2,846 ms
499,340 KB
testcase_26 AC 2,208 ms
499,368 KB
testcase_27 AC 2,339 ms
499,424 KB
testcase_28 MLE -
testcase_29 AC 34 ms
52,460 KB
testcase_30 AC 34 ms
52,268 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

## https://yukicoder.me/problems/no/2075

MOD = 998244353

def main():
    N = int(input())
    A = list(map(int, input().split()))

    # osa-k法による素因数となる素数を列挙
    max_a = max(A)
    if max_a == 1:
        print(N)
        return

    primes = [p for p in range(max_a + 1)]
    for p in range(2, max_a + 1):
        if primes[p] != p:
            continue

        x = 2 * p
        while x <= max_a:
            primes[x] = min(primes[x], p)
            x += p
    
    p_divisors = []
    for i in range(N):
        a = A[i]
        p_list = []
        while primes[a] != 1:
            if len(p_list) == 0 or p_list[-1] != primes[a]:
                p_list.append(primes[a])
            a //= primes[a]

        p_d_len = len(p_list)
        array = []
        for p_bit in range(1, 2 ** p_d_len):
            bit_count = 0
            a = 1
            for j in range(p_d_len):
                if p_bit & (1 << j) > 0:
                    a *= p_list[j]
                    bit_count += 1
            
            if bit_count % 2 == 1:
                array.append((a, 1))
            else:
                array.append((a, -1))
        p_divisors.append(array)

    # dpを使う
    primes_sum = [0] * (max_a + 1)
    answer = 0
    for i in range(N):
        ans = 1
        for a, bit_count in p_divisors[i]:
            ans += (primes_sum[a] * bit_count) % MOD
            ans %= MOD

        answer += ans
        answer %= MOD

        for a, _ in p_divisors[i]:
            primes_sum[a] += ans
            primes_sum[a] %= MOD

    print(answer)

if __name__ == "__main__":
    main()
0