結果

問題 No.2313 Product of Subsequence (hard)
ユーザー LyricalMaestro
提出日時 2026-04-26 02:41:43
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
WA  
実行時間 -
コード長 2,091 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 129 ms
コンパイル使用メモリ 85,376 KB
実行使用メモリ 190,596 KB
最終ジャッジ日時 2026-04-26 02:42:22
合計ジャッジ時間 9,231 ms
ジャッジサーバーID
（参考情報） 		judge1_1 / judge3_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 8 WA * 2 TLE * 1 -- * 16
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

import math

MOD = 998244353

def main():
    N, K = map(int, input().split())
    A = list(map(int, input().split()))
    if K == 1:
        return (pow(2, N, MOD) - 1) % MOD

    base_k = K

    # Kを素因数分解
    sqrt_k = int(math.sqrt(K))
    k_primes = {}
    for p in range(2, sqrt_k + 1):
        if K % p == 0:
            ans = 0
            while K % p == 0:
                ans += 1
                K //= p
            k_primes[p] = ans
    if K > 1:
        k_primes[K] = 1
    
    k_primes_array = [(p, value) for p, value in k_primes.items()]
    k_primes_array.sort(key=lambda x: x[0])

    # 状態を全て定義
    def dfs(k_primes_array, index, p, p_list, p_to_list_map):
        if index == len(k_primes_array):
            xx = tuple(p_list)
            p_to_list_map[p] = xx
            return

        p0, value = k_primes_array[index]
        for v in range(value + 1):
            p_list.append(v)
            dfs(k_primes_array, index + 1, p * (p0 ** v), p_list, p_to_list_map)
            p_list.pop()
        return

    p_to_list_map = {}
    dfs(k_primes_array, 0, 1, [], p_to_list_map)

    # 各要素について
    a_array = []
    for a in A:
        ans = []
        for p, value in k_primes_array:
            x = 0
            while a % p == 0:
                x += 1
                a //= p
            x = min(value, x)
            ans.append(x)
        a_array.append(tuple(ans))
    
    # dpで解く
    dp = {1: 1}
    for a in a_array:
        new_dp = dp.copy()
        for key_p, value in dp.items():
            key_p_list = p_to_list_map[key_p]
            q = 1
            for j in range(len(a)):
                p, v = k_primes_array[j]
                q *= p ** (min(v, key_p_list[j] + a[j]))

            if q not in new_dp:
                new_dp[q] = 0
            new_dp[q] += value
            new_dp[q] %= MOD            

        dp = new_dp
    
    if K not in dp:
        print(0)
    else:
        print(dp[base_k])



            






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