#include #define DEBUG 0 /* */ using namespace std; class combination { public: //Nの最大値がわかれば代入する combination(long long int combinationThreshold = 1e7) { (*this).combinationThreshold = combinationThreshold + 2; } //a^p % mod long long int powmod(long long int A, long long int P, long long int M = LLONG_MAX) { long long int ans = 1; long long int mul = A; for (; P > 0; P >>= 1, mul = (mul*mul) % M) { if ((P & 1) == 1) ans = (ans*mul) % M; } return ans; } //nCk mod m long long int cb(long long int N, long long int K, long long int M = LLONG_MAX) { if (N < combinationThreshold) {// 事前計算O(N)、答えるのはO(1) if (inv.size() == 0) { inv = vector(combinationThreshold); inv[1] = 1; for (long long int i = 2; i < combinationThreshold; i++) { inv[i] = M - (M / i)*inv[M%i] % M; } fact = vector(combinationThreshold); fact[0] = 1; revFact = vector (combinationThreshold); revFact[0] = 1; for (long long int i = 1; i < combinationThreshold; i++) { fact[i] = (fact[i - 1] * i) % M; revFact[i] = (revFact[i - 1] * inv[i]) % M; } } return (((fact[N] * revFact[K]) % M)*revFact[N - K]) % M; } else { return cbOnce(N, K, M); } } //一回だけなら最速 O(K + log(M)) long long int cbOnce(long long int N, long long int K, long long int M = LLONG_MAX) { if (K > N / 2) return cb(N, N - K, M); long long int ans = 1; long long int div = 1; for (long long int i = 0; i < K; i++) { ans *= N - i; ans %= M; div *= i + 1; div %= M; } ans *= powmod(div, M - 2, M); return ans%M; } private: vector>C;//3000以下の答え vectorfact;//階乗 vectorrevFact;//階乗の逆元 vectorinv;//mod pでの逆元 long long int combinationThreshold;//Nの最大値がわかれば入れる。 }; int main() { combination cb; long long int M; cin >> M; cout << (2017 + cb.powmod(2017 * 2017, 2017, M)) % M << endl; if (DEBUG) { int ewataeoprgjpsrteog; cin >> ewataeoprgjpsrteog; } }