結果

問題 No.854 公平なりんご分配
ユーザー tnakao0123tnakao0123
提出日時 2019-08-03 13:39:10
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 423 ms / 3,153 ms
コード長 2,485 bytes
コンパイル時間 1,574 ms
コンパイル使用メモリ 84,916 KB
実行使用メモリ 125,084 KB
最終ジャッジ日時 2023-09-19 05:13:20
合計ジャッジ時間 17,482 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,384 KB
testcase_01 AC 2 ms
4,384 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 2 ms
4,380 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 2 ms
4,384 KB
testcase_12 AC 2 ms
4,376 KB
testcase_13 AC 2 ms
4,380 KB
testcase_14 AC 2 ms
4,380 KB
testcase_15 AC 2 ms
4,376 KB
testcase_16 AC 2 ms
4,376 KB
testcase_17 AC 2 ms
4,376 KB
testcase_18 AC 2 ms
4,380 KB
testcase_19 AC 2 ms
4,380 KB
testcase_20 AC 2 ms
4,380 KB
testcase_21 AC 2 ms
4,384 KB
testcase_22 AC 3 ms
4,380 KB
testcase_23 AC 4 ms
4,380 KB
testcase_24 AC 5 ms
4,616 KB
testcase_25 AC 4 ms
4,380 KB
testcase_26 AC 5 ms
4,660 KB
testcase_27 AC 4 ms
4,652 KB
testcase_28 AC 4 ms
4,384 KB
testcase_29 AC 2 ms
4,376 KB
testcase_30 AC 3 ms
4,376 KB
testcase_31 AC 5 ms
4,488 KB
testcase_32 AC 49 ms
29,728 KB
testcase_33 AC 46 ms
17,416 KB
testcase_34 AC 87 ms
39,988 KB
testcase_35 AC 63 ms
31,760 KB
testcase_36 AC 32 ms
5,172 KB
testcase_37 AC 46 ms
27,792 KB
testcase_38 AC 37 ms
23,612 KB
testcase_39 AC 124 ms
41,060 KB
testcase_40 AC 50 ms
15,364 KB
testcase_41 AC 61 ms
25,764 KB
testcase_42 AC 91 ms
40,184 KB
testcase_43 AC 86 ms
27,672 KB
testcase_44 AC 92 ms
35,944 KB
testcase_45 AC 74 ms
11,228 KB
testcase_46 AC 114 ms
35,940 KB
testcase_47 AC 46 ms
21,556 KB
testcase_48 AC 74 ms
38,136 KB
testcase_49 AC 71 ms
38,628 KB
testcase_50 AC 44 ms
27,416 KB
testcase_51 AC 115 ms
38,128 KB
testcase_52 AC 61 ms
19,496 KB
testcase_53 AC 37 ms
19,512 KB
testcase_54 AC 63 ms
21,556 KB
testcase_55 AC 29 ms
19,436 KB
testcase_56 AC 29 ms
19,584 KB
testcase_57 AC 37 ms
17,376 KB
testcase_58 AC 62 ms
9,300 KB
testcase_59 AC 23 ms
15,636 KB
testcase_60 AC 45 ms
17,376 KB
testcase_61 AC 19 ms
7,132 KB
testcase_62 AC 52 ms
15,320 KB
testcase_63 AC 35 ms
15,352 KB
testcase_64 AC 17 ms
9,212 KB
testcase_65 AC 57 ms
35,948 KB
testcase_66 AC 31 ms
11,216 KB
testcase_67 AC 53 ms
23,692 KB
testcase_68 AC 44 ms
13,336 KB
testcase_69 AC 60 ms
41,188 KB
testcase_70 AC 25 ms
15,428 KB
testcase_71 AC 29 ms
16,440 KB
testcase_72 AC 41 ms
5,268 KB
testcase_73 AC 44 ms
28,188 KB
testcase_74 AC 76 ms
33,928 KB
testcase_75 AC 41 ms
17,664 KB
testcase_76 AC 52 ms
27,736 KB
testcase_77 AC 67 ms
38,132 KB
testcase_78 AC 65 ms
12,684 KB
testcase_79 AC 62 ms
23,584 KB
testcase_80 AC 63 ms
24,740 KB
testcase_81 AC 48 ms
25,836 KB
testcase_82 AC 276 ms
124,912 KB
testcase_83 AC 391 ms
124,992 KB
testcase_84 AC 389 ms
125,084 KB
testcase_85 AC 414 ms
124,940 KB
testcase_86 AC 83 ms
125,016 KB
testcase_87 AC 231 ms
124,940 KB
testcase_88 AC 247 ms
124,964 KB
testcase_89 AC 226 ms
124,992 KB
testcase_90 AC 244 ms
124,960 KB
testcase_91 AC 219 ms
125,080 KB
testcase_92 AC 423 ms
125,028 KB
testcase_93 AC 407 ms
124,912 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

/* -*- coding: utf-8 -*-
 *
 * 854.cc:  No.854 公平なりんご分配 - yukicoder
 */

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<string>
#include<vector>
#include<map>
#include<set>
#include<stack>
#include<list>
#include<queue>
#include<deque>
#include<algorithm>
#include<numeric>
#include<utility>
#include<complex>
#include<functional>
 
using namespace std;

/* constant */

const int MAX_N = 100000;
const int MAX_E2 = 1 << 18; // 262144
const int MAX_P = 2000;
const int PN = 310;

/* typedef */

typedef int vec[PN];

/* global variables */

bool primes[MAX_P + 1];
vec pnums, afss[MAX_N + 1];
int zss[MAX_N + 1];

/* subroutines */

int gen_primes(int maxp, vec pnums) {
  memset(primes, true, sizeof(primes));
  primes[0] = primes[1] = false;

  int p, pn = 0;
  for (p = 2; p * p <= maxp; p++)
    if (primes[p]) {
      pnums[pn++] = p;
      for (int q = p * p; q <= maxp; q += p) primes[q] = false;
    }
  for (; p <= maxp; p++)
    if (primes[p]) pnums[pn++] = p;
  return pn;
}

bool prime_decomp(int n, int pn, vec pnums, vec fs) {
  for (int i = 0; i < pn; i++) {
    int &pi = pnums[i];
    fs[i] = 0;
    while (n % pi == 0) n /= pi, fs[i]++;
  }
  return (n == 1);
}

void copyvec(int n, vec c, vec a) { copy(a, a + n, c); }

void addvec(int n, vec c, vec a, vec b) {
  for (int i = 0; i < n; i++) c[i] = a[i] + b[i];
}

void subvec(int n, vec c, vec a, vec b) {
  for (int i = 0; i < n; i++) c[i] = a[i] - b[i];
}

bool cmpvec(int n, vec a, vec b) {
  bool ok = true;
  for (int i = 0; ok && i < n; i++)
    ok = (a[i] >= b[i]);
  return ok;
}

/* main */

int main() {
  int pn = gen_primes(MAX_P, pnums);
  //printf("pn=%d\n", pn); return 0;
  
  int n;
  scanf("%d", &n);

  memset(afss[0], 0, sizeof(vec));
  zss[0] = 0;
  
  for (int i = 0; i < n; i++) {
    int ai;
    scanf("%d", &ai);
    if (ai == 0) {
      copyvec(pn, afss[i + 1], afss[i]);
      zss[i + 1] = zss[i] + 1;
    }
    else {
      vec fs;
      prime_decomp(ai, pn, pnums, fs);
      addvec(pn, afss[i + 1], afss[i], fs);
      zss[i + 1] = zss[i];
    }
  }

  int q;
  scanf("%d", &q);

  while (q--) {
    int p, l, r;
    scanf("%d%d%d", &p, &l, &r);
    l--;

    if (zss[r] > zss[l]) puts("Yes");
    else {
      vec fs, pfs;
      if (! prime_decomp(p, pn, pnums, pfs))
	puts("NO");
      else {
	subvec(pn, fs, afss[r], afss[l]);
	if (cmpvec(pn, fs, pfs)) puts("Yes");
	else puts("NO");
      }
    }
  }
  return 0;
}
0