#include #include #include #include using namespace std; typedef long long ll; typedef pair P; typedef pair T; long long int INF = 2e18; const ll fact_table = 1200008; double Pi = 3.1415926535897932384626; priority_queue pql; priority_queue

pqp; priority_queue

bag; //big priority queue priority_queue ,greater > pqls; priority_queue ,greater

> pqps; //small priority queue //top pop ll dx[8]={1,0,-1,0,1,1,-1,-1}; ll dy[8]={0,1,0,-1,1,-1,-1,1}; char dir[] = "DRUL"; //↓,→,↑,← #define endl "\n" #ifdef ENJAPMA #undef endl #endif #define p(x) cout< &a, bool rev = false) { ll N = a.size(), h = 0; for (ll i = 0; 1 << i < N; i++) h++; for (ll i = 0; i < N; i++) { ll j = 0; for (ll k = 0; k < h; k++) j |= (i >> k & 1) << (h - 1 - k); if (i < j) swap(a[i], a[j]); } for (ll i = 1; i < N; i *= 2) { ll base = mypow(root, (m - 1) / i / 2, m); if (rev) base = mypow(base, m - 2, m); ll w = 1; for (ll j = 0; j < i; j++) { for (ll k = 0; k < N; k += i*2) { ll s = a[j+k], t = a[j+k+i] * w % m; a[j+k+0] = (s + t) % m; a[j+k+i] = (s - t + m) % m; } w = w * base % m; } } if (rev) { ll tmp = mypow(N, m-2, m); for(ll &v : a) v = v * tmp % m; } } //(1012924417, 5), (924844033, 5) vector conv(ll _mod, ll root, const vector &a, const vector &b) { ll s = 1, t = a.size() + b.size() - 1; while(s < t) s *= 2; vector F(s), G(s); for(ll i = 0; i < (ll)a.size(); i++) F[i] = a[i]; for(ll i = 0; i < (ll)b.size(); i++) G[i] = b[i]; DFT(_mod, root, F); DFT(_mod, root, G); for(ll i = 0; i < s; i++) F[i] = (ll)F[i] * G[i] % _mod; DFT(_mod, root, F, true); return F; } }; ll fact[fact_table + 5],rfact[fact_table + 5]; void c3_init(){ fact[0] = rfact[0] = 1; for(ll i=1; i<=fact_table; i++){ fact[i] = (fact[i-1]*i) % mod; } rfact[fact_table] = mypow(fact[fact_table],mod - 2, mod); for(ll i=fact_table; i>=1; i--){ rfact[i-1] = rfact[i] * i; rfact[i-1] %= mod; } return;} ll c3(ll n,ll r){ assert(n >= r); assert(n >= 0); assert(r >= 0); return (((fact[n] * rfact[r]) % mod ) * rfact[n-r]) % mod;} bool multicase = false; ll n,m,num,a,b,c,d,e,q; ll x[100005], y[100005], z[100005]; ll k, ans; vector g; vector h; void solve(){ ll n, a1, a2; cin >> n >> a1 >> a2; if(a1 > a2)swap(a1, a2); assert(2 <= n && n <= 100000); assert(1 <= a1 && a1 <= n - 1); assert(1 <= a2 && a2 <= n - 1); ll A = a1, B = a2 - a1, C = n - a2; ll base = A * A; base %= mod; base += (B * B) % mod; base %= mod; base += (C * C) % mod; base %= mod; base *= fact[n-1]; base %= mod; ll res = 0; for(ll i=0;i> q; }else{ q = 1; } while(q--){ // pu("Case ");pu("#");pu(testcase);pu(": "); solve(); testcase++; } // solve(); return 0; }