結果
| 問題 |
No.1112 冥界の音楽
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2019-11-24 02:49:01 |
| 言語 | Java (openjdk 23) |
| 結果 |
AC
|
| 実行時間 | 121 ms / 2,000 ms |
| コード長 | 3,865 bytes |
| コンパイル時間 | 2,443 ms |
| コンパイル使用メモリ | 78,540 KB |
| 実行使用メモリ | 39,232 KB |
| 最終ジャッジ日時 | 2024-10-11 21:17:58 |
| 合計ジャッジ時間 | 6,080 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 34 |
ソースコード
// 制約侵害を確認
import java.io.FileNotFoundException;
import java.io.IOException;
import java.io.InputStream;
import java.util.Arrays;
import java.util.HashSet;
import java.util.NoSuchElementException;
public class Main {
final long MOD = (long) 1e9 + 7;
long[][] mul(long[][] a, long[][] b) {
long[][] ret = new long[a.length][b[0].length];
for (int i = 0; i < a.length; ++i) {
for (int j = 0; j < b[0].length; ++j) {
for (int k = 0; k < a[i].length; ++k) {
ret[i][j] += a[i][k] * b[k][j] % MOD;
ret[i][j] %= MOD;
}
}
}
return ret;
}
long[][] pow(long[][] a, long n) {
long[][] ret = new long[a.length][a.length];
for (int i = 0; i < a.length; ++i)
ret[i][i] = 1;
for (; n > 0; n >>= 1, a = mul(a, a)) {
if (n % 2 == 1) {
ret = mul(ret, a);
}
}
return ret;
}
// O(K^6 log(N))
void run() {
FastScanner sc = new FastScanner();
int K = sc.nextInt();
int M = sc.nextInt();
long N = sc.nextLong();
if (!(1 <= K && K <= 6 && 1 <= M && M <= K * K * K && 3 <= N && N <= 1e18))
throw new AssertionError();
int[] P = new int[M];
int[] Q = new int[M];
int[] R = new int[M];
long[][] mat = new long[K * K][K * K];
HashSet<Long> set = new HashSet<>();
for (int i = 0; i < M; ++i) {
P[i] = sc.nextInt();
Q[i] = sc.nextInt();
R[i] = sc.nextInt();
if (!(1 <= P[i] && P[i] <= K && 1 <= Q[i] && Q[i] <= K && 1 <= R[i] && R[i] <= K))
throw new AssertionError();
--P[i];
--Q[i];
--R[i];
long key = ((long) P[i] << 20) | ((long) Q[i] << 10) | ((long) R[i]);
if (set.contains(key))
throw new AssertionError();
else
set.add(key);
mat[K * R[i] + Q[i]][K * Q[i] + P[i]] = 1;
}
long[][] vec = new long[K * K][1];
for (int i = 0; i < K; ++i) {
vec[K * i + 0][0] = 1;
}
vec = mul(pow(mat, N - 2), vec);
long ans = 0;
for (int i = 0; i < K; ++i) {
ans += vec[K * 0 + i][0];
ans %= MOD;
}
System.out.println(ans);
}
public static void main(String[] args) throws FileNotFoundException {
new Main().run();
}
void tr(Object... objects) {
System.out.println(Arrays.deepToString(objects));
}
}
class FastScanner {
private final InputStream in = System.in;
private final byte[] buffer = new byte[1024];
private int ptr = 0;
private int buflen = 0;
private boolean hasNextByte() {
if (ptr < buflen) {
return true;
} else {
ptr = 0;
try {
buflen = in.read(buffer);
} catch (IOException e) {
e.printStackTrace();
}
if (buflen <= 0) {
return false;
}
}
return true;
}
private int readByte() {
if (hasNextByte())
return buffer[ptr++];
else
return -1;
}
private static boolean isPrintableChar(int c) {
return 33 <= c && c <= 126;
}
public boolean hasNext() {
while (hasNextByte() && !isPrintableChar(buffer[ptr]))
ptr++;
return hasNextByte();
}
public String next() {
if (!hasNext())
throw new NoSuchElementException();
StringBuilder sb = new StringBuilder();
int b = readByte();
while (isPrintableChar(b)) {
sb.appendCodePoint(b);
b = readByte();
}
return sb.toString();
}
public long nextLong() {
if (!hasNext())
throw new NoSuchElementException();
long n = 0;
boolean minus = false;
int b = readByte();
if (b == '-') {
minus = true;
b = readByte();
}
if (b < '0' || '9' < b) {
throw new NumberFormatException();
}
while (true) {
if ('0' <= b && b <= '9') {
n *= 10;
n += b - '0';
} else if (b == -1 || !isPrintableChar(b)) {
return minus ? -n : n;
} else {
throw new NumberFormatException();
}
b = readByte();
}
}
public int nextInt() {
long nl = nextLong();
if (nl < Integer.MIN_VALUE || nl > Integer.MAX_VALUE)
throw new NumberFormatException();
return (int) nl;
}
public double nextDouble() {
return Double.parseDouble(next());
}
}