結果

問題 No.854 公平なりんご分配
ユーザー chaemon
提出日時 2019-07-26 22:23:22
言語 Nim
(2.2.0)
結果
RE  
実行時間 -
コード長 2,005 bytes
コンパイル時間 3,843 ms
コンパイル使用メモリ 71,864 KB
実行使用メモリ 412,176 KB
最終ジャッジ日時 2024-07-02 07:44:28
合計ジャッジ時間 13,529 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 20 RE * 60 MLE * 2 -- * 10
権限があれば一括ダウンロードができます
コンパイルメッセージ
/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)
    ct = 0
  for p in 2..<2000:
    if pdiv[p]==p:
      pid[p] = ct
      ct += 1
  var
    t = newSeqWith(N,newSeq[int](ct))
  for i in 0..<N:
    var a = A[i]
    for p in 2..a:
      if p*p > a: break
      while a mod p==0:
        a = a div p
        t[i][pid[p]] += 1
    if a > 1:
      t[i][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
    for i in L..<R:
      for j in 0..<ct:
        tt[j] += t[i][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