n=int(input())

ans=n**2*(n+1)-n*(n+1)*(2*n+1)//3
mod=10**9+7
ans*=pow(n-1,mod-2,mod)
print(ans%mod)