結果
| 問題 | 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();
    }
}
            
            
            
        