結果
| 問題 | No.1731 Product of Subsequence |
| ユーザー |
 Navier_Boltzmann
|
| 提出日時 | 2022-01-30 19:17:27 |
| 言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 927 bytes |
| コンパイル時間 | 83 ms |
| コンパイル使用メモリ | 10,964 KB |
| 実行使用メモリ | 124,008 KB |
| 最終ジャッジ日時 | 2023-09-02 01:46:46 |
| 合計ジャッジ時間 | 6,651 ms |
|
ジャッジサーバーID （参考情報） |
judge16 / judge15 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | TLE | - |
| testcase_01 | -- | - |
| testcase_02 | -- | - |
| testcase_03 | -- | - |
| testcase_04 | -- | - |
| testcase_05 | -- | - |
| testcase_06 | -- | - |
| testcase_07 | -- | - |
| testcase_08 | -- | - |
| testcase_09 | -- | - |
| testcase_10 | -- | - |
| testcase_11 | -- | - |
| testcase_12 | -- | - |
| testcase_13 | -- | - |
| testcase_14 | -- | - |
| testcase_15 | -- | - |
| testcase_16 | -- | - |
| testcase_17 | -- | - |
| testcase_18 | -- | - |
| testcase_19 | -- | - |
| testcase_20 | -- | - |
| testcase_21 | -- | - |
| testcase_22 | -- | - |
| testcase_23 | -- | - |
| testcase_24 | -- | - |
| testcase_25 | -- | - |
| testcase_26 | -- | - |
| testcase_27 | -- | - |
| testcase_28 | -- | - |
| testcase_29 | -- | - |
| testcase_30 | -- | - |
| testcase_31 | -- | - |
| testcase_32 | -- | - |
| testcase_33 | -- | - |
| testcase_34 | -- | - |
ソースコード
def make_divisors(n):
lower_divisors , upper_divisors = [], []
i = 1
while i*i <= n:
if n % i == 0:
lower_divisors.append(i)
if i != n // i:
upper_divisors.append(n//i)
i += 1
return lower_divisors + upper_divisors[::-1]
from math import gcd
mod = 10**9 + 7
N,K = map(int,input().split())
A = list(map(int,input().split()))
if K == 1:
print((pow(2,N,mod)-1)%mod)
exit()
from collections import defaultdict
dp = [defaultdict(lambda:0) for _ in range(N)]
P = make_divisors(K)
P.sort()
n = len(P)
for i in range(n-1,-1,-1):
if A[0]%P[i] == 0:
dp[0][P[i]] = 1
break
for i in range(N):
dp[i][1] += 1
for i in range(N-1):
for p in P:
dp[i+1][p] = dp[i][p]
dp[i+1][p] %= mod
for p in P:
g = gcd(p*A[i+1],K)
dp[i+1][g] += dp[i][p]
dp[i+1][g] %= mod
print(dp[-1][P[-1]])