結果
問題 |
No.389 ロジックパズルの組み合わせ
|
ユーザー |
![]() |
提出日時 | 2016-07-09 00:36:35 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 175 ms / 2,000 ms |
コード長 | 1,125 bytes |
コンパイル時間 | 1,549 ms |
コンパイル使用メモリ | 170,604 KB |
実行使用メモリ | 51,428 KB |
最終ジャッジ日時 | 2024-10-13 07:53:39 |
合計ジャッジ時間 | 13,427 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 99 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define rep(i,a,b) for(int i=a;i<b;i++) typedef long long ll; ll mo = 1000000007; const int NUM_FAC = 2000001; ll aCb(ll a, ll b) { static ll fact[NUM_FAC + 1], factr[NUM_FAC + 1], inv[NUM_FAC + 1]; if (fact[0] == 0) { inv[1] = fact[0] = factr[0] = 1; for (int i = 2; i <= NUM_FAC; ++i) inv[i] = inv[mo % i] * (mo - mo / i) % mo; for (int i = 1; i <= NUM_FAC; ++i) fact[i] = fact[i - 1] * i%mo, factr[i] = factr[i - 1] * inv[i] % mo; } if (b<0 || b>a) return 0; return factr[b] * fact[a] % mo*factr[a - b] % mo; } ll nHk(ll n, ll k) { return aCb(n + k - 1, k - 1); } //----------------------------------------------------------------- int M; vector<int> H; //----------------------------------------------------------------- int main() { cin >> M; int x; while (cin >> x) H.push_back(x); if (H.size() == 1 && H[0] == 0) { printf("1\n"); return 0; } int sum = 0; rep(i, 0, H.size()) sum += H[i]; int a = M - sum - (H.size() - 1); int b = H.size() + 1; if (a < 0) { cout << "NA" << endl; return 0; } ll ans = nHk(a, b); cout << ans << endl; }