結果

問題 No.117 組み合わせの数
ユーザー fumiphysfumiphys
提出日時 2019-04-29 16:04:10
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 214 ms / 5,000 ms
コード長 3,002 bytes
コンパイル時間 1,759 ms
コンパイル使用メモリ 174,100 KB
実行使用メモリ 34,432 KB
最終ジャッジ日時 2024-06-06 20:07:12
合計ジャッジ時間 2,470 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 214 ms
34,432 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main(int, const char**)':
main.cpp:127:41: warning: 'k' may be used uninitialized [-Wmaybe-uninitialized]
  127 |     if(s[0] == 'C')cout << fac.comb(n, k) << endl;
      |                                         ^
main.cpp:118:11: note: 'k' was declared here
  118 |     ll n, k;
      |           ^
main.cpp:127:41: warning: 'n' may be used uninitialized [-Wmaybe-uninitialized]
  127 |     if(s[0] == 'C')cout << fac.comb(n, k) << endl;
      |                                         ^
main.cpp:118:8: note: 'n' was declared here
  118 |     ll n, k;
      |        ^

ソースコード

diff # 

// includes
#include <bits/stdc++.h>

// macros
#define ll long long int
#define pb emplace_back
#define mk make_pair
#define pq priority_queue
#define FOR(i, a, b) for(int i=(a);i<(b);++i)
#define rep(i, n) FOR(i, 0, n)
#define rrep(i, n) for(int i=((int)(n)-1);i>=0;i--)
#define all(x) (x).begin(),(x).end()
#define sz(x) ((int)(x).size())
#define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end())
#define FI first
#define SE second
using namespace std;

//  types
typedef pair<int, int> P;
typedef pair<ll, int> Pl;
typedef pair<ll, ll> Pll;
typedef pair<double, double> Pd;
 
// constants
const int inf = 1e9;
const ll linf = 1LL << 50;
const double EPS = 1e-10;
const int mod = 1e9 + 7;

// solve
template <class T>bool chmax(T &a, const T &b){if(a < b){a = b; return 1;} return 0;}
template <class T>bool chmin(T &a, const T &b){if(a > b){a = b; return 1;} return 0;}
template <typename T> istream &operator>>(istream &is, vector<T> &vec){for(auto &v: vec)is >> v; return is;}

template <typename T>
T power(T a, T n, T mod) {
  T res = 1;
  T tmp = n;
  T curr = a;
  while(tmp){
    if(tmp % 2 == 1){
      res = (T)((ll)res * curr % mod);
    }
    curr = (T)((ll)curr * curr % mod);
    tmp >>= 1;
  }
  return res;
}

struct Factorial{
  int n;
  const int MOD = 1e9 + 7;
  vector<long long> fac;
  vector<long long> inv_;
  Factorial(int n): n(n){
    fac.resize(n + 1);
    inv_.resize(n + 1);
    calc_factorial();
    calc_inv();
  }
  void calc_factorial(){
    fac[0] = 1;
    for(int i = 1; i <= n; i++){
      fac[i] = i * fac[i-1] % MOD;
    }
  }
  void calc_inv(){
    inv_[n] = power<long long>(fac[n], MOD - 2, MOD);
    for(int i = n - 1; i >= 0; i--){
      inv_[i] = (i + 1) * inv_[i+1] % MOD;
    }
  }
  long long& operator[](size_t i){
    if(i < 0 || i > n){
      cerr << "list index out of range" << endl;
      abort();
    }
    return fac[i];
  }
  long long inv(size_t i){
    if(i < 0 || i > n){
      cerr << "list index out of range" << endl;
      abort();
    }
    return inv_[i];
  }
  long long comb(int n, int k){
    if(n < k)return 0;
    long long res = fac[n];
    res = res * inv_[n-k] % MOD;
    res = res * inv_[k] % MOD;
    return res;
  }
  long long perm(int n, int k){
    if(n < k)return 0;
    long long res = fac[n];
    res = res * inv_[n-k] % MOD;
    return res;
  }
  long long h(int n, int k){
    if(n == 0 && k == 0)return 1;
    return comb(n + k - 1, k);
  }
};

int main(int argc, char const* argv[])
{
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  int t;
  cin >> t;
  Factorial fac(2000000);
  rep(i, t){
    string s;
    cin >> s;
    string t = s.substr(2, sz(s) - 3);
    ll n, k;
    for(int j = 0; j < sz(t); j++){
      if(t[j] == ','){
        n = stol(t.substr(0, j));
        k = stol(t.substr(j+1, sz(t) - j - 1));
        break;
      }
    }

    if(s[0] == 'C')cout << fac.comb(n, k) << endl;
    else if(s[0] == 'P')cout << fac.perm(n, k) << endl;
    else cout << fac.h(n, k) << endl;
  }
	return 0;
}
0