結果

問題 No.1112 冥界の音楽
コンテスト
ユーザー yakamoto
提出日時 2020-07-10 21:42:08
言語 Kotlin
(2.1.0)
結果
AC  
実行時間 436 ms / 2,000 ms
コード長 3,744 bytes
コンパイル時間 19,826 ms
コンパイル使用メモリ 466,812 KB
実行使用メモリ 55,844 KB
最終ジャッジ日時 2024-10-11 08:33:54
合計ジャッジ時間 31,568 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 34
権限があれば一括ダウンロードができます

ソースコード

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
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
}
class Solver(stream: InputStream, private val out: java.io.PrintWriter) {
  fun solve() {
    val K = ni()
    val M = ni()
    val N = nl()
    val A = Array(K * K){LongArray(K * K)}
    val s = Array(K * K){ LongArray(1) }
    for (i in 0 until M) {
      val p = ni() - 1
      val q = ni() - 1
      val r = ni() - 1
      A[q * K + r][p * K + q]++
      if (p == 0) {
        s[q * K + r][0]++
      }
    }

    if (N == 3L) {
      out.println(s.withIndex().filter { it.index % K == 0 }.map{it.value[0]}.sum() % MOD)
      return
    }

    val B = powMat(A, N - 3, MOD)
    debug{"A"}
    debugDim(A)
    debug{"B"}
    debugDim(B)
    debug{"s"}
    debugDim(s)
    val ans = mulMat(B, s, MOD)
    out.println(ans.withIndex().filter { it.index % K == 0 }.map{it.value[0]}.sum() % MOD)
  }













































  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