結果

問題 No.219 巨大数の概算
ユーザー むらため
提出日時 2019-01-25 17:13:27
言語 Nim
(2.2.0)
結果
AC  
実行時間 5 ms / 1,500 ms
コード長 1,359 bytes
コンパイル時間 2,416 ms
コンパイル使用メモリ 61,216 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-07-01 10:47:10
合計ジャッジ時間 5,053 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 1
other AC * 51
権限があれば一括ダウンロードができます

ソースコード

diff #

import math,algorithm
template times*(n:int,body) = (for _ in 0..<n: body)

proc getchar_unlocked():char {. importc:"getchar_unlocked",header: "<stdio.h>" .}
proc scan(): int =
  while true:
    let k = getchar_unlocked()
    if k < '0': break
    result = 10 * result + k.ord - '0'.ord


template div10(a:int32) : int32 = cast[int32]((0x1999999A * cast[int64](a)) shr 32)
template put(n:int) = putchar_unlocked("0123456789"[n])
proc putchar_unlocked(c:char){. importc:"putchar_unlocked",header: "<stdio.h>" .}
proc printInt(a:int) =
  if a == 0:
    putchar_unlocked('0')
    return
  var n = a
  var rev = a
  var cnt = 0
  while rev mod 10 == 0:
    cnt += 1
    rev = rev div 10
  rev = 0
  while n != 0:
    rev = rev * 10 + n mod 10
    n = n div 10
  while rev != 0:
    putchar_unlocked((rev mod 10 + '0'.ord).chr)
    rev = rev div 10
  while cnt != 0:
    putchar_unlocked('0')
    cnt -= 1

var P = newSeq[float](101)
for i in 10..100: P[i] = log10(i.float / 10.0)
scan().times:
  let a = scan()
  let b = scan()
  let xyz = b.float * log10(a.float)
  let z = xyz.int
  let p = xyz - z.float # 10^p == A => p*log10 = logA
  let i = cast[int32](P.upperBound(p) - 1)
  let x = i.div10
  let y = i - (x * 10)
  put(x)
  putchar_unlocked(' ')
  put(y)
  putchar_unlocked(' ')
  printInt(z)
  putchar_unlocked('\n')
  # b*(log(a)) == z * log10 + log(XY)
0