#include <iostream>
#include <vector>
#include <string>
#include <algorithm>

constexpr int MOD = 1e9 + 7;

struct modint {
	int n;
	modint(int n = 0) : n(n) {}
};

modint operator+(modint a, modint b) { return modint((a.n += b.n) >= MOD ? a.n - MOD : a.n); }
modint operator-(modint a, modint b) { return modint((a.n -= b.n) < 0 ? a.n + MOD : a.n); }
modint operator*(modint a, modint b) { return modint(1LL * a.n * b.n % MOD); }
modint &operator+=(modint &a, modint b) { return a = a + b; }
modint &operator-=(modint &a, modint b) { return a = a - b; }
modint &operator*=(modint &a, modint b) { return a = a * b; }

int main() {
	int n;
	std::cin >> n;
	std::vector<std::vector<modint>> dp(n + 10, std::vector<modint>(4));
	dp[0][0] = 1;
	for (int i = 0; i < n; i++) {
		for (int j = 0; j <= 3; j++) {
			for (int k = 1; k <= 3; k++) {
				if (j != k) {
					dp[i + k][k] += dp[i][j];
				}
			}
		}
	}

	std::cout << (dp[n][1] + dp[n][2] + dp[n][3]).n << std::endl;
}