結果
| 問題 |
No.823 Many Shifts Easy
|
| コンテスト | |
| ユーザー |
pockyny
|
| 提出日時 | 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;
}
pockyny