結果
問題 | No.2313 Product of Subsequence (hard) |
ユーザー |
![]() |
提出日時 | 2023-05-24 21:40:20 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 3,489 ms / 4,000 ms |
コード長 | 833 bytes |
コンパイル時間 | 656 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 138,056 KB |
最終ジャッジ日時 | 2024-07-21 11:13:36 |
合計ジャッジ時間 | 36,734 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
from math import gcd import sys input = sys.stdin.readline def divisor(n): divisors = [] i = 1 while i * i < n: if n % i == 0: divisors.append(i) divisors.append(n // i) i += 1 if i * i == n: divisors.append(i) divisors.sort() return divisors mod = 998244353 N, K = map(int, input().split()) A = list(map(int, input().split())) if K == 1: print((pow(2, N, mod) - 1) % mod) exit() div = divisor(K) idx = {d: i for i, d in enumerate(div)} D = len(div) prod_idx = [[0]*D for _ in range(D)] for i in range(D): for j in range(D): prod_idx[i][j] = idx[gcd(div[i]*div[j], K)] dp = [0]*D dp[0] = 1 for x in A: i = idx[gcd(x, K)] for j in range(D)[::-1]: k = prod_idx[i][j] dp[k] = (dp[k]+dp[j]) % mod print(dp[-1])