結果

問題 No.498 ワープクリスタル (給料日編)
ユーザー はまやんはまやんはまやんはまやん
提出日時 2017-03-24 23:07:00
言語 C++14
(gcc 11.1.0 + boost 1.77.0)
結果
AC  
実行時間 136 ms / 2,000 ms
コード長 1,455 Byte
コンパイル時間 1,375 ms
使用メモリ 5,172 KB
最終ジャッジ日時 2021-09-27 07:42:45
合計ジャッジ時間 5,974 ms
ジャッジサーバーID
(参考情報)
judge15 / judge12
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 122 ms
5,016 KB
testcase_01 AC 121 ms
5,004 KB
testcase_02 AC 122 ms
5,020 KB
testcase_03 AC 121 ms
5,068 KB
testcase_04 AC 135 ms
4,984 KB
testcase_05 AC 122 ms
5,028 KB
testcase_06 AC 135 ms
5,016 KB
testcase_07 AC 136 ms
5,172 KB
testcase_08 AC 136 ms
5,032 KB
testcase_09 AC 136 ms
5,004 KB
testcase_10 AC 122 ms
5,020 KB
testcase_11 AC 123 ms
4,936 KB
testcase_12 AC 122 ms
5,000 KB
testcase_13 AC 123 ms
4,984 KB
testcase_14 AC 122 ms
5,000 KB
testcase_15 AC 121 ms
4,928 KB
testcase_16 AC 123 ms
5,072 KB
testcase_17 AC 124 ms
4,984 KB
testcase_18 AC 135 ms
4,984 KB
testcase_19 AC 134 ms
4,928 KB
testcase_20 AC 135 ms
5,072 KB
testcase_21 AC 135 ms
5,016 KB
testcase_22 AC 136 ms
5,016 KB
testcase_23 AC 135 ms
5,016 KB
testcase_24 AC 134 ms
4,984 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=a;i<b;i++)





int mod = 1000000007;
int add(int x, int y) { return (x += y) >= mod ? x - mod : x; }
int mul(int x, int y) { return 1LL * x * y % mod; }
int modpow(int a, long long b) {
	int ret = 1;
	while (b > 0) {
		if (b & 1) ret = 1LL * ret * a % mod;
		a = 1LL * a * a % mod; b >>= 1;
	}
	return ret;
}
int modinv(int a) { return modpow(a, mod - 2); }
int fac[201010], ifac[201010];
void initfac() {
	fac[0] = ifac[0] = 1;
	rep(i, 1, 201010) fac[i] = 1LL * i * fac[i - 1] % mod;
	rep(i, 1, 201010) ifac[i] = modinv(fac[i]);
}
//-----------------------------------------------------------------------------------
int GX, GY, K;
int X[10], Y[10], N[10];
//-----------------------------------------------------------------------------------
int ans = 0;
int use[10];
void rec(int idx) {
	if (idx == K) {
		int x = 0, y = 0;
		rep(i, 0, K) x += X[i] * use[i];
		rep(i, 0, K) y += Y[i] * use[i];
		if (x != GX) return;
		if (y != GY) return;

		int sm = 0;
		rep(i, 0, K) sm += use[i];
		int com = fac[sm];
		rep(i, 0, K) if(use[i] != 0) com = mul(com, ifac[use[i]]);
		ans = add(ans, com);
		return;
	}

	rep(i, 0, N[idx] + 1) {
		use[idx] = i;
		rec(idx + 1);
	}
}
//-----------------------------------------------------------------------------------
int main() {
	initfac();

	cin >> GX >> GY >> K;
	rep(i, 0, K) cin >> X[i] >> Y[i] >> N[i];

	rec(0);

	cout << ans << endl;
}
0