結果

問題 No.498 ワープクリスタル (給料日編)
ユーザー kotamanegi
提出日時 2017-03-24 23:46:22
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,183 bytes
コンパイル時間 819 ms
コンパイル使用メモリ 105,152 KB
実行使用メモリ 81,988 KB
最終ジャッジ日時 2024-07-06 01:43:12
合計ジャッジ時間 3,730 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 19 WA * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define  _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <algorithm>
#include <utility>
#include <functional>
#include <cstring>
#include <queue>
#include <stack>
#include <math.h>
#include <iterator>
#include <vector>
#include <string>
#include <set>
#include <math.h>
#include <iostream> 
#include<map>
#include <iomanip>
#include <time.h>
#include <random>
#include <stdlib.h>
#include <list>
#include <typeinfo>
#include <list>
#include <set>
using namespace std;
#define LONG_INF 10000000000000000
#define MAX_MOD 1000000007
#define REP(i,n) for(long long i = 0;i < n;++i)
long long kaijo[10000] = {};
long long solve(long long n) {
	if (kaijo[n]) return kaijo[n];
	return kaijo[n] = (n*solve(n - 1))%MAX_MOD;
}
long long gcd(long long a, long long b) {
	if (b == 0) return a;
	return gcd(b, a%b);
}
long long lcm(long long a, long long b) {
	return a*b / gcd(a, b);
}
long long ex_gcd(long long a, long long b, long long &c, long long &d) {
	long long g = a;
	c = 1;
	d = 0;
	if (b != 0) g = ex_gcd(b, a%b,d,c), d -= (a / b) * c;
	return g;
}
long long gyaku(long long n) {
	long long x, y;
	if (ex_gcd(n, MAX_MOD, x, y) == 1) return (x + MAX_MOD) % MAX_MOD;
	return 0;
}
int main() {
#define int long long
	kaijo[1] = 1;
	int gx, gy, k;
	cin >> gx >> gy >> k;
	vector<vector<int>> wow[100];
	vector<int> hoge;
	wow[0].push_back(hoge);
	vector<pair<int, int>> magic;
	for (int i = 0;i < k;++i) {
		int x, y, n;
		cin >> x >> y >> n;
		magic.push_back(make_pair(x, y));
		for (int j = 0;j < wow[i].size();++j) {
			for (int q = 0;q <= n;++q) {
				vector<int> hogehoge = wow[i][j];
				hogehoge.push_back(q);
				wow[i + 1].push_back(hogehoge);
			}
		}
	}
	long long ans = 0;
	for (int i = 0;i < wow[k].size();++i) {
		int xs=0, ys=0,sum=0;
		for (int j = 0;j < wow[k][i].size();++j) {
			xs += wow[k][i][j] * magic[j].first;
			ys += wow[k][i][j] * magic[j].second;
			sum += wow[k][i][j];
		}
		if (xs == gx&&ys == gy) {
			long long nya = solve(sum);
			for (int j = 0;j < wow[k][i].size();++j) {
				if(wow[k][i][j] != 0)
				nya *= gyaku(solve(wow[k][i][j])%MAX_MOD);
				nya %= MAX_MOD;
			}
			ans += nya;
			ans %= MAX_MOD;
		}
	}
	cout << ans << endl;
	return 0;
}
0