結果
問題 | No.1731 Product of Subsequence |
ユーザー | miscalc |
提出日時 | 2021-11-05 21:54:47 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,018 ms / 2,000 ms |
コード長 | 1,914 bytes |
コンパイル時間 | 2,311 ms |
コンパイル使用メモリ | 219,344 KB |
実行使用メモリ | 13,960 KB |
最終ジャッジ日時 | 2023-08-07 23:05:24 |
合計ジャッジ時間 | 12,021 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge11 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 638 ms
13,280 KB |
testcase_01 | AC | 2 ms
4,376 KB |
testcase_02 | AC | 1 ms
4,376 KB |
testcase_03 | AC | 2 ms
4,376 KB |
testcase_04 | AC | 1 ms
4,380 KB |
testcase_05 | AC | 2 ms
4,376 KB |
testcase_06 | AC | 2 ms
4,376 KB |
testcase_07 | AC | 2 ms
4,376 KB |
testcase_08 | AC | 53 ms
4,376 KB |
testcase_09 | AC | 3 ms
4,376 KB |
testcase_10 | AC | 1,018 ms
13,904 KB |
testcase_11 | AC | 853 ms
12,348 KB |
testcase_12 | AC | 5 ms
4,380 KB |
testcase_13 | AC | 22 ms
4,376 KB |
testcase_14 | AC | 788 ms
11,220 KB |
testcase_15 | AC | 1 ms
4,376 KB |
testcase_16 | AC | 1 ms
4,376 KB |
testcase_17 | AC | 2 ms
4,376 KB |
testcase_18 | AC | 679 ms
13,876 KB |
testcase_19 | AC | 14 ms
4,376 KB |
testcase_20 | AC | 560 ms
11,704 KB |
testcase_21 | AC | 966 ms
13,908 KB |
testcase_22 | AC | 139 ms
13,960 KB |
testcase_23 | AC | 3 ms
4,376 KB |
testcase_24 | AC | 25 ms
4,376 KB |
testcase_25 | AC | 12 ms
4,376 KB |
testcase_26 | AC | 180 ms
5,348 KB |
testcase_27 | AC | 242 ms
5,972 KB |
testcase_28 | AC | 485 ms
8,660 KB |
testcase_29 | AC | 191 ms
5,440 KB |
testcase_30 | AC | 576 ms
13,604 KB |
testcase_31 | AC | 138 ms
13,552 KB |
testcase_32 | AC | 2 ms
4,376 KB |
testcase_33 | AC | 11 ms
4,376 KB |
testcase_34 | AC | 40 ms
4,376 KB |
ソースコード
#include <bits/stdc++.h> #include <atcoder/modint> using namespace std; using namespace atcoder; //using mint = modint998244353; using mint = modint1000000007; using ll = long long; using ld = long double; using pll = pair<ll, ll>; using tlll = tuple<ll, ll, ll>; constexpr ll INF = 1LL << 60; template<class T> bool chmin(T& a, T b) {if (a > b) {a = b; return true;} return false;} template<class T> bool chmax(T& a, T b) {if (a < b) {a = b; return true;} return false;} ll safemod(ll A, ll M) {return (A % M + M) % M;} ll divfloor(ll A, ll B) {if (B < 0) {return divfloor(-A, -B);} return (A - safemod(A, B)) / B;} ll divceil(ll A, ll B) {if (B < 0) {return divceil(-A, -B);} return divfloor(A + B - 1, B);} #define FINALANS(A) do {cout << (A) << '\n'; exit(0);} while (false) int main() { ll N, K; cin >> N >> K; vector<ll> A(N); for (ll i = 0; i < N; i++) { cin >> A.at(i); } if (K == 1) FINALANS((mint(2).pow(N) - 1).val()); vector<ll> D; for (ll d = 1; d * d <= K; d++) { if (K % d == 0) { D.push_back(d); if (d != K / d) D.push_back(K / d); } } sort(D.begin(), D.end()); ll M = D.size(); map<ll, ll> mp; for (ll i = 0; i < M; i++) { mp[D.at(i)] = i; //cerr << D.at(i) << endl; } vector<vector<mint>> dp(N + 1, vector<mint>(M, 0)); dp.at(0).at(0) = 1; for (ll i = 0; i < N; i++) { for (ll j = 0; j < M; j++) { ll d = D.at(j); ll nd = __gcd(d * __gcd(A.at(i), K), K); ll nj = mp[nd]; dp.at(i + 1).at(nj) += dp.at(i).at(j); dp.at(i + 1).at(j) += dp.at(i).at(j); //cerr << j << " " << d << " " << A.at(i) << " " << nd << " " << nj << endl; } } mint ans = dp.at(N).at(M - 1); cout << ans.val() << endl; /* for (ll i = 0; i < N + 1; i++) { for (ll j = 0; j < M; j++) { cerr << dp.at(i).at(j).val() << " "; } cerr << endl; } //*/ }