import sys
input=lambda:sys.stdin.readline().rstrip()
mod=pow(10,9)+7
N=int(input())
dp=[[0 for i in range(1001)] for j in range(1001)]
for i in range(1001):
	dp[i][3]=int(i%3==0)
for i in range(4,1001):
	for j in range(1001):
		for k in range(j//i+1):
			dp[j][i]+=dp[j-i*k][i-1]
			dp[j][i]%=mod
print(dp[N][1000])