import sys
stdin=sys.stdin

ip=lambda: int(sp())
fp=lambda: float(sp())
lp=lambda:list(map(int,stdin.readline().split()))
tp=lambda:tuple(map(int,stdin.readline().split()))
sp=lambda:stdin.readline().rstrip()
yp=lambda:print('Yes')
np=lambda:print('No')

n=ip()
mod=10**9+7
dp=[[[[0 for _ in range(2)] for _ in range(2)] for _ in range(2)] for _ in range(10**6+10)]

if n>2:
    dp[3][0][0][0]=0
    dp[3][0][0][1]=1
    dp[3][0][1][0]=1
    dp[3][0][1][1]=0
    dp[3][1][0][0]=0
    dp[3][1][0][1]=0
    dp[3][1][1][1]=0
    dp[3][1][1][0]=0

    for i in range(4,10**6+10):
        dp[i][0][0][0]=0
        dp[i][0][0][1]=(dp[i-1][0][0][0]+dp[i-1][1][0][0])%mod
        dp[i][0][1][0]=(dp[i-1][0][0][1]+dp[i-1][1][0][1])%mod
        dp[i][0][1][1]=0
        dp[i][1][0][0]=(dp[i-1][0][1][0]+dp[i-1][0][1][1])%mod
        dp[i][1][0][1]=(dp[i-1][0][1][0]+dp[i-1][1][1][0])%mod
        dp[i][1][1][1]=0
        dp[i][1][1][0]=0

else:
    if n==1:
        print(1)
    else:
        print(2)
ans=0
for i in range(2):
    for j in range(2):
        for t in range(2):
            ans+=dp[n][i][j][t]

print(ans%mod)