結果
| 問題 | No.1731 Product of Subsequence | 
| コンテスト | |
| ユーザー |  ぷら | 
| 提出日時 | 2021-11-05 21:28:48 | 
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) | 
| 結果 | 
                                WA
                                 
                             | 
| 実行時間 | - | 
| コード長 | 975 bytes | 
| コンパイル時間 | 2,370 ms | 
| コンパイル使用メモリ | 198,556 KB | 
| 最終ジャッジ日時 | 2025-01-25 11:57:40 | 
| ジャッジサーバーID (参考情報) | judge2 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 4 | 
| other | AC * 27 WA * 4 | 
ソースコード
#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;
        }
    }
    cout << dp[N][0] << endl;
}
            
            
            
        