結果

問題 No.117 組み合わせの数
ユーザー masamasa
提出日時 2016-07-11 18:18:02
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 326 ms / 5,000 ms
コード長 1,504 bytes
コンパイル時間 655 ms
コンパイル使用メモリ 73,072 KB
実行使用メモリ 34,440 KB
最終ジャッジ日時 2024-04-21 12:51:33
合計ジャッジ時間 1,548 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 326 ms
34,440 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:76:23: warning: 'ans' may be used uninitialized [-Wmaybe-uninitialized]
   76 |                 printf("%d\n", ans);
      |                 ~~~~~~^~~~~~~~~~~~~
main.cpp:63:22: note: 'ans' was declared here
   63 |         int t, n, k, ans;
      |                      ^~~

ソースコード

diff # 

#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <utility>

using namespace std;

const int MOD = 1e9 + 7;
const int N_MAX = 2e6;

vector<long long> factorial;
vector<long long> inv_factorial;

void make_factorials() {
	factorial.assign(N_MAX + 1, 1);
	for (int i = 2; i <= N_MAX; i++) {
		factorial[i] = factorial[i-1] * i % MOD;
	}

	int p2 = MOD - 2;
	inv_factorial.assign(N_MAX + 1, 1);
	for (int i = 2; i <= N_MAX; i++) {
		long long tmp = factorial[i];
		for (int j = 0; j <= 30; j++) {
			if ((p2 >> j) & 1) {
				inv_factorial[i] = inv_factorial[i] * tmp % MOD;
			}
			tmp = tmp * tmp % MOD;
		}
	}

	// for (int i = 0; i <= 10; i++) {
	// 	printf("%2d, %10lld, %10lld\n", i, factorial[i], inv_factorial[i]);
	// }
}

int calc_p(int n, int k) {
	if (n < k) {
		return 0;
	}
	return factorial[n] * inv_factorial[n - k] % MOD;
}

int calc_c(int n, int k) {
	if (n < k) {
		return 0;
	}
	return factorial[n] * inv_factorial[k] % MOD * inv_factorial[n - k] % MOD;
}

int calc_h(int n, int k) {
	if (n == 0 && k == 0) {
		return 1;
	}
	return calc_c(n + k - 1, k);
}


int main() {
	make_factorials();

	int t, n, k, ans;
	char c;

	scanf("%d\n", &t);
	for (int i = 0; i < t; i++) {
		scanf("%c(%d,%d)\n", &c, &n, &k);
		switch (c) {
		case 'P': ans = calc_p(n, k); break;
		case 'C': ans = calc_c(n, k); break;
		case 'H': ans = calc_h(n, k); break;
		default: break;
		}
		// printf("%d: %d%c%d = %d\n", i, n, c, k, ans);
		printf("%d\n", ans);
	}
	return 0;
}
0