結果

問題 No.854 公平なりんご分配
ユーザー hedwig100
提出日時 2020-05-10 12:03:35
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,792 bytes
コンパイル時間 1,651 ms
コンパイル使用メモリ 174,940 KB
実行使用メモリ 41,344 KB
最終ジャッジ日時 2024-07-07 12:25:35
合計ジャッジ時間 9,562 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 50 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 = 3; i * i <= n; i += 2) {
        if (prime[i]) {
            int j = 2;
            while (i * j <= n) {
                prime[i * j] = false;
                j ++;
            }
        }
    }
    vector<int> ans = {2};
    for (int i = 3; 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(2002);
    int Length = primes.size();
    vector<vector<int>> fact(N + 1);
    fact[0] = vector<int>(Length,0);
    rep(i,N) fact[i + 1] = factorization(A[i],primes);
    rep(i,N)rep(j,Length) {
        fact[i + 1][j] += fact[i][j];
    }
    
    int Q; cin >> Q;
    rep(i,Q) {
        int P,L,R; cin >> P >> L >> R;
        rep(j,Length) {
            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