結果
| 問題 | No.1731 Product of Subsequence | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2022-01-30 19:16:18 | 
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) | 
| 結果 | 
                                TLE
                                 
                             | 
| 実行時間 | - | 
| コード長 | 873 bytes | 
| コンパイル時間 | 133 ms | 
| コンパイル使用メモリ | 12,928 KB | 
| 実行使用メモリ | 108,416 KB | 
| 最終ジャッジ日時 | 2024-06-11 08:24:35 | 
| 合計ジャッジ時間 | 6,615 ms | 
| ジャッジサーバーID (参考情報) | judge2 / judge4 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | -- * 4 | 
| other | TLE * 1 -- * 30 | 
ソースコード
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
N,K = map(int,input().split())
A = list(map(int,input().split()))
from collections import defaultdict
dp = [defaultdict(lambda:0) for _ in range(N)]
P = make_divisors(K)
P.sort()
mod = 10**9 + 7
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]])
            
            
            
        