結果
| 問題 | No.1112 冥界の音楽 |
| ユーザー |
 もりを
|
| 提出日時 | 2020-07-10 22:51:09 |
| 言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
| 結果 |
TLE
(最新)
AC
(最初)
|
| 実行時間 | - |
| コード長 | 11,543 bytes |
| コンパイル時間 | 4,151 ms |
| コンパイル使用メモリ | 243,208 KB |
| 実行使用メモリ | 39,956 KB |
| 最終ジャッジ日時 | 2024-11-15 04:55:01 |
| 合計ジャッジ時間 | 22,433 ms |
|
ジャッジサーバーID （参考情報） |
judge5 / judge2 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | AC | 417 ms
39,808 KB |
| testcase_01 | AC | 22 ms
39,936 KB |
| testcase_02 | AC | 21 ms
39,808 KB |
| testcase_03 | AC | 21 ms
39,808 KB |
| testcase_04 | AC | 21 ms
39,936 KB |
| testcase_05 | AC | 21 ms
39,808 KB |
| testcase_06 | TLE | - |
| testcase_07 | AC | 22 ms
39,848 KB |
| testcase_08 | AC | 26 ms
39,920 KB |
| testcase_09 | AC | 439 ms
39,936 KB |
| testcase_10 | AC | 21 ms
39,916 KB |
| testcase_11 | AC | 431 ms
39,860 KB |
| testcase_12 | AC | 21 ms
39,936 KB |
| testcase_13 | AC | 1,786 ms
39,808 KB |
| testcase_14 | TLE | - |
| testcase_15 | AC | 443 ms
39,860 KB |
| testcase_16 | AC | 21 ms
39,808 KB |
| testcase_17 | AC | 459 ms
39,872 KB |
| testcase_18 | AC | 459 ms
39,900 KB |
| testcase_19 | AC | 21 ms
39,832 KB |
| testcase_20 | AC | 81 ms
39,852 KB |
| testcase_21 | TLE | - |
| testcase_22 | TLE | - |
| testcase_23 | AC | 447 ms
39,852 KB |
| testcase_24 | AC | 21 ms
39,900 KB |
| testcase_25 | AC | 22 ms
39,924 KB |
| testcase_26 | AC | 21 ms
39,860 KB |
| testcase_27 | AC | 22 ms
39,748 KB |
| testcase_28 | AC | 26 ms
39,936 KB |
| testcase_29 | AC | 21 ms
39,936 KB |
| testcase_30 | AC | 457 ms
39,808 KB |
| testcase_31 | AC | 27 ms
39,896 KB |
| testcase_32 | AC | 451 ms
39,956 KB |
| testcase_33 | AC | 21 ms
39,804 KB |
| testcase_34 | AC | 21 ms
39,936 KB |
| testcase_35 | AC | 21 ms
39,804 KB |
| testcase_36 | AC | 1,783 ms
39,772 KB |
ソースコード
#pragma GCC target ("avx")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx")
#include <bits/stdc++.h>
using namespace std;
// #define int long long
// #define endl '\n'
#pragma region TEMPLATE
/* TYPE */
typedef long long ll; typedef long double ld;
typedef pair<int, int> pii; typedef pair<ll, ll> pll;
typedef vector<pii> vpii; typedef vector<pll> vpll;
typedef vector<int> vi; typedef vector<ll> vl;
typedef vector<string> vst; typedef vector<bool> vb;
typedef vector<ld> vld; typedef vector<vector<int>> vvi;
template<typename T, typename Cmp = less<>> using prique = priority_queue<T, vector<T>, Cmp>;
template<typename T> using prique_r = prique<T, greater<>>;
/* CONSTANT */
#define ln '\n'
const int INF = 1 << 30; const ll INFF = 1LL << 60; const string ALPHABET = "ABCDEFGHIJKLMNOPQRSTUVWXYZ";
const int MOD = 1e9 + 7; const int MODD = 998244353; const string alphabet = "abcdefghijklmnopqrstuvwxyz";
const double EPS = 1e-9; const ld PI = 3.14159265358979323846264338327950288;
const int dx[] = { 1, 0, -1, 0, 1, -1, -1, 1, 0 };
const int dy[] = { 0, 1, 0, -1, -1, -1, 1, 1, 0 };
/* CONTAINER */
#define PB emplace_back
#define ALL(v) (v).begin(), (v).end()
#define RALL(v) (v).rbegin(), (v).rend()
#define SORT(v) sort(ALL(v))
#define RSORT(v) sort(RALL(v))
#define LESS(x, val) (lower_bound(x.begin(), x.end(), val) - x.begin())
#define LEQ(x, val) (upper_bound(x.begin(), x.end(), val) - x.begin())
#define GREATER(x, val) (int)(x).size() - LEQ((x), (val))
#define GEQ(x, val) (int)(x).size() - LESS((x), (val))
#define UNIQUE(v) sort(ALL(v)); (v).erase(unique(ALL(v)), (v).end())
template<typename T> vector<T> make_v(size_t a) { return vector<T>(a); }
template<typename T, typename... Ts> auto make_v(size_t a, Ts... ts) { return vector<decltype(make_v<T>(ts...))>(a, make_v<T>(ts...)); }
template<typename T, typename U, typename... V> enable_if_t<is_same<T, U>::value != 0> fill_v(U &u, const V... v) { u = U(v...); }
template<typename T, typename U, typename... V> enable_if_t<is_same<T, U>::value == 0> fill_v(U &u, const V... v) { for (auto &e : u) fill_v<T>(e, v...); }
/* LOOP */
#define _overload3(_1, _2, _3, name, ...) name
#define _REP(i, n) REPI(i, 0, n)
#define REPI(i, a, b) for (ll i = (ll)a; i < (ll)b; ++i)
#define REP(...) _overload3(__VA_ARGS__, REPI, _REP,)(__VA_ARGS__)
#define _RREP(i, n) RREPI(i, n, 0)
#define RREPI(i, a, b) for (ll i = (ll)a; i >= (ll)b; --i)
#define RREP(...) _overload3(__VA_ARGS__, RREPI, _RREP,)(__VA_ARGS__)
#define EACH(e, v) for (auto& e : v)
#define PERM(v) sort(ALL(v)); for (bool c##p = true; c##p; c##p = next_permutation(ALL(v)))
/* INPUT */
template<typename T> void SSS(T& t) { cin >> t; }
template<typename Head, typename... Tail> void SSS(Head&& head, Tail&&... tail) { cin >> head; SSS(tail...); }
#define SS(T, ...) T __VA_ARGS__; SSS(__VA_ARGS__);
#define SV(T, v, n) vector<T> v(n); for (auto& i : v) cin >> i;
#define SVV(T, v, n, m) vector<vector<T>> v(n, vector<T>(m)); for (auto& r : v) for (auto& i : r) cin >> i;
/* OUTPUT */
// Yes / No
inline int YES(bool x) { cout << (x ? "YES" : "NO") << endl; return 0; }
inline int Yes(bool x) { cout << (x ? "Yes" : "No") << endl; return 0; }
inline int yes(bool x) { cout << (x ? "yes" : "no") << endl; return 0; }
inline int yES(bool x) { cout << (x ? "yES" : "nO") << endl; return 0; }
inline int Yay(bool x) { cout << (x ? "Yay!" : ":(") << endl; return 0; }
// PROTOTYPE DECLARATION
template<typename T, typename U> ostream &operator<<(ostream &os, const pair<T, U> &j);
template<typename... T> ostream &operator<<(ostream &os, const tuple<T...> &t);
template<class C, enable_if_t<!is_same<C, string>::value, decltype(declval<const C &>().begin(), nullptr)> = nullptr> ostream& operator<<(ostream &os, const C &c);
template<typename T> ostream &operator<<(ostream &os, const stack<T> &j);
template<typename T> ostream &operator<<(ostream &os, const queue<T> &j);
template<typename T, typename C, typename Cmp> ostream &operator<<(ostream &os, const priority_queue<T, C, Cmp> &j);
// IMPLEMENTATION
template<typename T, typename U> ostream &operator<<(ostream &os, const pair<T, U> &j) { return os << '{' << j.first << ", " << j.second << '}'; }
template<size_t num = 0, typename... T> enable_if_t<num == sizeof...(T)> PRINT_TUPLE(ostream &os, const tuple<T...> &t) {}
template<size_t num = 0, typename... T> enable_if_t<num < sizeof...(T)> PRINT_TUPLE(ostream &os, const tuple<T...> &t) { os << get<num>(t); if (num + 1 < sizeof...(T)) os << ", "; PRINT_TUPLE<num + 1>(os, t); }
template<typename... T> ostream &operator<<(ostream &os, const tuple<T...> &t) { PRINT_TUPLE(os << '{', t); return os << '}'; }
template<class C, enable_if_t<!is_same<C, string>::value, decltype(declval<const C &>().begin(), nullptr)>> ostream& operator<<(ostream &os, const C &c) { os << '{'; for (auto it = begin(c); it != end(c); it++) { if (begin(c) != it) os << ", "; os << *it; } return os << '}'; }
template<typename T> ostream &operator<<(ostream &os, const stack<T> &j) { deque<T> d; for (auto c = j; !c.empty(); c.pop()) d.push_front(c.top()); return os << d; }
template<typename T> ostream &operator<<(ostream &os, const queue<T> &j) { deque<T> d; for (auto c = j; !c.empty(); c.pop()) d.push_back(c.front()); return os << d; }
template<typename T, typename C, typename Cmp> ostream &operator<<(ostream &os, const priority_queue<T, C, Cmp> &j) { deque<T> d; for (auto c = j; !c.empty(); c.pop()) d.push_front(c.top()); return os << d; }
// OUTPUT FUNCTION
template<typename T> int PV(T &v) { int sz = v.size(); for (int i = 0; i < sz; ++i) cout << v[i] << " \n"[i == sz - 1]; return 0; }
inline int print() { cout << endl; return 0; }
template<typename Head> int print(Head&& head){ cout << head; return print(); }
template<typename Head, typename... Tail> int print(Head&& head, Tail&&... tail) { cout << head << " "; return print(forward<Tail>(tail)...); }
#ifdef LOCAL
inline void dump() { cerr << endl; }
template<typename Head> void dump(Head&& head) { cerr << head; dump(); }
template<typename Head, typename... Tail> void dump(Head&& head, Tail&&... tail) { cerr << head << ", "; dump(forward<Tail>(tail)...); }
#define debug(...) do {cerr << __LINE__ << ":\t" << #__VA_ARGS__ << " = "; dump(__VA_ARGS__); } while (false)
#else
#define dump(...)
#define debug(...)
#endif
/* OTHER */
#define fi first
#define se second
#define MP make_pair
#define MT make_tuple
#define tmax(x, y, z) max((x), max((y), (z)))
#define tmin(x, y, z) min((x), min((y), (z)))
template<typename T, typename A, typename B> inline bool between(T x, A a, B b) { return ((a <= x) && (x < b)); }
template<typename A, typename B> inline bool chmax(A &a, const B &b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A &a, const B &b) { if (a > b) { a = b; return true; } return false; }
inline ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }
inline ll lcm(ll a, ll b) { return a / gcd(a, b) * b; }
inline ll POW(ll a, ll b) { ll r = 1; do { if (b & 1) r *= a; a *= a; } while (b >>= 1); return r; }
struct abracadabra {
abracadabra() {
cin.tie(nullptr); ios::sync_with_stdio(false);
cout << fixed << setprecision(20);
cerr << fixed << setprecision(5);
};
} ABRACADABRA;
#pragma endregion
#pragma region math modint
/**
* @brief ModInt
* @docs docs/math/modint.md
*/
template< int MODULO > struct ModInt {
using i32 = int;
using i64 = long long;
using u32 = unsigned int;
using u64 = unsigned long long;
u64 x; ModInt() : x(0) {}
ModInt(i64 y) : x(set(y % MODULO + MODULO)) {}
static u64 set(const i64 &y) { return (y < MODULO) ? y : y - MODULO; }
ModInt operator+(const ModInt &m) const { return ModInt(set(x + m.x)); }
ModInt operator-(const ModInt &m) const { return ModInt(set(x + MODULO - m.x)); }
ModInt operator*(const ModInt &m) const { return ModInt(x * m.x % MODULO); }
ModInt operator/(const ModInt &m) const { return ModInt(x) * ~ModInt(m.x); }
ModInt &operator+=(const ModInt &m) { return *this = *this + m; }
ModInt &operator-=(const ModInt &m) { return *this = *this - m; }
ModInt &operator*=(const ModInt &m) { return *this = *this * m; }
ModInt &operator/=(const ModInt &m) { return *this = *this / m; }
ModInt &operator^=(const u64 &y) { return *this = *this ^ y; }
ModInt operator~ () const { return *this ^ (MODULO - 2); }
ModInt operator- () const { return ModInt(set(MODULO - x)); }
ModInt operator++() { return *this = ModInt(set(x + 1)); }
ModInt operator--() { return *this = ModInt(set(x + MODULO - 1)); }
bool operator==(const ModInt &m) const { return x == m.x; }
bool operator!=(const ModInt &m) const { return x != m.x; }
bool operator< (const ModInt &m) const { return x < m.x; }
bool operator<=(const ModInt &m) const { return x <= m.x; }
bool operator> (const ModInt &m) const { return x > m.x; }
bool operator>=(const ModInt &m) const { return x >= m.x; }
explicit operator u64() const { return x; }
ModInt operator^(i64 y) const { return pow(x, y); }
static ModInt pow(i64 x, i64 y) {
bool neg = false;
if (y < 0) y = -y, neg = true;
ModInt u(1), t(x);
while (y) { if (y & 1) u *= t; t *= t; y >>= 1; }
return neg ? ModInt(1) / u : u;
}
friend ostream &operator<<(ostream &os, const ModInt< MODULO > &m) { return os << m.x; }
friend istream &operator>>(istream &is, ModInt< MODULO > &m) { u64 y; is >> y; m = ModInt(y); return is; }
};
constexpr int MODULO = (int)1e9 + 7;
using modint = ModInt< MODULO >;
#pragma endregion
using mint = modint;
mint dp[100][6][6][6][6][6][6];
mint init[6][6][6], nxt[6][6][6];
bool ok[6][6][6];
signed main() {
SS(int, K, M);
SS(ll, N);
vector<tuple<int, int, int>> frm, to;
REP(i, M) {
SS(int, P, Q, R);
--P, --Q, --R;
ok[P][Q][R] = true;
if (P == 0) {
frm.emplace_back(P, Q, R);
init[P][Q][R] = 1;
}
if (R == 0) to.emplace_back(P, Q, R);
}
// 2^0個すすめるとき
REP(i, K) REP(j, K) REP(k, K) {
if (not ok[i][j][k]) continue;
REP(l, K) {
if (not ok[j][k][l]) continue;
dp[0][i][j][k][j][k][l] = 1;
}
}
REP(m, 1, 61) {
// 2^m個進めるとき
// -> 2^(m-1)進めて, 2^(m-1)進める
REP(p1, K) REP(q1, K) REP(r1, K) {
REP(p2, K) REP(q2, K) REP(r2, K) {
REP(p3, K) REP(q3, K) REP(r3, K) {
dp[m][p1][q1][r1][p3][q3][r3] += dp[m - 1][p1][q1][r1][p2][q2][r2] * dp[m - 1][p2][q2][r2][p3][q3][r3];
}
}
}
}
N -= 3;
RREP(i, 60) if (N >> i & 1) {
REP(p, 6) REP(q, 6) REP(r, 6) nxt[p][q][r] = 0;
REP(p1, K) REP(q1, K) REP(r1, K) {
REP(p2, K) REP(q2, K) REP(r2, K) {
nxt[p2][q2][r2] += dp[i][p1][q1][r1][p2][q2][r2] * init[p1][q1][r1];
}
}
REP(p, 6) REP(q, 6) REP(r, 6) init[p][q][r] = nxt[p][q][r];
}
mint res = 0;
for (auto& [p, q, r]: to) {
res += init[p][q][r];
}
// REP(p, 6) REP(q, 6) REP(r, 6) debug(init[p][q][r]);
print(res);
}