結果

問題 No.854 公平なりんご分配
ユーザー hedwig100
提出日時 2020-05-10 11:29:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,833 bytes
コンパイル時間 1,622 ms
コンパイル使用メモリ 174,928 KB
実行使用メモリ 61,568 KB
最終ジャッジ日時 2024-07-07 11:20:00
合計ジャッジ時間 10,194 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 17 WA * 33 TLE * 1 -- * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i,n) for (int i = 0; i < (n); i ++)
using namespace std;
typedef long long ll;
typedef pair<ll,ll> PL;
typedef pair<int,int> P;
const ll INF = 1e9;
const ll MOD = 1e9 + 7;

vector<int> sieve(int n) {
    vector<bool> prime(n + 1,true);
    prime[0] = false;
    prime[1] = false;
    for (int i = 2; i * i <= n; i += 2) {
        if (prime[i]) {
            int j = 2;
            while (i * j <= n) {
                prime[i * j] = false;
                j ++;
            }
        }
    }
    vector<int> ans;
    for (int i = 2; i < n + 1; i += 2) {
        if (prime[i]) ans.push_back(i);
    }
    return ans; 
}

vector<int> factorization(int N,vector<int>& primes) {
    vector<int> ans(primes.size(),0);
    for (int i = 0;i < primes.size();i ++) {
        while (N%primes[i] == 0) {
            N /= primes[i]; ans[i] ++;
        }
    }
    if (N > 1) {
        primes.push_back(N);
        ans.push_back(1);
    }
    return ans;
}

int main() {
    int N; cin >> N;
    vector<int> A(N);
    rep(i,N) cin >> A[i];

    vector<int> primes = sieve(60);
    vector<vector<int>> fact(N + 1);
    rep(i,N) fact[i + 1] = factorization(A[i],primes);
    rep(i,N)rep(j,primes.size()) {
        if (j >= fact[i].size()) fact[i].push_back(0);
        if (j >= fact[i + 1].size()) fact[i + 1].push_back(0);
        fact[i + 1][j] += fact[i][j];
    }
    
    int Q; cin >> Q;
    rep(i,Q) {
        int P,L,R; cin >> P >> L >> R;
        rep(j,primes.size()) {
            int divisor = fact[R][j] - fact[L - 1][j];
            while (P%primes[j] == 0 && divisor > 0) {
                P /= primes[j];
                divisor --;
            }
            if (P%primes[j] == 0) break;
        }
        if (P > 1) cout << "NO" << endl;
        else cout << "Yes" << endl;
    }
    return 0;
}
0