結果
問題 | No.1073 無限すごろく |
ユーザー |
![]() |
提出日時 | 2020-06-05 21:52:34 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,200 bytes |
コンパイル時間 | 1,502 ms |
コンパイル使用メモリ | 173,932 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-12-17 14:23:41 |
合計ジャッジ時間 | 2,443 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 WA * 10 |
ソースコード
#include <bits/stdc++.h> using namespace std; long long MOD = 1000000007; long long p = 166666668; vector<vector<long long>> matrix_multiplication(vector<vector<long long>> &A, vector<vector<long long>> &B){ int a = A.size(); int b = A[0].size(); int c = B[0].size(); vector<vector<long long>> ans(a, vector<long long>(c, 0)); for (int i = 0; i < a; i++){ for (int j = 0; j < b; j++){ for (int k = 0; k < c; k++){ ans[i][k] = (ans[i][k] + A[i][j] * B[j][k]) % MOD; } } } return ans; } vector<vector<long long>> matrix_exponentiation(vector<vector<long long>> &mat, int p){ int N = mat.size(); vector<vector<long long>> ans(6, vector<long long>(6, 0)); for (int i = 0; i < 6; i++){ ans[i][i] = 1; } while (p > 0){ if (p % 2 == 1){ ans = matrix_multiplication(ans, mat); } mat = matrix_multiplication(mat, mat); p = p / 2; } return ans; } int main(){ long long N; cin >> N; vector<vector<long long>> mat(6, vector<long long>(6, 0)); for (int i = 0; i < 6; i++){ mat[0][i] = p; } for (int i = 1; i < 6; i++){ mat[i][i - 1] = 1; } mat = matrix_exponentiation(mat, N); cout << mat[0][0] << endl; }