import numpy as np
p = 10**9+7
I = np.array([[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1]],dtype=np.int64)
def pow(A,n):
    if n==0:
        return I
    if n==1:
        return A
    if n%2==0:
        return (np.matmul(pow(A,n//2),pow(A,n//2)))%p
    else:
        return np.matmul(A,np.matmul(pow(A,(n-1)//2),pow(A,(n-1)//2))%p)%p
A = np.array([[0,1,1,1],[1,0,1,1],[1,1,0,1],[1,1,1,0]],dtype=np.int64)
N = int(input())
a = np.matmul(np.array([1,0,0,0],dtype=np.int64),pow(A,N))
print(a[0])