結果
| 問題 |
No.675 ドットちゃんたち
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2019-02-03 22:05:46 |
| 言語 | Nim (2.2.0) |
| 結果 |
AC
|
| 実行時間 | 46 ms / 2,000 ms |
| コード長 | 3,274 bytes |
| コンパイル時間 | 3,160 ms |
| コンパイル使用メモリ | 77,704 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-07-01 11:15:45 |
| 合計ジャッジ時間 | 4,907 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 8 |
コンパイルメッセージ
/home/judge/data/code/Main.nim(4, 1) Warning: template 'stopwatch' is implicitly redefined; this is deprecated, add an explicit .redefine pragma [ImplicitTemplateRedefinition] /home/judge/data/code/Main.nim(1, 8) Warning: imported and not used: 'sequtils' [UnusedImport]
ソースコード
import sequtils,times
proc printf(formatstr: cstring){.header: "<stdio.h>", varargs.}
template stopwatch(body) = (let t1 = cpuTime();body;echo "TIME:",(cpuTime() - t1) * 1000,"ms")
template stopwatch(body) = body
template useMatrix =
type Matrix[T] = ref object
w,h:int
data: seq[T]
proc `[]`[T](m:Matrix[T],x,y:int):T = m.data[x + y * m.w]
proc `[]=`[T](m:var Matrix[T],x,y:int,val:T) = m.data[x + y * m.w] = val
proc newMatrix[T](w,h:int):Matrix[T] =
new(result)
result.w = w
result.h = h
result.data = newSeq[T](w * h)
proc identityMatrix[T](n:int):Matrix[T] =
result = newMatrix[T](n,n)
for i in 0..<n: result[i,i] = 1
proc newMatrix[T](arr:seq[seq[T]]):Matrix[T] =
new(result)
result.w = arr[0].len
result.h = arr.len
result.data = newSeq[T](result.w * result.h)
for x in 0..<result.w:
for y in 0..<result.h:
result[x,y] = arr[y][x]
proc `*`[T](a,b:Matrix[T]): Matrix[T] =
assert a.w == b.h
result = newMatrix[T](a.h,b.w)
for y in 0..<a.h:
for x in 0..<b.w:
var n : T
for k in 0..<a.w:
n += a[k,y] * b[x,k]
result[x,y] = n
var unsafeMat = identityMatrix[int](3)
proc `*=`[T](a:var Matrix[T],b:Matrix[T]) =
assert a.w == b.h
for y in 0..<a.h:
for x in 0..<b.w:
var n : T
for k in 0..<a.w:
n += a[k,y] * b[x,k]
unsafeMat[x,y] = n
swap(unsafeMat,a)
proc `^`[T](m:Matrix[T],n:int) : Matrix[T] =
assert m.w == m.h
if n <= 0 : return identityMatrix[T](m.w)
if n == 1 : return m
let m2 = m^(n div 2)
if n mod 2 == 0 : return m2 * m2
return m2 * m2 * m
proc `$`[T](m:Matrix[T]) : string =
result = ""
for y in 0..<m.h:
result &= "["
for x in 0..<m.w:
result &= $m[x,y]
if x != m.w - 1 : result &= ","
result &= "]"
if y != m.h - 1 : result &= "\n"
result &= "\n"
proc `*`[T](a:Matrix,b:seq[T]):seq[T] =
assert a.w == b.len
result = newSeq[T](b.len)
for x in 0..<a.h:
var n : T
for k in 0..<a.w:
n += a[k,x] * b[k]
result[x] = n
useMatrix()
proc getchar_unlocked():char {. importc:"getchar_unlocked",header: "<stdio.h>" .}
proc scan(): int =
var minus = false
block:
let k = getchar_unlocked()
if k == '-' : minus = true
else: result = 10 * result + k.ord - '0'.ord
while true:
let k = getchar_unlocked()
if k < '0' or k > '9': break
result = 10 * result + k.ord - '0'.ord
if minus: result *= -1
stopwatch:
let n = scan()
let px = scan()
let py = scan()
var C = newSeq[tuple[q,v:int]](n)
for i in 0..<n:
C[i].q = scan()
if C[i].q == 3 : continue
C[i].v = scan()
var M = identityMatrix[int](3)
var answers = newSeq[tuple[x,y:int]](n)
let m90 = newMatrix(@[@[0,1,0],@[-1,0,0],@[0,0,1]])
var my = newMatrix(@[@[1,0,0],@[0,1,9],@[0,0,1]])
var mx = newMatrix(@[@[1,0,9],@[0,1,0],@[0,0,1]])
stopwatch:
let p = @[px,py,1]
for i in (n-1).countdown(0):
let (q,v) = C[i]
if q == 3: M *= m90
elif q == 2:
my[2,1] = v
M *= my
else:
mx[2,0] = v
M *= mx
let mp = M * p
answers[i] = (mp[0],mp[1])
stopwatch:
for ans in answers: printf("%d %d\n",ans.x,ans.y)