#include #define MOD 666013 using namespace std; int n, sol; // aranj[n][K] = n * (n - 1) * ... * (n - K + 1) int main() { cin >> n; int aranj = (2 * n) % MOD; for(int i = 2 * n - 1; i > 0; --i) { aranj = (1LL * aranj * (i % MOD)) % MOD; sol += aranj; if(sol >= MOD) sol -= MOD; } aranj = n % MOD; for(int i = n - 1; i > 0; --i) { aranj = (1LL * aranj * (i % MOD)) % MOD; sol -= aranj; if(sol < 0) sol += MOD; sol -= aranj; if(sol < 0) sol += MOD; } cout << sol << "\n"; return 0; }