結果
問題 |
No.1112 冥界の音楽
|
ユーザー |
|
提出日時 | 2021-02-09 18:43:44 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,091 bytes |
コンパイル時間 | 1,006 ms |
コンパイル使用メモリ | 104,900 KB |
実行使用メモリ | 13,760 KB |
最終ジャッジ日時 | 2024-07-07 00:46:07 |
合計ジャッジ時間 | 5,846 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 WA * 2 TLE * 1 |
ソースコード
#include <iostream> #include <string> #include <algorithm> #include <vector> #include <iomanip> #include <cmath> #include <stdio.h> #include <queue> #include <deque> #include <cstdio> #include <set> #include <map> #include <bitset> #include <stack> #include <cctype> using namespace std; int p[220], q[220], r[220]; int main() { long long mod = 1000000007; int m_size, m; cin >> m_size >> m; m_size = m; vector<vector<long long>> mat(220); long long y; cin >> y; y -= 3; vector<vector<long long>> ans(220); for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { mat[i].emplace_back(0); } } int co = 0; for (int i = 0; i < m; i++) { cin >> p[i] >> q[i] >> r[i]; if (y == 0 && p[i] == 1 && q[i] == 1) { co++; } } if (y == 0) { cout << co << endl; return 0; } for (int i = 0; i < m; i++) { for (int j = 0; j < m; j++) { if (q[i] == p[j] && r[i] == q[j]) { mat[j][i]++; } } } for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { ans[i].emplace_back(mat[i][j]); } } y--; while (y > 0) { if ((y & 1) == 1) { vector<vector<long long>> ans_mem(220); for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { ans_mem[i].emplace_back(0); } } for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { for (int k = 0; k < m_size; k++) { ans_mem[i][j] += ans[i][k] * mat[k][j] % mod; ans_mem[i][j] %= mod; } } } ans = ans_mem; } vector<vector<long long>> mat_mem(220); for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { mat_mem[i].emplace_back(0); } } for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { for (int k = 0; k < m_size; k++) { mat_mem[i][j] += mat[i][k] * mat[k][j] % mod; mat_mem[i][j] %= mod; } } } mat = mat_mem; y >>= 1; } long long ans1 = 0; for (int i = 0; i < m; i++) {; for (int j = 0; j < m; j++) { if (r[i] == 1 && p[j] == 1) { ans1 += ans[i][j]; ans1 %= mod; } } } cout << ans1 << endl; }