結果
問題 | No.1330 Multiply or Divide |
ユーザー |
|
提出日時 | 2024-04-02 03:11:32 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,228 bytes |
コンパイル時間 | 2,173 ms |
コンパイル使用メモリ | 196,828 KB |
最終ジャッジ日時 | 2025-02-20 19:33:04 |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 43 WA * 3 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/modint>using namespace atcoder;using namespace std;using ll = long long;using mint = modint998244353;int main() {int N;ll M, P;cin >> N >> M >> P;int K = 12;vector<ll> A(K + 1, 0);bool ret = false;for (int i = 1;i <= N;i++) {ll a;cin >> a;if (a > M) {cout << 1 << endl;return 0;}int c = 0;while (a % P == 0) {c++;a /= P;}A[c] = max(A[c], a);if (a != 1) {ret = true;}}if (!ret) {cout << -1 << endl;return 0;}vector<ll> dp(K + 1, 1);ll ep = 1;for (int i = 0;i <= K;i++) {if (A[i] == 0) {if (ep > M) continue;ep *= P;continue;}vector<ll> ndp(K + 1, 1);ll x = A[i];ll y = ep * A[i];for (int j = 1;j <= K;j++) {//iをつかわないndp[j] = dp[j];if (dp[j] == M + 1) continue;//iをつかう//いっかいかけるでoverになるばあいif (y * ndp[j - 1] > M) {ndp[j] = M + 1;continue;}//ならなかったばあいif (i + 1 <= j) {ndp[j] = max(ndp[j], min(M + 1, x * ndp[j - i - 1]));}}swap(dp, ndp);ep *= P;}for (int i = 1;i <= K;i++) {if (dp[i] == M + 1) {cout << i << endl;return 0;}}}