結果
問題 | No.793 うし数列 2 |
ユーザー |
![]() |
提出日時 | 2019-02-22 22:53:24 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 920 bytes |
コンパイル時間 | 1,918 ms |
コンパイル使用メモリ | 174,996 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-11-25 19:05:29 |
合計ジャッジ時間 | 2,750 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 21 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define REP(i,a,n) for(int i=(a); i<(int)(n); i++)#define rep(i,n) REP(i,0,n)#define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it)#define ALLOF(c) (c).begin(), (c).end()typedef long long ll;typedef unsigned long long ull;const ll MOD = 1000000007LL;typedef vector<ll> vec;typedef vector<vec> mat;mat mul(mat &A, mat &B){mat C(A.size(), vec(B[0].size()));rep(i,A.size()) rep(k,B.size()){rep(j,B[0].size()){C[i][j] = (C[i][j]+A[i][k]*B[k][j])%MOD;}}return C;}mat pow(mat A, ll n){mat B(A.size(), vec(A.size()));rep(i,A.size()) B[i][i] = 1;while(n>0){if(n&1) B=mul(B,A);A=mul(A,A);n>>=1;}return B;}int main(){ll N;cin >> N;mat A = {{10,3},{0,1}};mat x = {{1,0},{1,0}};mat p = pow(A, N);mat ret = mul(p,x);cout << ret[0][0] << endl;return 0;}