結果

問題 No.1073 無限すごろく
ユーザー yakamotoyakamoto
提出日時 2020-06-05 22:10:08
言語 Kotlin
(1.9.10)
結果
AC  
実行時間 264 ms / 2,000 ms
コード長 4,329 bytes
コンパイル時間 17,749 ms
コンパイル使用メモリ 438,000 KB
実行使用メモリ 50,768 KB
最終ジャッジ日時 2023-08-22 15:48:06
合計ジャッジ時間 27,181 ms
ジャッジサーバーID
(参考情報)
judge11 / judge15
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 259 ms
50,540 KB
testcase_01 AC 259 ms
50,584 KB
testcase_02 AC 257 ms
50,584 KB
testcase_03 AC 257 ms
50,408 KB
testcase_04 AC 254 ms
50,408 KB
testcase_05 AC 257 ms
50,268 KB
testcase_06 AC 255 ms
50,356 KB
testcase_07 AC 255 ms
50,516 KB
testcase_08 AC 253 ms
50,448 KB
testcase_09 AC 255 ms
50,320 KB
testcase_10 AC 255 ms
50,488 KB
testcase_11 AC 264 ms
50,356 KB
testcase_12 AC 256 ms
50,576 KB
testcase_13 AC 257 ms
50,768 KB
testcase_14 AC 257 ms
50,324 KB
testcase_15 AC 257 ms
50,440 KB
testcase_16 AC 258 ms
50,248 KB
testcase_17 AC 259 ms
50,268 KB
testcase_18 AC 256 ms
50,348 KB
testcase_19 AC 254 ms
50,328 KB
testcase_20 AC 253 ms
50,744 KB
testcase_21 AC 253 ms
50,424 KB
testcase_22 AC 256 ms
50,328 KB
testcase_23 AC 258 ms
50,300 KB
testcase_24 AC 259 ms
50,412 KB
testcase_25 AC 256 ms
50,276 KB
testcase_26 AC 253 ms
50,340 KB
testcase_27 AC 256 ms
50,568 KB
testcase_28 AC 254 ms
50,344 KB
testcase_29 AC 257 ms
50,292 KB
testcase_30 AC 257 ms
50,336 KB
testcase_31 AC 256 ms
50,560 KB
testcase_32 AC 257 ms
50,524 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