結果

問題 No.117 組み合わせの数
ユーザー pianonekopianoneko
提出日時 2018-09-24 11:29:26
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 3,034 bytes
コンパイル時間 1,223 ms
コンパイル使用メモリ 118,428 KB
実行使用メモリ 19,964 KB
最終ジャッジ日時 2024-09-21 18:13:36
合計ジャッジ時間 3,124 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <iterator>
#include <limits>
#include <list>
#include <locale>
#include <map>
#include <memory>
#include <new>
#include <numeric>
#include <ostream>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <typeinfo>
#include <utility>
#include <valarray>
#include <vector>
/* 
PascalCase  クラス名
camelCase   メソッド名，関数名，変数名
snake_case  ファイル名，名前空間，std のメソッド名，変数名，型名に準じる場合
SNAKE_CASE	マクロ名，static const 定数，列挙体
*/

#define REP(i, n) for (int i = 0; i < (n); i++)
#define FOR(i, init, n) for(int i = init; i < (n); i++)
#define ALL(obj) (obj).begin(), (obj).end()
#define RALL(obj) (obj).rbegin(), (obj).rend()
#define Cout(obj) cout << obj << endl
#define Size(obj) (int)(obj).size()
#define fcout cout << fixed << setprecision(10)
#define fi first
#define se second

using namespace std;
using lint = long long int;
using vvlint = vector<vector<long long int>>;
using vlint = vector<long long int>;
using vvint = vector<vector<int>>;
using vint = vector<int>;
const int MOD = 1e9 + 7;
const int iINF = 1e9;
const long long int llINF = 1e18;
const double PI = acos(-1.0);

const int dx[4] = {1, 0, -1, 0};
const int dy[4] = {0, 1, 0, -1};

const int MAX = 1000000;

vector<lint> fac(MAX + 1);
vector<lint> modFac(MAX + 1);

lint powMod(lint n, lint p)
{
    lint ans = 1;

    while (p != 0)
    {
        if (p & 1)
            ans = ans * n % MOD;

        n = n * n % MOD;

        p = p >> 1;
    }

    return ans;
}

lint comb(lint n, lint r)
{
    if(n < r) 
        return 0;

    if (n < 0 || r < 0)
        return 0;

    return (fac[n] * modFac[r]) % MOD * modFac[n - r] % MOD;
}


lint dupc(lint n, lint r)
{
    if (n == 0 && r == 0)
        return 1;

    if (n <= 0 || r < 0)
        return 0;

    return ((fac[n + r - 1] * modFac[n - 1]) % MOD * modFac[r]) % MOD;
}

lint perm(lint n, lint r)
{
    if(r < 0 || n < r)
        return 0;

    return (fac[n] * modFac[n - r]) % MOD;
}

int main()
{
    cin.tie(0);
    ios::sync_with_stdio(false);

    fac[0] = 1;
    modFac[0] = 1;

    REP(i, MAX)
    {
        fac[i + 1] = (fac[i] * (i + 1)) % MOD;
        modFac[i + 1] = modFac[i] * powMod(i + 1, MOD - 2) % MOD;
    }

    int T;
    cin >> T;

    vector<char> com(T);
    vector<int> N(T), K(T);

    REP(i, T)
    {
        scanf(" %c(%d,%d)", &com[i], &N[i], &K[i]);
    }

    REP(i, T)
    {
        if(com[i] == 'C')
        {
            cout << comb(N[i], K[i]);
        }
        else if(com[i] == 'P')
        {
            cout << perm(N[i], K[i]);
        }
        else
        {
            cout << dupc(N[i], K[i]);
        }

        cout << endl;
    }
    return 0;
}
0