#include using namespace std; #define _p(...) (void)printf(__VA_ARGS__) #define forr(x,arr) for(auto&& x:arr) #define _overload3(_1,_2,_3,name,...) name #define _rep2(i,n) _rep3(i,0,n) #define _rep3(i,a,b) for(int i=int(a);i=int(a);i--) #define rrep(...) _overload3(__VA_ARGS__,_rrep3,_rrep2,)(__VA_ARGS__) #define ALL(x) (x).begin(), (x).end() #define BIT(n) (1LL<<(n)) #define SZ(x) ((int)(x).size()) #define fst first #define snd second typedef vector vi;typedef vector vvi;typedef pair pii;typedef vector vpii; typedef long long ll; typedef vector Vec; typedef vector Mat; const long long mod = 1000; // C = AB // O(N^3) Mat matMul(const Mat &l, const Mat &r) { assert(l[0].size() == r.size()); // l.col == r.row int nrow = (int) l.size(); int ncol = (int) r[0].size(); int x = (int) r.size(); Mat ret(nrow, Vec(ncol)); for (int i = 0; i < nrow; i++) { for (int k = 0; k < x; k++) { for (int j = 0; j < ncol; j++) { //(ret[i][j] += (l[i][k] * r[k][j]) % mod + mod) %= mod; //ret[i][j] += l[i][k] * r[k][j]; // mod不要 (ret[i][j] += l[i][k] * r[k][j]) %= mod; } } } return ret; } // B = A^n // O(N^3 log m) Mat matPow(Mat A, long long m) { assert(A.size() == A[0].size()); // 正方行列 int size = A.size(); Mat B(size, Vec(size, 0)); for (int i = 0; i < size; i++) B[i][i] = 1; while (m > 0) { if (m & 1) B = matMul(B, A); A = matMul(A, A); m >>= 1; } return B; } void Main() { int n; scanf("%d", &n); Mat A = {{1,3},{1,1}}; Mat B = matPow(A, n); int a = B[0][0]; a = a * 2 % mod; if (n % 2 == 0) a = (a - 1 + mod) % mod; _p("%d\n", a); } int main() { cin.tie(0); ios::sync_with_stdio(false); Main(); return 0; }