結果

問題 No.854 公平なりんご分配
ユーザー rpy3cpp
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0