結果

問題 No.1073 無限すごろく
ユーザー siman
提出日時 2022-01-25 07:22:33
言語 Ruby
(3.4.1)
結果
AC  
実行時間 85 ms / 2,000 ms
コード長 783 bytes
コンパイル時間 494 ms
コンパイル使用メモリ 7,296 KB
実行使用メモリ 12,288 KB
最終ジャッジ日時 2024-12-15 15:04:56
合計ジャッジ時間 4,189 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます
コンパイルメッセージ
Syntax OK

ソースコード

diff # 

class Integer
  def mod_inverse(mod)
    self.pow(mod - 2, mod)
  end
end

MOD = 10 ** 9 + 7

def mul(a, b)
  ret = Array.new(a.size) { Array.new(b[0].size, 0) }

  a.size.times do |i|
    b.size.times do |k|
      b[0].size.times do |j|
        ret[i][j] = (ret[i][j] + a[i][k] * b[k][j]) % MOD
      end
    end
  end

  ret
end

def pow(a, n)
  h = a.size
  w = a[0].size
  ret = Array.new(h) { Array.new(w, 0) }

  h.times do |i|
    ret[i][i] = 1
  end

  while n > 0
    ret = mul(ret, a) if n % 2 == 1
    a = mul(a, a)
    n /= 2
  end

  ret
end

A = [
  [1, 1, 1, 1, 1, 1],
  [6, 0, 0, 0, 0, 0], 
  [0, 6, 0, 0, 0, 0], 
  [0, 0, 6, 0, 0, 0], 
  [0, 0, 0, 6, 0, 0], 
  [0, 0, 0, 0, 6, 0], 
]

N = gets.to_i
B = pow(A, N)

puts 6.pow(N, MOD).mod_inverse(MOD) * B[0][0] % MOD
0