import numpy as np MOD = 10 ** 9 + 7 M = np.array([[0, 0, 1], [1, 0, 0], [1, 1, 0]]) doubling = [M] for i in range(20): doubling.append(doubling[-1] @ doubling[-1]) doubling[-1] = np.mod(doubling[-1], MOD) N = int(input()) - 1 V = np.array([1, 0, 0]) for i in range(20)[::-1]: if N >> i & 1: V = doubling[i] @ V V = np.mod(V, MOD) print(sum(V) % MOD)