結果
| 問題 |
No.2313 Product of Subsequence (hard)
|
| ユーザー |
sotanishy
|
| 提出日時 | 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])
sotanishy