#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; typedef pair p_ii; const int INF = 1e9; const double PI = acos(-1.0); const ll MOD = 1e9 + 7; const ll limit = 1e12; ll dp[51]; int main() { ll N; cin>>N; dp[0]=1; for (int i = 0; i < 50; i++) { dp[i+1]=dp[i]*(i+1)%limit; } if(N<50){ cout<