// {{{ Templates #include #define show(x) cerr << #x << " = " << x << endl using namespace std; using ll = long long; using pii = pair; using vi = vector; template ostream& operator<<(ostream& os, const vector& v) { os << "sz:" << v.size() << "\n["; for (const auto& p : v) { os << p << ","; } os << "]\n"; return os; } template ostream& operator<<(ostream& os, const pair& p) { os << "(" << p.first << "," << p.second << ")"; return os; } constexpr ll MOD = (ll)1e3; template constexpr T INF = numeric_limits::max() / 100; // }}} struct Matrix { Matrix() { for (int i = 0; i < 2; i++) { for (int j = 0; j < 2; j++) { table[i][j] = 0LL; } } } Matrix operator*(const Matrix& mat) const { Matrix ret; for (int i = 0; i < 2; i++) { for (int j = 0; j < 2; j++) { for (int k = 0; k < 2; k++) { ret.table[i][j] += table[i][k] * mat.table[k][j]; } ret.table[i][j] = ret.table[i][j] % MOD; } } return ret; } static Matrix unit() { Matrix ret; for (int i = 0; i < 2; i++) { ret.table[i][i] = 1LL; } return ret; } ll table[2][2]; }; Matrix power(const Matrix& mat, const ll n) { if (n == 0) { return Matrix::unit(); } if (n % 2 == 1) { return power(mat, n - 1) * mat; } else { const Matrix p = power(mat, n / 2); return p * p; } } int main() { cin.tie(0); ios::sync_with_stdio(false); ll n; cin >> n; if (n == 0) { cout << 1 << endl; return 0; } Matrix mat; mat.table[0][0] = 1; mat.table[0][1] = 3; mat.table[1][0] = 1; mat.table[1][1] = 1; const Matrix po = power(mat, n); const ll an = ((po.table[0][0] * 2) % MOD + MOD - ((n % 2 == 0) ? 1 : 0)) % MOD; cout << an << endl; return 0; }