結果
| 問題 | 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 | 
| コンパイル時間 | 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 | 
ソースコード
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])
            
            
            
        