結果
問題 |
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; } }