#pragma GCC optimize ("O3") #pragma GCC target ("avx") #include "bits/stdc++.h" // define macro "/D__MAI" using namespace std; typedef long long int ll; #define xprintf(fmt,...) fprintf(stderr,fmt,__VA_ARGS__) #define debugv(v) {printf("L%d %s > ",__LINE__,#v);for(auto e:v){cout< ",__LINE__,#m);for(int x=0;x<(w);x++){cout<<(m)[x]<<" ";}cout<\n",__LINE__,#m);for(int y=0;y<(h);y++){for(int x=0;x<(w);x++){cout<<(m)[y][x]<<" ";}cout< ostream& operator <<(ostream &o, const pair p) { o << "(" << p.first << ":" << p.second << ")"; return o; } template inline size_t argmin(iterator begin, iterator end) { return distance(begin, min_element(begin, end)); } template inline size_t argmax(iterator begin, iterator end) { return distance(begin, max_element(begin, end)); } template T& maxset(T& to, const T& val) { return to = max(to, val); } template T& minset(T& to, const T& val) { return to = min(to, val); } mt19937_64 randdev(8901016); inline ll rand_range(ll l, ll h) { return uniform_int_distribution(l, h)(randdev); } #ifdef __MAI #define getchar_unlocked getchar #define putchar_unlocked putchar #endif #ifdef __VSCC #define getchar_unlocked _getchar_nolock #define putchar_unlocked _putchar_nolock #endif namespace { #define isvisiblechar(c) (0x21<=(c)&&(c)<=0x7E) class MaiScanner { public: template void input_integer(T& var) { var = 0; T sign = 1; int cc = getchar_unlocked(); for (; cc<'0' || '9'>(int& var) { input_integer(var); return *this; } inline MaiScanner& operator>>(long long& var) { input_integer(var); return *this; } inline MaiScanner& operator>>(string& var) { int cc = getchar_unlocked(); for (; !isvisiblechar(cc); cc = getchar_unlocked()); for (; isvisiblechar(cc); cc = getchar_unlocked()) var.push_back(cc); return *this; } template void in(IT begin, IT end) { for (auto it = begin; it != end; ++it) *this >> *it; } }; } MaiScanner scanner; // power // x^p 対数mod // MEMO : 逆元...powm(n,E107-2,E107) template T powm(T x, T p, T mod = 1000000007ll) { T y = 1; x = x%mod; while (0 T combim(T n, T r, T mod = 1000000007ll) { if (n / 2 < r) return combim(n, n - r); T p = 1; for (; 0> n; { // create combination table ll p = 1; cbm[0] = 1; for (ll r=1; r<=n; ++r) p = cbm[r] = (((p*(n - r + 1)) % MD)*powm(r, MD - 2)) % MD; } // とある集合Sを選ぶ. // x in S : a[x] = x // y in ~S : a[y] = t : t in S // これを数え上げれば良さそう. // |S| ごとに数え上げられそう // あるSを考えると, // |S| ** (N - |S|) // これが,|S|の combination(N, |S|) 個ある ll ans = 0; iterate(c, 1ll, n + 1) { ans += powm(c, n - c)*cbm[c]; ans %= MD; } cout << ans << endl; return 0; }