結果
問題 | No.823 Many Shifts Easy |
ユーザー |
![]() |
提出日時 | 2019-04-28 02:39:13 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 10 ms / 2,000 ms |
コード長 | 744 bytes |
コンパイル時間 | 634 ms |
コンパイル使用メモリ | 67,008 KB |
最終ジャッジ日時 | 2025-01-07 03:16:13 |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 10 |
ソースコード
#include <iostream> using namespace std; typedef long long ll; const int MAX = 100010; ll inv[MAX],fact[MAX],fi[MAX],mod = 1000000007; void inverse(){ int i; inv[1] = 1; for(i=2;i<MAX;i++){ inv[i] = mod - (mod/i)*inv[mod%i]%mod; } fact[0] = fi[0] = 1; for(i=1;i<MAX;i++){ fact[i] = fact[i-1]*i%mod; fi[i] = (fi[i-1]*inv[i])%mod; } } ll npk(ll n,ll k){ if(n<0 || k<0 || n<k) return 0; return fact[n]*fi[n - k]%mod; } ll nck(ll n, ll k){ if(n<0 || k<0 || n<k) return 0; return fact[n]*fi[k]%mod*fi[n-k]%mod; } int main(){ ll i,n,k,x = 0; cin >> n >> k; inverse(); for(i=1;i<n;i++){ (x += i*(npk(n - 1,k) + nck(n - 2,k - 2)*fact[k]%mod*(mod + 1)/2%mod)) %= mod; } (x += n*(npk(n - 1,k))) %= mod; cout << x << endl; }