結果

問題 No.1810 RGB Biscuits
コンテスト
ユーザー 👑 obakyan
提出日時 2022-01-15 23:04:07
言語 Lua
(LuaJit 2.1.1734355927)
結果
AC  
実行時間 69 ms / 2,000 ms
コード長 1,838 bytes
コンパイル時間 236 ms
コンパイル使用メモリ 5,248 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-22 05:02:59
合計ジャッジ時間 1,731 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

local ffi = require("ffi")
local C = ffi.C
ffi.cdef[[
long long atoll(const char*);
]]

local function lltonumber(str)
  return C.atoll(str)
end

local mod = 1000000007
local mfl = math.floor

local function bmul(x, y)
  local x1, y1 = mfl(x / 31623), mfl(y / 31623)
  local x0, y0 = x - x1 * 31623, y - y1 * 31623
  return (x1 * y1 * 14122 + (x1 * y0 + x0 * y1) * 31623 + x0 * y0) % mod
end

local function badd(x, y)
  return (x + y) % mod
end

local n = 3
local tmpmat = {}
for i = 1, n do
  tmpmat[i] = {}
end
local tmpvec = {}

local function matmat(mat2, mat)
  for i = 1, n do
    for j = 1, n do
      local v = 0
      for k = 1, n do
        v = badd(v, bmul(mat2[i][k], mat[k][j]))
      end
      tmpmat[i][j] = v
    end
  end
  for i = 1, n do for j = 1, n do
    mat[i][j] = tmpmat[i][j]
  end end
end
local function matpow(mat, p)
  local ret = {}
  local z = {}
  for i = 1, n do
    ret[i] = {}
    z[i] = {}
    for j = 1, n do
      ret[i][j] = i == j and 1 or 0
      z[i][j] = mat[i][j]
    end
  end
  while 0LL < p do
    if p % 2LL == 1LL then
      matmat(z, ret)
    end
    matmat(z, z)
    p = p / 2LL
  end
  return ret
end

local function matvec(mat, vec)
  for i = 1, n do
    local v = 0
    for k = 1, n do
      v = badd(v, bmul(mat[i][k], vec[k]))
    end
    tmpvec[i] = v
  end
  for i = 1, n do
    vec[i] = tmpvec[i]
  end
end

-- RGB
local a, b = io.read("*n", "*n")
local t1 = {
  {1, 0, 0},
  {0, 1, 0},
  {a, b, 1}
}
local t2 = {
  {0, 0, 1},
  {1, 0, 0},
  {0, 0, 0}
}
local t3 = {
  {1, 0, 0},
  {0, 1, 0},
  {a, b, 1}
}
matmat(t2, t3)

local q = io.read("*n", "*l")
for iq = 1, q do
  local v = {1, 1, 0}
  local x = lltonumber(io.read())
  local hx = x / 2LL
  local m = matpow(t3, hx)
  if x % 2LL == 1LL then
    matmat(t1, m)
  end
  matvec(m, v)
  print((v[1] + v[2] + v[3]) % mod)
end
0