結果
問題 | No.1704 Many Bus Stops (easy) |
ユーザー |
![]() |
提出日時 | 2024-04-13 13:45:07 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 177 ms / 2,000 ms |
コード長 | 1,181 bytes |
コンパイル時間 | 2,658 ms |
コンパイル使用メモリ | 205,096 KB |
最終ジャッジ日時 | 2025-02-21 01:01:23 |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 41 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/modint>using namespace std;using mint = atcoder::modint1000000007;void fast_io() {ios_base::sync_with_stdio(false);cin.tie(nullptr);}using vec = vector<mint>;using mat = vector<vec>;mint inv3 = mint(3).inv();mat A = {{inv3, 0, 0, 0, inv3, inv3}, {0, inv3, 0, inv3, 0, inv3},{0, 0, inv3, inv3, inv3, 0}, {1, 0, 0, 0, 0, 0},{0, 1, 0, 0, 0, 0}, {0, 0, 1, 0, 0, 0},};mat mul(mat a, mat b) {mat c(6, vec(6));for (int i = 0; i < 6; i++) {for (int j = 0; j < 6; j++) {for (int k = 0; k < 6; k++) {c[i][k] += a[i][j] * b[j][k];}}}return c;}mat pow(mat a, long long n) {mat b(6, vec(6));for (int i = 0; i < 6; i++) {b[i][i] = 1;}while (n > 0) {if (n & 1) {b = mul(b, a);}a = mul(a, a);n >>= 1;}return b;}int main() {fast_io();int t;cin >> t;for (; t--;) {int n;cin >> n;mat B = pow(A, n - 1);mint ans = B[0][0] * inv3 + B[0][3];cout << ans.val() << '\n';}}