結果

問題 No.344 ある無理数の累乗
ユーザー pekempeypekempey
提出日時 2016-02-13 00:08:11
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,233 bytes
コンパイル時間 3,292 ms
コンパイル使用メモリ 170,180 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-22 03:45:11
合計ジャッジ時間 3,175 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,348 KB
testcase_01 AC 1 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 1 ms
4,348 KB
testcase_04 AC 2 ms
4,348 KB
testcase_05 AC 2 ms
4,348 KB
testcase_06 AC 1 ms
4,348 KB
testcase_07 AC 2 ms
4,348 KB
testcase_08 AC 1 ms
4,348 KB
testcase_09 AC 1 ms
4,348 KB
testcase_10 AC 1 ms
4,348 KB
testcase_11 AC 1 ms
4,348 KB
testcase_12 AC 1 ms
4,348 KB
testcase_13 AC 2 ms
4,348 KB
testcase_14 AC 2 ms
4,348 KB
testcase_15 AC 2 ms
4,348 KB
testcase_16 AC 2 ms
4,348 KB
testcase_17 AC 1 ms
4,348 KB
testcase_18 AC 2 ms
4,348 KB
testcase_19 AC 2 ms
4,348 KB
testcase_20 AC 1 ms
4,348 KB
testcase_21 AC 2 ms
4,348 KB
testcase_22 AC 1 ms
4,348 KB
testcase_23 AC 1 ms
4,348 KB
testcase_24 AC 1 ms
4,348 KB
testcase_25 AC 2 ms
4,348 KB
testcase_26 AC 2 ms
4,348 KB
testcase_27 AC 1 ms
4,348 KB
testcase_28 AC 1 ms
4,348 KB
testcase_29 AC 2 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0