結果

問題 No.498 ワープクリスタル (給料日編)
ユーザー Tawara
提出日時 2016-03-14 19:26:44
言語 Python2
(2.7.18)
結果
AC  
実行時間 568 ms / 2,000 ms
コード長 932 bytes
コンパイル時間 725 ms
コンパイル使用メモリ 6,944 KB
実行使用メモリ 148,860 KB
最終ジャッジ日時 2024-10-15 19:43:32
合計ジャッジ時間 8,271 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import log
def solve():
	p = 10**9+7; F = [1]; FI = [0]*16
	for i in xrange(75): F.append((i+1)*F[i]%p)
	FI[15] = pow(F[15],p-2,p)
	for i in xrange(15,0,-1): FI[i-1] = i*FI[i]%p

	Gx,Gy,K = map(int,raw_input().split())
	G = (Gx,Gy); dp = [None]*(1<<K*4); dp[0] = (0,0)
	for i in xrange(K):
		x,y,N = map(int,raw_input().split())
		mul = 1; shift = i*4
		for j in xrange(int(log(N,2))):
			mx,my = x*mul,y*mul
			for k in xrange((1<<j+shift)-1,-1,-1):
				if dp[k] is None: continue
				dp[k+(mul<<shift)]=(dp[k][0]+mx,dp[k][1]+my)
			N -= mul; mul *= 2
		if N > 0:
			mx,my = x*N,y*N
			for k in xrange((1<<4+shift)-1,-1,-1):
				if dp[k] is None: continue
				dp[k+(N<<shift)]=(dp[k][0]+mx,dp[k][1]+my)
	ans = 0
	for k in xrange(1<<K*4):
		if dp[k] != G:continue
		total = 0; deno = 1
		for j in xrange(K):
			dup_n = k%16; k>>=4
			total += dup_n; deno=deno*FI[dup_n]%p
		ans = (ans + F[total]*deno)%p
	print ans
solve()
0