結果
問題 | No.1731 Product of Subsequence |
ユーザー |
![]() |
提出日時 | 2021-11-05 21:57:21 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,581 bytes |
コンパイル時間 | 3,021 ms |
コンパイル使用メモリ | 205,880 KB |
最終ジャッジ日時 | 2025-01-25 12:15:06 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 27 WA * 4 |
ソースコード
#include <bits/stdc++.h>#define rep(i,n) for(int i = 0; i < (int)(n); i++)using namespace std;using LL = long long;using P = pair<int,int>;using vv = vector<vector<int>>;const int INF = (int)1e9;const LL LINF = (LL)1e18;long long const mod = 1000000007;struct mint{long long val;mint(long long val = 0): val(val % mod) {}mint& operator += (const mint n){val += n.val;if(val >= mod) val -= mod;return *this;}mint& operator -= (const mint n){val -= n.val;if(val < 0) val += mod;return *this;}mint& operator *= (const mint n){val *= n.val;val %= mod;if(val < 0) val += mod;return *this;}mint operator + (const mint n) const{mint res(*this);return res += n;}mint operator - (const mint n) const{mint res(*this);return res -= n;}mint operator * (const mint n) const{mint res(*this);return res *= n;}mint pow(long long n) const{if(n == 0) return 1;mint m = pow(n >> 1);m *= m;if(n & 1) m *= *this;return m;}// mint division for prime modmint inv() const{return pow(mod - 2);}mint& operator /= (const mint n){return (*this) *= n.inv();}mint operator / (const mint n) const{mint res(*this);return res /= n;}};int main(){int N;LL K;cin >> N >> K;vector<LL> A(N);rep(i,N) cin >> A[i];vector<map<LL,mint>> dp(N+1);dp[0][K] = 1;rep(i,N){for(auto itr : dp[i]){LL res = itr.first;mint cnt = itr.second;dp[i+1][res] += cnt;LL g = gcd(res, A[i]);dp[i+1][res / g] += cnt;}}cout << dp[N][1].val << endl;return 0;}