結果

問題 No.1731 Product of Subsequence
ユーザー ぷらぷら
提出日時 2021-11-05 21:30:46
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 521 ms / 2,000 ms
コード長 1,018 bytes
コンパイル時間 2,425 ms
コンパイル使用メモリ 198,364 KB
最終ジャッジ日時 2025-01-25 11:58:26
ジャッジサーバーID
（参考情報） 		judge4 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

long long gcd(long long a,long long b) {
    if(a % b == 0) {
        return b;
    }
    else {
        return gcd(b,a % b);
    }
}

long long dp[2005][5005];
constexpr int mod = 1000000007;

int main() {
    int N,K;
    cin >> N >> K;
    vector<int>tmp;
    for(int i = 1; i*i <= K; i++) {
        if(K%i == 0) {
            tmp.push_back(i);
            if(i*i != K) {
                tmp.push_back(K/i);
            }
        }
    }
    sort(tmp.begin(),tmp.end());
    dp[0][tmp.size()-1] = 1;
    for(int i = 0; i < N; i++) {
        long long A;
        cin >> A;
        for(int j = 0; j < tmp.size(); j++) {
            long long nxt = tmp[j]/gcd(A,tmp[j]);
            int it = lower_bound(tmp.begin(),tmp.end(),nxt)-tmp.begin();
            dp[i+1][j] += dp[i][j];
            dp[i+1][j] %= mod;
            dp[i+1][it] += dp[i][j];
            dp[i+1][it] %= mod;
        }
    }
    if(K == 1) {
        dp[N][0]--;
    }
    cout << dp[N][0] << endl;
}
0