結果

問題 No.854 公平なりんご分配
ユーザー kibuna
提出日時 2019-07-26 22:25:13
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 2,162 bytes
コンパイル時間 1,864 ms
コンパイル使用メモリ 180,516 KB
実行使用メモリ 813,964 KB
最終ジャッジ日時 2024-07-02 07:47:54
合計ジャッジ時間 47,336 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 34 WA * 7 RE * 29 MLE * 11 -- * 11
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
using ll = long long;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
using vi = vector<int>;
using vl = vector<ll>;
using vvi = vector<vi>;
using vvl = vector<vl>;
const ll INF = 1LL << 60;
const ll MOD = 1000000007;
template <class T>
bool chmax(T &a, const T &b) {
    return (a < b) ? (a = b, 1) : 0;
}
template <class T>
bool chmin(T &a, const T &b) {
    return (b < a) ? (a = b, 1) : 0;
}

template <typename T = int>
void primeFactors(T a, map<T, int> &facs) {
    double sqrtA = sqrt(a);
    for (int i = 2; i <= sqrtA + 1e-10; ++i) {
        while (a % i == 0) {
            facs[i]++;
            a /= i;
        }
    }
    if (a > sqrtA)
        facs[a]++;
    return;
}

template <typename T = int>
void primeFactorsP(T a, map<T, int> &facs) {
    for (int i = 2; i <= 2000; ++i) { // don't check p > 2000
        while (a % i == 0) {
            facs[i]++;
            a /= i;
        }
    }
    if (a > 1)
        facs[a]++;
    return;
}

int main() {
    int n;
    cin >> n;
    vi a(n);
    for (int i = 0; i < n; ++i) {
        cin >> a[i];
    }
    int q;
    cin >> q;
    vi p(n), l(n), r(n);
    for (int i = 0; i < q; ++i) {
        cin >> p[i] >> l[i] >> r[i];
        l[i]--;
        r[i]--;
    }

    vector<map<int, int>> facs(n);
    for (int i = 0; i < n; ++i) {
        primeFactors(a[i], facs[i]);
    }
    for (int i = 0; i < n - 1; ++i) {
        for (auto &ff : facs[i]) {
            facs[i + 1][ff.first] += ff.second;
        }
    }
    for (int i = 0; i < q; ++i) {
        map<int, int> facsToUse = facs[r[i]];
        if (l[i] > 0) {
            for (auto &ff : facs[l[i] - 1]) {
                facsToUse[ff.first] -= ff.second;
            }
        }
        map<int, int> pprimes;
        primeFactorsP(p[i], pprimes);
        bool ok = true;
        for (auto &pp : pprimes) {
            if (facsToUse[pp.first] < pp.second) {
                ok = false;
                break;
            }
        }
        if (ok)
            cout << "Yes"
                 << "\n";
        else
            cout << "NO"
                 << "\n";
    }

    return 0;
}
0