結果
問題 | No.823 Many Shifts Easy |
ユーザー | pekempey |
提出日時 | 2019-04-26 21:57:37 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 47 ms / 2,000 ms |
コード長 | 1,585 bytes |
コンパイル時間 | 595 ms |
コンパイル使用メモリ | 72,444 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-11-25 03:31:51 |
合計ジャッジ時間 | 1,204 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 5 ms
6,816 KB |
testcase_01 | AC | 5 ms
6,820 KB |
testcase_02 | AC | 5 ms
6,816 KB |
testcase_03 | AC | 47 ms
6,816 KB |
testcase_04 | AC | 5 ms
6,816 KB |
testcase_05 | AC | 25 ms
6,820 KB |
testcase_06 | AC | 40 ms
6,816 KB |
testcase_07 | AC | 5 ms
6,816 KB |
testcase_08 | AC | 40 ms
6,820 KB |
testcase_09 | AC | 10 ms
6,816 KB |
ソースコード
#include <iostream>#include <algorithm>#include <vector>using namespace std;const int MOD = 1e9 + 7;struct mint {int n;mint(int n_ = 0) : n(n_) {}};mint operator+(mint a, mint b) { a.n += b.n; if (a.n >= MOD) a.n -= MOD; return a; }mint operator-(mint a, mint b) { a.n -= b.n; if (a.n < 0) a.n += MOD; return a; }mint operator*(mint a, mint b) { return (long long)a.n * b.n % MOD; }mint &operator+=(mint &a, mint b) { return a = a + b; }mint &operator-=(mint &a, mint b) { return a = a - b; }mint &operator*=(mint &a, mint b) { return a = a * b; }mint modinv(mint n) {int a = n.n;int b = MOD;int s = 1;int t = 0;while (b != 0) {int q = a / b;a -= q * b;s -= q * t;swap(a, b);swap(s, t);}if (s < 0) s += MOD;return s;}mint operator/(mint a, mint b) { return a * modinv(b); }mint F[202020] = {1, 1};mint R[202020] = {1, 1};mint I[202020] = {0, 1};mint C(int n, int r) {if (n < 0 || r < 0 || n < r) return 0;return F[n] * R[n - r] * R[r];}mint P(int n, int r) {if (n < 0 || r < 0 || n < r) return 0;return F[n] * R[n - r];}void init() {for (int i = 2; i < 202020; i++) {I[i] = I[MOD % i] * (MOD - MOD / i);F[i] = F[i - 1] * i;R[i] = R[i - 1] * I[i];}}int main() {init();int N, K;cin >> N >> K;mint ans = 0;for (int i = 1; i <= N; i++) {mint p = 0;// not i// P(N-1,K) / P(N,K)p += P(N-1,K)/P(N,K);if (i < N) {p += P(N-2,K-2) * C(K,2) / P(N,K);}ans += i*p;}ans *= P(N,K);cout << ans.n << endl;}