結果
問題 | No.891 隣接3項間の漸化式 |
ユーザー |
👑 |
提出日時 | 2020-04-10 21:23:23 |
言語 | Lua (LuaJit 2.1.1734355927) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,353 bytes |
コンパイル時間 | 449 ms |
コンパイル使用メモリ | 5,248 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-15 16:47:27 |
合計ジャッジ時間 | 1,784 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
local bls, brs = bit.lshift, bit.rshiftlocal mod = 1000000007local mfl = math.floorlocal function bmul(x, y)local x1, y1 = mfl(x / 31623), mfl(y / 31623)local x0, y0 = x - x1 * 31623, y - y1 * 31623return (x1 * y1 * 14122 + (x1 * y0 + x0 * y1) * 31623 + x0 * y0) % modendlocal function badd(x, y)return (x + y) % modendlocal function bsub(x, y)return x < y and x - y + mod or x - yendlocal function modpow(src, pow)local res = 1while 0 < pow doif pow % 2 == 1 thenres = bmul(res, src)pow = pow - 1endsrc = bmul(src, src)pow = mfl(pow / 2)endreturn resendlocal function modinv(src)return modpow(src, mod - 2)endlocal a, b, n = io.read("*n", "*n", "*n")if n <= 1 thenprint(n)elseif a == 0 and b == 0 thenprint(0)elseif b == 0 thenprint(modpow(a, n - 1))elseif a == 0 thenif n % 2 == 0 then print(0)elseprint(modpow(b, brs((n - 1), 1)))endelselocal z, r = bmul(a, 500000004), 1local sq = badd(b, bmul(z, z))local resz, resr = 1, 0local pow = nwhile 0 < pow doif pow % 2 == 1 thenresz, resr = badd(bmul(resz, z), bmul(sq, bmul(resr, r))), badd(bmul(resz, r), bmul(resr, z))pow = pow - 1endz, r = badd(bmul(z, z), bmul(sq, bmul(r, r))), bmul(badd(z, z), r)pow = brs(pow, 1)endprint(resr)end