結果

問題 No.1731 Product of Subsequence
コンテスト
ユーザー srjywrdnprkt
提出日時 2023-06-03 01:14:13
言語 C++17
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 1,030 ms / 2,000 ms
コード長 1,046 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,152 ms
コンパイル使用メモリ 120,816 KB
最終ジャッジ日時 2025-02-13 21:55:58
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <vector>
#include <cmath>
#include <map>
#include <set>
#include <iomanip>
#include <queue>
#include <algorithm>
#include <numeric>
#include <deque>
#include <complex>
#include <cassert>

using namespace std;
using ll = long long;

vector<ll> factor;

void all_factor(ll n){
    for (ll i = 1; i*i <= n; i++){
        if (n % i == 0){
            factor.push_back(i);
            if (i*i != n) factor.push_back(n / i);
        }
    }
    sort(factor.begin(), factor.end());
}

const ll modc = 1e9+7;

int main(){

    ll N, K, g;
    cin >> N >> K;
    all_factor(K);
    vector<ll> A(N+1);
    vector<map<ll, ll>> dp(N+1);
    for (int i=1; i<=N; i++){
        cin >> A[i];
        A[i] = gcd(A[i], K);
    }
    dp[0][1] = 1;
    for (int i=1; i<=N; i++){
        for (auto [x, y] : dp[i-1]){
            dp[i][x] += y;
            dp[i][x] %= modc;
            g = gcd(A[i]*x, K);
            dp[i][g] += y;
            dp[i][g] %= modc;
        }
    }

    cout << dp[N][K] - (K == 1) << endl;

    return 0;
}
0