# include <bits/stdc++.h>
using namespace std;
const int mod=1e9+7,_2=(mod+1)/2;
typedef long long ll;
int n,prime[2000010],tot;
bool nopri[20000010];
ll iv1,iv2,pw[20000010],fac[10000010],inv[10000010];
ll C(int n,int m){return fac[n]*inv[n-m]%mod*inv[m]%mod;}
ll power(ll a,int b)
{
	ll ans=1;
	for(;b;b>>=1,a=a*a%mod)
		if(b&1) ans=ans*a%mod;
	return ans;
}
ll calc(int m,int a)// \int_0^m x^{2*n-1}(x-a) dx
{
	return (m*iv2+(mod-a)*iv1)%mod*pw[m]%mod;
}
int main()
{
	int ans=0;
	cin>>n;
	fac[0]=fac[1]=inv[0]=inv[1]=1;
	for(int i=2;i<=2*n+1;i++) fac[i]=fac[i-1]*i%mod;
	inv[2*n+1]=power(fac[2*n+1],mod-2);
	for(int i=2*n+1;i>=2;i--) inv[i-1]=inv[i]*i%mod;
	nopri[0]=nopri[1]=1;pw[1]=1;
	for(int i=2;i<=4*n;i++)
	{
		if(!nopri[i]) prime[++tot]=i,pw[i]=power(i,2*n);
		for(int j=1;j<=tot && i*prime[j]<=4*n;j++)
		{
			nopri[i*prime[j]]=1;
			pw[i*prime[j]]=pw[i]*pw[prime[j]]%mod;
			if(i%prime[j]==0) break;
		}
	}
	iv1=power(2*n,mod-2);iv2=power(2*n+1,mod-2);
	ll w=inv[2*n-1]*power(_2,2*n+1)%mod;
	for(int i=0;i<=2*n;i++)
	{
		ll t1=calc(max(0,2*n-2*i),2*n-2*i),t2=calc(4*n-2*i,2*n-2*i);
		ans=(ans+(i&1?mod-C(2*n,i):C(2*n,i))*(t2-2*t1+2*mod))%mod;
	}
	cout<<(ll)ans*w%mod<<endl;
	return 0;
}