結果
| 問題 |
No.1731 Product of Subsequence
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-06-29 00:03:16 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 1,424 ms / 2,000 ms |
| コード長 | 1,690 bytes |
| コンパイル時間 | 1,966 ms |
| コンパイル使用メモリ | 176,552 KB |
| 実行使用メモリ | 169,088 KB |
| 最終ジャッジ日時 | 2024-06-29 00:03:32 |
| 合計ジャッジ時間 | 15,712 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 31 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define all(x) x.begin(),x.end()
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define rep_r(i,a,b) for(int i=a;i>=b;i--)
#define each(a,x) for (auto& x : a)
using pi = pair<int,int>;
using pl = pair<ll,ll>;
using vi = vector<int>;
using vl = vector<ll>;
#define sz(x) int(x.size())
#define so(x) sort(all(x))
#define so_r(x) sort(all(x),greater<int>())
#define lb lower_bound
#define ub upper_bound
const char nl = '\n';
int dx[4] = {1,-1,0,0};
int dy[4] = {0,0,1,-1};
int bit_cnt(int x){
return __builtin_popcount(x);
}
ll bex(ll a, ll b, ll mod = 1e9 + 7){ll res = 1LL; while(b){ if (b&1) res = res * a % mod; a = a * a % mod; b >>= 1;} return res;}
template<class t,class u> bool chmax(t&a,u b){if(a<b){a=b;return true;}else return false;}
template<class t,class u> bool chmin(t&a,u b){if(b<a){a=b;return true;}else return false;}
const int MOD = 1e9+7;
const int N = 5e5 + 5;
const ll INF = 1e16;
int n,k;
ll a[2005];
map<ll,ll> dp[2005];
void solve(){
cin >> n >> k;
rep(i,1,n) cin >> a[i];
dp[0][1] = 1;
rep(i,1,n){
ll G = __gcd(a[i], 1LL * k);
each(dp[i-1],p){
ll uoc = p.first;
ll g = __gcd(G * uoc, 1LL * k);
// choose a[i]
dp[i][g] += dp[i-1][uoc];
dp[i][g] %= MOD;
// not choose a[i]
dp[i][uoc] += dp[i-1][uoc];
dp[i][uoc] %= MOD;
}
}
if (k == 1) dp[n][k] = (dp[n][k] + MOD - 1) % MOD;
cout << dp[n][k] << nl;
}
int main(){
ios_base::sync_with_stdio(false);
cin.tie(NULL);
int t = 1;
// cin >> t;
while (t--){
solve();
}
}