#include <bits/stdc++.h>
using namespace std;
long long int MOD(long long int a, long long int b) {
  if (a > b)
    return a - b;
  else
    return b - a;
}
long long int power(long long int x, long long int y, long long int p) {
  long long int res = 1;
  x = x % p;
  while (y > 0) {
    if (y & 1) res = (res * x) % p;
    y = y >> 1;
    x = (x * x) % p;
  }
  return res;
}
long long int logg(long long int a) {
  long long int x = 0;
  while (a > 1) {
    x++;
    a /= 2;
  }
  return x;
}
void matmul(long long int C[2][2], long long int A[2][2]) {
  long long int B[2][2];
  for (int i = 0; i < 2; i++) {
    for (int j = 0; j < 2; j++) {
      B[i][j] = 0;
      for (int k = 0; k < 2; k++)
        B[i][j] = (B[i][j] + ((C[i][k] * A[k][j]) % 1000000007)) % 1000000007;
    }
  }
  for (int i = 0; i < 2; i++)
    for (int j = 0; j < 2; j++) C[i][j] = B[i][j];
}
void matexp(long long int A[2][2], long long int n) {
  long long int C[2][2] = {{1, 0}, {0, 1}};
  while (n > 0) {
    if (n % 2 == 1) matmul(C, A);
    matmul(A, A);
    n /= 2;
  }
  for (int i = 0; i < 2; i++)
    for (int j = 0; j < 2; j++) A[i][j] = C[i][j];
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  cout.tie(NULL);
  ;
  long long int n;
  cin >> n;
  if (n == 0) {
    cout << 1;
    exit(0);
  }
  long long int A[2][2];
  A[0][0] = A[1][1] = 3;
  A[0][1] = A[1][0] = 1;
  matexp(A, n);
  cout << A[0][0];
}
