結果

問題 No.1731 Product of Subsequence
コンテスト
ユーザー tktk_snsn
提出日時 2021-11-06 09:03:13
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
TLE  
実行時間 -
コード長 975 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 91 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 24,912 KB
最終ジャッジ日時 2024-11-07 03:38:06
合計ジャッジ時間 6,659 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 4
other TLE * 1 -- * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

from collections import defaultdict


def prime_factorization(n):  # 素因数分解
    res = []
    for i in range(2, int(n ** 0.5) + 1):
        if n % i == 0:
            cnt = 0
            while n % i == 0:
                cnt += 1
                n //= i
            res.append((i, cnt))
    if n > 1:
        res.append((n, 1))
    return sorted(res)


mod = 998244353
N, K = map(int, input().split())
A = list(map(int, input().split()))
P, F = zip(*prime_factorization(K))
M = len(P)

dp = defaultdict(int)
dp[(0,) * M] = 1
for a in A:
    cnt = [0] * M
    for i in range(M):
        while a % P[i] == 0:
            a //= P[i]
            cnt[i] += 1
    ndp = defaultdict(int)
    for key, val in dp.items():
        ndp[key] += val
        ndp[key] %= mod

        nkey = cnt[:]
        for i, k in enumerate(key):
            nkey[i] = min(F[i], k + nkey[i])
        nkey = tuple(nkey)
        ndp[nkey] += val
        ndp[nkey] %= mod
    dp = ndp
print(dp[F])
0