#include using namespace std; #define rep(i,n) for (int i=0;i<(n);i++) #define rep2(i,a,b) for (int i=(a);i<(b);i++) #define rrep(i,n) for (int i=(n)-1;i>=0;i--) #define rrep2(i,a,b) for (int i=(b)-1;i>=(a);i--) #define all(a) (a).begin(),(a).end() typedef long long ll; typedef pair P; typedef vector vi; typedef vector

vp; typedef vector vll; typedef vector Vec; typedef vector Mat; const ll mod = 1000; Mat mulMatMat(Mat A, Mat B) { int n = A.size(); Mat C(n, Vec(n, 0)); rep(i, n) { rep(j, n) { rep(k, n) { C[i][j] += A[i][k] * B[k][j]; C[i][j] %= mod; } } } return C; } Vec mulMatVec(Mat A, Vec vec) { int n = A.size(); Vec res(n, 0); // reverse(all(vec)); // 添え字の大きい要素を上に持ってくるため rep(i, n) { rep(j, n) { res[i] += A[i][j] * vec[j]; res[i] %= mod; } } return res; } Mat powMat(Mat A, ll p) { int n = A.size(); Mat R(n, Vec(n, 0)); rep(i, n) R[i][i] = 1; // 単位行列 for (; p >= 1; p /= 2) { if (p & 1) { R = mulMatMat(A, R); } A = mulMatMat(A, A); } return R; } signed main() { std::ios::sync_with_stdio(false); std::cin.tie(0); ll n; cin >> n; if (n == 0) { cout << 1 << endl; return 0; } Mat A(2, Vec(2, 1)); A[0][1] = 3; A = powMat(A, n); auto vec = mulMatVec(A, {1, 0}); cout << (2 * vec[0] - (n % 2 == 0)) % mod << endl; }