#include<bits/stdc++.h>
using namespace std;

const long long int MOD=1e9+7;
int main(){
  long long int n;
  cin>>n;
  vector<long long int> dp(n,0); //dp[i]=(rest i+1 char);
  dp[0]=n;
  long long int sum=dp[0];
  if(n!=1) dp[1]=n*n,sum+=dp[1];
  for(int i=2;i<n;i++){
    dp[i]=sum-(i-2>=0 ? dp[i-2] : 0)+((n-1)*(n-1))%MOD+n-i+1;
    dp[i]%=MOD;
    sum+=dp[i];
    sum%=MOD;
  }
  dp[n-1]=(dp[n-1]+MOD)%MOD;
  cout<<dp[n-1]<<endl;
  return 0;
}
