結果

問題 No.854 公平なりんご分配
ユーザー chaemonchaemon
提出日時 2019-07-26 22:43:10
言語 Nim
(2.0.2)
結果
RE  
実行時間 -
コード長 2,283 bytes
コンパイル時間 3,598 ms
コンパイル使用メモリ 70,804 KB
実行使用メモリ 406,592 KB
最終ジャッジ日時 2024-07-02 08:34:58
合計ジャッジ時間 15,258 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,812 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 1 ms
6,940 KB
testcase_08 AC 1 ms
6,940 KB
testcase_09 AC 1 ms
6,944 KB
testcase_10 AC 1 ms
6,944 KB
testcase_11 AC 1 ms
6,944 KB
testcase_12 AC 2 ms
6,940 KB
testcase_13 AC 2 ms
6,944 KB
testcase_14 AC 2 ms
6,940 KB
testcase_15 AC 2 ms
6,940 KB
testcase_16 AC 1 ms
6,940 KB
testcase_17 AC 2 ms
6,940 KB
testcase_18 AC 1 ms
6,940 KB
testcase_19 AC 1 ms
6,940 KB
testcase_20 AC 2 ms
6,940 KB
testcase_21 AC 2 ms
6,944 KB
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
testcase_38 RE -
testcase_39 RE -
testcase_40 RE -
testcase_41 RE -
testcase_42 RE -
testcase_43 RE -
testcase_44 RE -
testcase_45 RE -
testcase_46 RE -
testcase_47 RE -
testcase_48 RE -
testcase_49 RE -
testcase_50 RE -
testcase_51 RE -
testcase_52 RE -
testcase_53 RE -
testcase_54 RE -
testcase_55 RE -
testcase_56 RE -
testcase_57 RE -
testcase_58 RE -
testcase_59 RE -
testcase_60 RE -
testcase_61 RE -
testcase_62 RE -
testcase_63 RE -
testcase_64 RE -
testcase_65 RE -
testcase_66 RE -
testcase_67 RE -
testcase_68 RE -
testcase_69 RE -
testcase_70 RE -
testcase_71 RE -
testcase_72 RE -
testcase_73 RE -
testcase_74 RE -
testcase_75 RE -
testcase_76 RE -
testcase_77 RE -
testcase_78 RE -
testcase_79 RE -
testcase_80 RE -
testcase_81 RE -
testcase_82 MLE -
testcase_83 MLE -
testcase_84 MLE -
testcase_85 MLE -
testcase_86 MLE -
testcase_87 MLE -
testcase_88 MLE -
testcase_89 MLE -
testcase_90 MLE -
testcase_91 MLE -
testcase_92 MLE -
testcase_93 MLE -
権限があれば一括ダウンロードができます
コンパイルメッセージ
/home/judge/data/code/Main.nim(2, 57) Warning: imported and not used: 'strutils' [UnusedImport]
/home/judge/data/code/Main.nim(2, 45) Warning: imported and not used: 'math' [UnusedImport]
/home/judge/data/code/Main.nim(2, 29) Warning: imported and not used: 'tables' [UnusedImport]
/home/judge/data/code/Main.nim(2, 37) Warning: imported and not used: 'macros' [UnusedImport]

ソースコード

diff # 

#{{{ header
import algorithm, sequtils, tables, macros, math, sets, strutils
when defined(MYDEBUG):
  import header

proc scanf(formatstr: cstring){.header: "<stdio.h>", varargs.}
proc getchar(): char {.header: "<stdio.h>", varargs.}
proc nextInt(): int = scanf("%lld",addr result)
proc nextFloat(): float = scanf("%lf",addr result)
proc nextString(): string =
  var get = false
  result = ""
  while true:
    var c = getchar()
    if int(c) > int(' '):
      get = true
      result.add(c)
    else:
      if get: break
      get = false
template `max=`*(x,y:typed):void = x = max(x,y)
template `min=`*(x,y:typed):void = x = min(x,y)
template infty(T): untyped = ((T(1) shl T(sizeof(T)*8-2)) - 1)
#}}}

var pdiv: seq[int]

proc sieve_of_eratosthenes(n:int):void =
  pdiv.setLen(n)
  for i in 2..<n:
    pdiv[i] = i
  for i in 2..<n:
    if i * i >= n: break
    if pdiv[i] == i:
      for j in countup(i*i,n-1,i):
        pdiv[j] = i

proc is_prime(n:int):bool =
  return n!=1 and pdiv[n]==n

proc main():void =
  sieve_of_eratosthenes(2000)
  var
    N = nextInt()
    A = newSeqWith(N,nextInt())
    pid = newSeq[int](2000)
    zero_ct = newSeqWith(N+1,0)
    ct = 0
  for p in 2..<2000:
    if pdiv[p]==p:
      pid[p] = ct
      ct += 1
  var
    t = newSeqWith(N+1,newSeq[int](ct))
  zero_ct[0] = 0
  for i in 0..<N:
    var a = A[i]
    zero_ct[i+1] = zero_ct[i]
    for j in 0..<ct:
      t[i+1][j] = t[i][j]
    if a == 0:
      zero_ct[i] += 1
    else:
      for p in 2..a:
        if p*p > a: break
        while a mod p==0:
          a = a div p
          t[i+1][pid[p]] += 1
      if a > 1:
        t[i+1][pid[a]] += 1
  var
    Q = nextInt()
  for _ in 0..<Q:
    var
      P = nextInt()
      L = nextInt()
      R = nextInt()
      tt = newSeq[int](ct)
      u = newSeq[int](ct)
    L -= 1
    if zero_ct[R] - zero_ct[L] > 0:
      echo "Yes"
    else:
      for j in 0..<ct:
        tt[j] += t[R][j] - t[L][j]
      var a = P
      for p in 2..a:
        if p*p > a: break
        while a mod p==0:
          a = a div p
          u[pid[p]] += 1
      if a > 1:
        u[pid[a]] += 1
      var valid = true
      for i in 0..<ct:
        if tt[i] < u[i]:
          valid = false
      if valid:
        echo "Yes"
      else:
        echo "NO"
  discard


main()
0