結果
問題 | No.1879 How many matchings? |
ユーザー |
![]() |
提出日時 | 2022-03-18 22:07:58 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3 ms / 2,000 ms |
コード長 | 1,439 bytes |
コンパイル時間 | 2,463 ms |
コンパイル使用メモリ | 200,140 KB |
最終ジャッジ日時 | 2025-01-28 10:24:24 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 15 |
ソースコード
#include <bits/stdc++.h>using namespace std;const long long MOD = 1000000007;vector<vector<long long>> matmul(vector<vector<long long>> A, vector<vector<long long>> B){int N = A.size();vector<vector<long long>> ans(N, vector<long long>(N, 0));for (int i = 0; i < N; i++){for (int j = 0; j < N; j++){for (int k = 0; k < N; k++){ans[i][k] += A[i][j] * B[j][k];ans[i][k] %= MOD;}}}return ans;}vector<vector<long long>> matexp(vector<vector<long long>> A, long long b){int N = A.size();vector<vector<long long>> ans(N, vector<long long>(N, 0));for (int i = 0; i < N; i++){ans[i][i] = 1;}while (b > 0){if (b % 2 == 1){ans = matmul(ans, A);}A = matmul(A, A);b /= 2;}return ans;}int main(){long long n;cin >> n;vector<vector<long long>> M(8, vector<long long>(8, 0));for (int i = 0; i < 2; i++){for (int j = 0; j < 2; j++){for (int k = 0; k < 2; k++){int p = i * 4 + j * 2 + k;if (i == 1){int p2 = j * 4 + k;M[p][p2]++;} else {if (j == 1){int p2 = k;M[p][p2]++;}if (k == 0){int p2 = j * 4 + 1;M[p][p2]++;}int p2 = j * 4 + 2 + k;M[p][p2]++;}}}}M = matexp(M, n);if (n % 2 == 0){cout << M[0][0] << endl;} else {cout << M[0][1] << endl;}}