結果

問題 No.854 公平なりんご分配
ユーザー hedwig100hedwig100
提出日時 2020-05-10 11:46:19
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 1,839 bytes
コンパイル時間 2,109 ms
コンパイル使用メモリ 174,168 KB
実行使用メモリ 55,788 KB
最終ジャッジ日時 2023-09-21 18:21:15
合計ジャッジ時間 11,218 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,384 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 2 ms
4,384 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 1 ms
4,380 KB
testcase_08 AC 1 ms
4,376 KB
testcase_09 AC 1 ms
4,376 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 2 ms
4,376 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 1 ms
4,376 KB
testcase_14 AC 1 ms
4,380 KB
testcase_15 AC 2 ms
4,376 KB
testcase_16 AC 1 ms
4,376 KB
testcase_17 AC 2 ms
4,376 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 2 ms
4,376 KB
testcase_20 AC 2 ms
4,376 KB
testcase_21 AC 1 ms
4,376 KB
testcase_22 AC 4 ms
4,376 KB
testcase_23 AC 3 ms
4,376 KB
testcase_24 AC 6 ms
4,768 KB
testcase_25 AC 3 ms
4,376 KB
testcase_26 AC 6 ms
4,700 KB
testcase_27 AC 4 ms
4,680 KB
testcase_28 AC 3 ms
4,380 KB
testcase_29 AC 3 ms
4,380 KB
testcase_30 AC 3 ms
4,376 KB
testcase_31 AC 5 ms
4,416 KB
testcase_32 AC 81 ms
42,908 KB
testcase_33 AC 75 ms
27,696 KB
testcase_34 AC 139 ms
50,188 KB
testcase_35 AC 99 ms
39,096 KB
testcase_36 AC 46 ms
5,620 KB
testcase_37 AC 76 ms
40,616 KB
testcase_38 AC 57 ms
30,240 KB
testcase_39 AC 199 ms
53,688 KB
testcase_40 AC 82 ms
23,876 KB
testcase_41 AC 97 ms
36,584 KB
testcase_42 AC 129 ms
55,788 KB
testcase_43 AC 145 ms
45,508 KB
testcase_44 AC 148 ms
45,280 KB
testcase_45 AC 121 ms
13,656 KB
testcase_46 AC 184 ms
46,592 KB
testcase_47 AC 76 ms
30,356 KB
testcase_48 AC 117 ms
46,620 KB
testcase_49 AC 114 ms
49,756 KB
testcase_50 AC 72 ms
39,760 KB
testcase_51 AC 181 ms
44,484 KB
testcase_52 TLE -
testcase_53 -- -
testcase_54 -- -
testcase_55 -- -
testcase_56 -- -
testcase_57 -- -
testcase_58 -- -
testcase_59 -- -
testcase_60 -- -
testcase_61 -- -
testcase_62 -- -
testcase_63 -- -
testcase_64 -- -
testcase_65 -- -
testcase_66 -- -
testcase_67 -- -
testcase_68 -- -
testcase_69 -- -
testcase_70 -- -
testcase_71 -- -
testcase_72 -- -
testcase_73 -- -
testcase_74 -- -
testcase_75 -- -
testcase_76 -- -
testcase_77 -- -
testcase_78 -- -
testcase_79 -- -
testcase_80 -- -
testcase_81 -- -
testcase_82 -- -
testcase_83 -- -
testcase_84 -- -
testcase_85 -- -
testcase_86 -- -
testcase_87 -- -
testcase_88 -- -
testcase_89 -- -
testcase_90 -- -
testcase_91 -- -
testcase_92 -- -
testcase_93 -- -
権限があれば一括ダウンロードができます

ソースコード

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(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