結果
問題 | No.1731 Product of Subsequence |
ユーザー |
![]() |
提出日時 | 2021-11-20 01:20:27 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 418 ms / 2,000 ms |
コード長 | 966 bytes |
コンパイル時間 | 527 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 87,680 KB |
最終ジャッジ日時 | 2024-11-08 05:14:34 |
合計ジャッジ時間 | 6,756 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 31 |
ソースコード
import math from sys import stdin n, m, *a = map(int, stdin.read().split()) p = 10**9+7 divideset = set() for i in range(1,int(math.floor(math.sqrt(m)))+2): if m % i == 0: divideset.add(i) divideset.add(m//i) dividelist = list(divideset) dividelist.sort() dividedict = dict() nd = len(dividelist) for i in range(nd): dividedict[dividelist[i]] = i check = True for i in range(n): kari = math.gcd(a[i],m) a[i] = kari if kari != m: check = False if check: ans = 1 for i in range(n): ans = ans * 2 % p ans = (ans - 1) % p print("{}".format(ans)) exit() dp = [[0 for i in range(nd)] for j in range(n+1)] dp[0][0] = 1 for i in range(n): for j in range(nd): kari = a[i] * dividelist[j] val = math.gcd(kari,m) dp[i+1][j] = (dp[i+1][j] + dp[i][j]) % p dp[i+1][dividedict[val]] = (dp[i+1][dividedict[val]] + dp[i][j]) % p ans = dp[n][nd-1] print("{}".format(ans))