結果

問題 No.1073 無限すごろく
ユーザー yakamotoyakamoto
提出日時 2020-06-05 22:10:08
言語 Kotlin
(1.9.23)
結果
AC  
実行時間 316 ms / 2,000 ms
コード長 4,329 bytes
コンパイル時間 15,240 ms
コンパイル使用メモリ 465,840 KB
実行使用メモリ 50,188 KB
最終ジャッジ日時 2024-05-09 21:32:07
合計ジャッジ時間 25,758 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 279 ms
50,060 KB
testcase_01 AC 284 ms
50,056 KB
testcase_02 AC 286 ms
50,144 KB
testcase_03 AC 285 ms
50,156 KB
testcase_04 AC 281 ms
49,976 KB
testcase_05 AC 280 ms
50,152 KB
testcase_06 AC 279 ms
50,156 KB
testcase_07 AC 283 ms
50,188 KB
testcase_08 AC 280 ms
50,124 KB
testcase_09 AC 282 ms
50,048 KB
testcase_10 AC 281 ms
50,160 KB
testcase_11 AC 276 ms
50,004 KB
testcase_12 AC 273 ms
50,024 KB
testcase_13 AC 275 ms
50,160 KB
testcase_14 AC 279 ms
50,068 KB
testcase_15 AC 279 ms
50,044 KB
testcase_16 AC 279 ms
50,052 KB
testcase_17 AC 285 ms
50,072 KB
testcase_18 AC 275 ms
49,972 KB
testcase_19 AC 277 ms
50,040 KB
testcase_20 AC 275 ms
50,076 KB
testcase_21 AC 278 ms
50,140 KB
testcase_22 AC 276 ms
49,972 KB
testcase_23 AC 316 ms
49,992 KB
testcase_24 AC 276 ms
50,000 KB
testcase_25 AC 275 ms
50,036 KB
testcase_26 AC 279 ms
49,964 KB
testcase_27 AC 281 ms
49,928 KB
testcase_28 AC 281 ms
50,008 KB
testcase_29 AC 277 ms
50,188 KB
testcase_30 AC 275 ms
50,028 KB
testcase_31 AC 277 ms
49,916 KB
testcase_32 AC 273 ms
50,104 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.io.BufferedReader
import java.io.InputStream
import java.io.InputStreamReader
import java.util.*
import kotlin.math.abs
import kotlin.math.max
import kotlin.math.min

val MOD = 1_000_000_007
class Comb(n: Int, val mod: Int) {

  val F = LongArray(n + 1)
  val I = LongArray(n + 1)
  init {
    F[0] = 1
    for (i in 1..n) {
      F[i] = F[i - 1] * i % mod
    }
    I[n] = powMod(F[n], (mod - 2).toLong(), mod)
    for (i in n - 1 downTo 0) {
      I[i] = I[i + 1] * (i + 1) % mod
    }
  }

  private fun powMod(a: Long, n: Long, mod: Int): Long {
    if (n == 0L) return 1
    val res = powMod(a * a % mod, n / 2, mod)
    return if (n % 2 == 1L) res * a % mod else res
  }

  fun comb(n: Int, k: Int): Long {
    return F[n] * I[k] % mod * I[n - k] % mod
  }

  fun inv(x: Int): Long {
    return I[x] * F[x - 1] % mod
  }

  /**
   * nのグループからk回重複ありで選ぶ組み合わせ数
   * n - 1のしきりとkの○で考える
   */
  fun H(n: Int, k: Int) = comb(n + k - 1, k)
}
class Solver(stream: InputStream, private val out: java.io.PrintWriter) {
  fun solve() {
    val N = nl()
    val comb = Comb(10, MOD)
    val A = arrayOf(
      longArrayOf(0, 1, 0, 0, 0, 0),
      longArrayOf(0, 0, 1, 0, 0, 0),
      longArrayOf(0, 0, 0, 1, 0, 0),
      longArrayOf(0, 0, 0, 0, 1, 0),
      longArrayOf(0, 0, 0, 0, 0, 1),
      longArrayOf(comb.inv(6), comb.inv(6), comb.inv(6), comb.inv(6), comb.inv(6), comb.inv(6))
    )

    val pow = powMat(A, N, MOD)
    debugDim(pow)
    out.println(pow[5][5])
  }

  fun powMat(a: Array<LongArray>, n: Long, mod: Int): Array<LongArray> {
    if (n == 1L) return a
    val res = powMat(mulMat(a, a, mod), n / 2, mod)
    return if (n % 2 == 1L) mulMat(res, a, mod) else res
  }
  fun mulMat(a: Array<LongArray>, b: Array<LongArray>, mod: Int): Array<LongArray> {
    assert(a[0].size == b.size)
    val len = a[0].size
    val h = a.size
    val w = b[0].size
    val res = Array(h){LongArray(w)}
    val th = 7e18.toLong()
    for (i in 0 until h) {
      for (j in 0 until w) {
        var v = 0L
        for (k in 0 until len) {
          v += a[i][k] * b[k][j]
          if (v > th) v %= mod
        }
        res[i][j] = if (v >= mod) v % mod else v
      }
    }
    return res
  }










































  private val isDebug = try {
    // なんか本番でエラーでる
    System.getenv("MY_DEBUG") != null
  } catch (t: Throwable) {
    false
  }

  private var tokenizer: StringTokenizer? = null
  private val reader = BufferedReader(InputStreamReader(stream), 32768)
  private fun next(): String {
    while (tokenizer == null || !tokenizer!!.hasMoreTokens()) {
      tokenizer = StringTokenizer(reader.readLine())
    }
    return tokenizer!!.nextToken()
  }

  private fun ni() = next().toInt()
  private fun nl() = next().toLong()
  private fun ns() = next()
  private fun na(n: Int, offset: Int = 0): IntArray {
    return map(n) { ni() + offset }
  }
  private fun nal(n: Int, offset: Int = 0): LongArray {
    val res = LongArray(n)
    for (i in 0 until n) {
      res[i] = nl() + offset
    }
    return res
  }

  private fun na2(n: Int, offset: Int = 0): Array<IntArray> {
    val a  = Array(2){IntArray(n)}
    for (i in 0 until n) {
      for (e in a) {
        e[i] = ni() + offset
      }
    }
    return a
  }

  private inline fun map(n: Int, f: (Int) -> Int): IntArray {
    val res = IntArray(n)
    for (i in 0 until n) {
      res[i] = f(i)
    }
    return res
  }

  private inline fun debug(msg: () -> String) {
    if (isDebug) System.err.println(msg())
  }

  private fun debug(a: LongArray) {
    debug { a.joinToString(" ") }
  }

  private fun debug(a: IntArray) {
    debug { a.joinToString(" ") }
  }

  private fun debug(a: BooleanArray) {
    debug { a.map { if (it) 1 else 0 }.joinToString("") }
  }

  private fun debugDim(A: Array<LongArray>) {
    if (isDebug) {
      for (a in A) {
        debug(a)
      }
    }
  }
  private fun debugDim(A: Array<IntArray>) {
    if (isDebug) {
      for (a in A) {
        debug(a)
      }
    }
  }

  /**
   * 勝手にimport消されるのを防ぎたい
   */
  private fun hoge() {
    min(1, 2)
    max(1, 2)
    abs(-10)
  }
}

fun main() {
  val out = java.io.PrintWriter(System.out)
  Solver(System.`in`, out).solve()
  out.flush()
}
0