結果
問題 | No.389 ロジックパズルの組み合わせ |
ユーザー | hamray |
提出日時 | 2021-06-25 14:19:00 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 79 ms / 2,000 ms |
コード長 | 5,230 bytes |
コンパイル時間 | 1,721 ms |
コンパイル使用メモリ | 173,948 KB |
実行使用メモリ | 26,892 KB |
最終ジャッジ日時 | 2024-06-25 07:41:02 |
合計ジャッジ時間 | 6,181 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 37 ms
18,432 KB |
testcase_05 | AC | 7 ms
5,376 KB |
testcase_06 | AC | 22 ms
12,032 KB |
testcase_07 | AC | 52 ms
25,088 KB |
testcase_08 | AC | 24 ms
12,800 KB |
testcase_09 | AC | 47 ms
22,144 KB |
testcase_10 | AC | 11 ms
6,400 KB |
testcase_11 | AC | 48 ms
22,272 KB |
testcase_12 | AC | 7 ms
5,376 KB |
testcase_13 | AC | 63 ms
26,772 KB |
testcase_14 | AC | 13 ms
6,656 KB |
testcase_15 | AC | 10 ms
6,528 KB |
testcase_16 | AC | 23 ms
11,520 KB |
testcase_17 | AC | 34 ms
16,384 KB |
testcase_18 | AC | 55 ms
20,212 KB |
testcase_19 | AC | 2 ms
5,376 KB |
testcase_20 | AC | 2 ms
5,376 KB |
testcase_21 | AC | 2 ms
5,376 KB |
testcase_22 | AC | 2 ms
5,376 KB |
testcase_23 | AC | 2 ms
5,376 KB |
testcase_24 | AC | 2 ms
5,376 KB |
testcase_25 | AC | 2 ms
5,376 KB |
testcase_26 | AC | 2 ms
5,376 KB |
testcase_27 | AC | 2 ms
5,376 KB |
testcase_28 | AC | 2 ms
5,376 KB |
testcase_29 | AC | 2 ms
5,376 KB |
testcase_30 | AC | 2 ms
5,376 KB |
testcase_31 | AC | 2 ms
5,376 KB |
testcase_32 | AC | 2 ms
5,376 KB |
testcase_33 | AC | 2 ms
5,376 KB |
testcase_34 | AC | 2 ms
5,376 KB |
testcase_35 | AC | 2 ms
5,376 KB |
testcase_36 | AC | 2 ms
5,376 KB |
testcase_37 | AC | 2 ms
5,376 KB |
testcase_38 | AC | 2 ms
5,376 KB |
testcase_39 | AC | 2 ms
5,376 KB |
testcase_40 | AC | 2 ms
5,376 KB |
testcase_41 | AC | 2 ms
5,376 KB |
testcase_42 | AC | 2 ms
5,376 KB |
testcase_43 | AC | 2 ms
5,376 KB |
testcase_44 | AC | 2 ms
5,376 KB |
testcase_45 | AC | 2 ms
5,376 KB |
testcase_46 | AC | 2 ms
5,376 KB |
testcase_47 | AC | 2 ms
5,376 KB |
testcase_48 | AC | 2 ms
5,376 KB |
testcase_49 | AC | 59 ms
20,908 KB |
testcase_50 | AC | 13 ms
7,808 KB |
testcase_51 | AC | 6 ms
5,376 KB |
testcase_52 | AC | 47 ms
21,632 KB |
testcase_53 | AC | 14 ms
7,936 KB |
testcase_54 | AC | 56 ms
25,344 KB |
testcase_55 | AC | 10 ms
6,656 KB |
testcase_56 | AC | 7 ms
5,376 KB |
testcase_57 | AC | 35 ms
15,128 KB |
testcase_58 | AC | 65 ms
25,272 KB |
testcase_59 | AC | 50 ms
20,088 KB |
testcase_60 | AC | 21 ms
9,344 KB |
testcase_61 | AC | 79 ms
26,492 KB |
testcase_62 | AC | 35 ms
15,360 KB |
testcase_63 | AC | 57 ms
21,292 KB |
testcase_64 | AC | 50 ms
23,040 KB |
testcase_65 | AC | 12 ms
6,784 KB |
testcase_66 | AC | 17 ms
8,064 KB |
testcase_67 | AC | 33 ms
15,488 KB |
testcase_68 | AC | 55 ms
24,832 KB |
testcase_69 | AC | 39 ms
19,328 KB |
testcase_70 | AC | 28 ms
14,720 KB |
testcase_71 | AC | 18 ms
10,112 KB |
testcase_72 | AC | 46 ms
22,656 KB |
testcase_73 | AC | 20 ms
10,880 KB |
testcase_74 | AC | 19 ms
10,496 KB |
testcase_75 | AC | 19 ms
10,496 KB |
testcase_76 | AC | 27 ms
14,464 KB |
testcase_77 | AC | 28 ms
14,464 KB |
testcase_78 | AC | 12 ms
7,296 KB |
testcase_79 | AC | 16 ms
9,344 KB |
testcase_80 | AC | 11 ms
7,040 KB |
testcase_81 | AC | 4 ms
5,376 KB |
testcase_82 | AC | 43 ms
21,248 KB |
testcase_83 | AC | 22 ms
11,776 KB |
testcase_84 | AC | 12 ms
7,552 KB |
testcase_85 | AC | 49 ms
23,936 KB |
testcase_86 | AC | 44 ms
21,632 KB |
testcase_87 | AC | 6 ms
5,376 KB |
testcase_88 | AC | 2 ms
5,376 KB |
testcase_89 | AC | 32 ms
13,788 KB |
testcase_90 | AC | 64 ms
26,892 KB |
testcase_91 | AC | 60 ms
21,716 KB |
testcase_92 | AC | 14 ms
7,936 KB |
testcase_93 | AC | 49 ms
19,472 KB |
testcase_94 | AC | 36 ms
18,176 KB |
testcase_95 | AC | 39 ms
14,812 KB |
testcase_96 | AC | 62 ms
26,308 KB |
testcase_97 | AC | 54 ms
23,068 KB |
testcase_98 | AC | 50 ms
22,912 KB |
ソースコード
#include <bits/stdc++.h> //#include <atcoder/all> //using namespace atcoder; #pragma GCC target ("avx2") #pragma GCC optimization ("O3") #pragma GCC optimization ("unroll-loops") using namespace std; typedef vector<int> VI; typedef vector<VI> VVI; typedef vector<string> VS; typedef pair<int, int> PII; typedef pair<int, int> pii; typedef pair<long long, long long> PLL; typedef pair<int, PII> TIII; typedef long long ll; typedef long double ld; typedef unsigned long long ull; #define FOR(i, s, n) for (int i = s; i < (int)n; ++i) #define REP(i, n) FOR(i, 0, n) #define rep(i, a, b) for (int i = a; i < (b); ++i) #define trav(a, x) for (auto &a : x) #define all(x) x.begin(), x.end() #define MOD 1000000007 template<class T1, class T2> inline bool chmax(T1 &a, T2 b) {if (a < b) {a = b; return true;} return false;} template<class T1, class T2> inline bool chmin(T1 &a, T2 b) {if (a > b) {a = b; return true;} return false;} const double EPS = 1e-12, PI = acos(-1); const double pi = 3.141592653589793238462643383279; //ここから編集 typedef string::const_iterator State; ll GCD(ll a, ll b){ return (b==0)?a:GCD(b, a%b); } ll LCM(ll a, ll b){ return a/GCD(a, b) * b; } template< int mod > struct ModInt { int x; ModInt() : x(0) {} ModInt(int64_t y) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {} ModInt &operator+=(const ModInt &p) { if((x += p.x) >= mod) x -= mod; return *this; } ModInt &operator-=(const ModInt &p) { if((x += mod - p.x) >= mod) x -= mod; return *this; } ModInt &operator*=(const ModInt &p) { x = (int) (1LL * x * p.x % mod); return *this; } ModInt &operator/=(const ModInt &p) { *this *= p.inverse(); return *this; } ModInt operator-() const { return ModInt(-x); } ModInt operator+(const ModInt &p) const { return ModInt(*this) += p; } ModInt operator-(const ModInt &p) const { return ModInt(*this) -= p; } ModInt operator*(const ModInt &p) const { return ModInt(*this) *= p; } ModInt operator/(const ModInt &p) const { return ModInt(*this) /= p; } bool operator==(const ModInt &p) const { return x == p.x; } bool operator!=(const ModInt &p) const { return x != p.x; } ModInt inverse() const { int a = x, b = mod, u = 1, v = 0, t; while(b > 0) { t = a / b; swap(a -= t * b, b); swap(u -= t * v, v); } return ModInt(u); } ModInt pow(int64_t n) const { ModInt ret(1), mul(x); while(n > 0) { if(n & 1) ret *= mul; mul *= mul; n >>= 1; } return ret; } friend ostream &operator<<(ostream &os, const ModInt &p) { return os << p.x; } friend istream &operator>>(istream &is, ModInt &a) { int64_t t; is >> t; a = ModInt< mod >(t); return (is); } static int get_mod() { return mod; } }; using modint = ModInt< 1000000007 >; template< typename T > struct Combination { vector< T > _fact, _rfact, _inv; Combination(int sz) : _fact(sz + 1), _rfact(sz + 1), _inv(sz + 1) { _fact[0] = _rfact[sz] = _inv[0] = 1; for(int i = 1; i <= sz; i++) _fact[i] = _fact[i - 1] * i; _rfact[sz] /= _fact[sz]; for(int i = sz - 1; i >= 0; i--) _rfact[i] = _rfact[i + 1] * (i + 1); for(int i = 1; i <= sz; i++) _inv[i] = _rfact[i] * _fact[i - 1]; } inline T fact(int k) const { return _fact[k]; } inline T rfact(int k) const { return _rfact[k]; } inline T inv(int k) const { return _inv[k]; } T P(int n, int r) const { if(r < 0 || n < r) return 0; return fact(n) * rfact(n - r); } T C(int p, int q) const { if(q < 0 || p < q) return 0; return fact(p) * rfact(q) * rfact(p - q); } T H(int n, int r) const { if(n < 0 || r < 0) return (0); return r == 0 ? 1 : C(n + r - 1, r); } }; ll modpow(ll x, ll n, ll mod) { ll res = 1; while(n) { if(n&1) res = (res * x) % mod; x = (x * x) % mod; n >>= 1; } return res; } inline long long mod(long long a, long long m) { return (a % m + m) % m; } template<typename T> struct BIT{ int N; std::vector<T> node[2]; BIT(int n){ N = n; node[0].resize(N+10); node[1].resize(N+10); } /* a: 1-indexed */ void add(int a, T x, int b){ for(int i=a; i<(int)node[b].size(); i += i & -i) node[b][i] += x; } void rangeadd(int l, int r, T x) { add(l, x, 0); add(r+1, -x, 0); add(l, x*(l-1), 1); add(r+1, -x*r, 1); } /* [1, a] */ T sum(int a, int b){ T res = 0; for(int x=a; x>0; x -= x & -x) res += node[b][x]; return res; } T prefixsum(int i) { return sum(i, 0) * i - sum(i, 1); } /* [l, r] */ T rangesum(int l, int r){ return prefixsum(r) - prefixsum(l-1); } }; int main() { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); int M; cin >> M; Combination<modint> comb(2*M+1); vector<int> a; int t; while(cin >> t) { a.push_back(t); } if(a.size() == 1 && a[0] == 0) { cout << 1 << endl; return 0; } if(accumulate(all(a), 0) + (int)a.size()-1 > M) cout << "NA" << endl; else{ int block = M - (accumulate(all(a), 0) + (int)a.size()-1); cout << comb.C(block+a.size(),block) << endl; } return 0; }