結果

問題 No.117 組み合わせの数
ユーザー youse_8zlyouse_8zl
提出日時 2017-12-28 23:09:42
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 2,331 bytes
コンパイル時間 1,787 ms
コンパイル使用メモリ 168,752 KB
最終ジャッジ日時 2024-06-01 05:59:12
合計ジャッジ時間 2,169 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。（リジャッジ後のコンパイルエラーは公開されます）
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:94:17: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized]
   94 |         cout << res << endl;
      |                 ^~~
main.cpp:87:26: note: 'res' was declared here
   87 |     char c; int n, k, t, res;
      |                          ^~~
/tmp/ccz4v81q.o: in function `_GLOBAL__sub_I_fact':
main.cpp:(.text.startup+0x184): relocation truncated to fit: R_X86_64_PC32 against `.bss'
collect2: error: ld returned 1 exit status

ソースコード

diff # 

#include <bits/stdc++.h>
#define INF 1e9
#define MOD (int)1e9+7
#define debug(x) cerr << #x << ": " << x << "\n";
using namespace std;
using ll = long long;
using P = pair<int, int>;

int fact[MOD];

// return n! mod p
int factorial(int n, int p = INT_MAX){
    if(fact[n] > 0) return fact[n];
    if(n == 0) return fact[n] = 1;
    int i = 0; while(fact[i] > 0) i++;
    for(int j = i; j <= n; j++){
        if(j == 0){
            fact[j] = 1;
        }else{
            fact[j] = (ll) fact[j - 1] * j % p;
        }
    }
    return fact[n];
}

//solve ax+by=gcd(a, b)
//return gcd(a, b)
int gcd(int a, int b, int* px = NULL, int* py = NULL){
    if(py == NULL){
        if(b == 0) return a;
        return gcd(b, a % b);
    }
    int d = a;
    if( b != 0){
        d = gcd(b, a % b, py, px);
        *py -= (a / b) * *px;
    }else{
        *px = 1; *py = 0;
    }
    return d;
}

// return a^(-1) mod m
int mod_inverse(int a, int m = MOD){
    int x, y;
    gcd(a, m, &x, &y);
    return (m + x % m) % m;
}

// return a mod p (when n!=a*p^e), o(log_p n)
int mod_fact(int n, int& e, int p = MOD){
    e = 0;
    if(n == 0) return 1;

    int res = mod_fact(n / p, e, p);
    e += n / p;

    if((n / p % 2 != 0)) return (ll) res * (p - factorial(n % p, p)) % p;
    return (ll) res * factorial(n % p, p) % p;
}

// return nPk mod p
int perm(int n, int k, int p = MOD){
    if(n < 0 or k < 0 or n < k) return 0;
    int e1, e2;
    int a1 = mod_fact(n, e1, p), a2 = mod_fact(n - k, e2, p);
    return (ll) a1 * mod_inverse(a2, p) % p;
}

// return nCk mod p, o(log_p n)
int comb(int n, int k, int p = MOD){
    if(n < 0 or k < 0 or n < k) return 0;
    int e1, e2, e3;
    int a1 = mod_fact(n, e1, p), a2 = mod_fact(k, e2, p), a3 = mod_fact(n - k, e3, p);
    if(e1 > e2 + e3) return 0;
    return (ll) a1 * mod_inverse((ll)a2 * a3 % p, p) % p;
}

// return nHk mod p
int hcomb(int n, int k, int p = MOD){
    if(n < 0 or k < 0) return 0;
    if(n == 0 and k == 0) return 1;
    return comb(n + k - 1, k, p);
}

int main(void){
    char c; int n, k, t, res;
    cin >> t;
    for(int i = 0; i < t; i++){
        scanf(" %c(%d, %d)", &c, &n, &k);
        if(c == 'C') res = comb(n, k);
        else if(c == 'H') res = hcomb(n, k);
        else if(c == 'P') res = perm(n, k);
        cout << res << endl;
    }

    return 0;
}
0