結果
問題 | No.854 公平なりんご分配 |
ユーザー |
|
提出日時 | 2019-09-09 23:18:30 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,257 bytes |
コンパイル時間 | 1,699 ms |
コンパイル使用メモリ | 179,060 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-06-28 18:33:28 |
合計ジャッジ時間 | 7,964 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 11 WA * 70 RE * 11 |
ソースコード
#include <bits/stdc++.h>using namespace std;vector<pair<int, int>> factorize(int n){vector<pair<int, int>> fs;int k = 0;for (; n % 2 == 0; n /= 2, ++k);if (k > 0) fs.emplace_back(2, k);for (int p = 3; p * p <= n; p += 2){k = 0;for(; n % p == 0; n /= p, ++k);if (k > 0) fs.emplace_back(p, k);}if (n > 1) fs.emplace_back(n, 1);return fs;}int main(){cin.tie(0);ios::sync_with_stdio(false);int N;cin >> N;vector<int> As(N, 0);for (auto & a : As) cin >> a;vector<vector<int>> cum(2001, vector<int>(1, 0));vector<vector<int>> idx(2001, vector<int>(1, 0));for (int i = 0; i < N; ++i){int a = As[i];if (a == 0){cum[0].push_back(cum[0].back() + 1);idx[0].push_back(i + 1);continue;}else if (a == 1){continue;}auto factors = factorize(a);for (auto pk : factors){auto p = pk.first;auto k = pk.second;cum[p].push_back(cum[p].back() + k);idx[p].push_back(i + 1);}}int Q;cin >> Q;int P, L, R;while (Q--){cin >> P >> L >> R;auto itR = upper_bound(idx[0].begin(), idx[0].end(), R);auto itL = upper_bound(idx[0].begin(), idx[0].end(), L - 1);--itR;--itL;if (cum[0][*itR] > cum[0][*itL]){cout << "Yes\n";continue;} else if (P == 1){cout << "Yes\n";continue;}auto factors = factorize(P);bool flag = true;for (auto pk : factors){auto p = pk.first;auto k = pk.second;if (p > 2000){flag = false;break;}auto itR = upper_bound(idx[p].begin(), idx[p].end(), R);auto itL = upper_bound(idx[p].begin(), idx[p].end(), L - 1);--itR;--itL;if (cum[p][*itR] - cum[p][*itL] < k){flag = false;break;}}if (flag){cout << "Yes\n";}else{cout << "NO\n";}}return 0;}