結果
問題 | No.344 ある無理数の累乗 |
ユーザー |
|
提出日時 | 2016-02-13 00:08:11 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,233 bytes |
コンパイル時間 | 2,326 ms |
コンパイル使用メモリ | 169,016 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-22 05:08:40 |
合計ジャッジ時間 | 2,480 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 30 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i, a) for (int i = 0; i < (a); i++) #define rep2(i, a, b) for (int i = (a); i < (b); i++) #define repr(i, a) for (int i = (a) - 1; i >= 0; i--) #define repr2(i, a, b) for (int i = (b) - 1; i >= (a); i--) template<class T1, class T2> bool chmin(T1 &a, T2 b) { return b < a && (a = b, true); } template<class T1, class T2> bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); } typedef long long ll; const ll mod = 1000; struct matrix { vector<vector<ll>> dat; int n; matrix(int n) : n(n), dat(n, vector<ll>(n)) {} static matrix I(int n) { matrix ret(n); rep(i, n) ret[i][i] = 1; return ret; } matrix pow(ll b) { matrix ret = I(n); for (matrix A = *this; b > 0; A = A * A, b /= 2) if (b & 1) ret = ret * A; return ret; } matrix operator*(matrix &b) { matrix ret(n); rep(i, n) rep(j, n) rep(k, n) (ret[i][j] += dat[i][k] * b[k][j]) %= mod; return ret; } vector<ll> &operator[](int i) { return dat[i]; } }; int main() { vector<int> a(100); matrix A(2); A.dat = { { 0, 1 }, { 2, 2 }, }; ll n; cin >> n; A = A.pow(n); ll ans = (A[0][0] + A[0][1]) * 2 % mod; if (n % 2 == 0) (ans += mod - 1) %= mod; cout << ans << endl; return 0; }