結果

問題 No.498 ワープクリスタル (給料日編)
ユーザー TawaraTawara
提出日時 2016-03-16 22:02:36
言語 PyPy2
(7.3.15)
結果
AC  
実行時間 148 ms / 2,000 ms
コード長 1,302 bytes
コンパイル時間 1,798 ms
コンパイル使用メモリ 76,340 KB
実行使用メモリ 104,704 KB
最終ジャッジ日時 2024-04-23 18:57:52
合計ジャッジ時間 5,180 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 71 ms
75,804 KB
testcase_01 AC 71 ms
75,436 KB
testcase_02 AC 72 ms
75,536 KB
testcase_03 AC 76 ms
76,668 KB
testcase_04 AC 147 ms
104,320 KB
testcase_05 AC 71 ms
76,584 KB
testcase_06 AC 147 ms
103,716 KB
testcase_07 AC 148 ms
103,808 KB
testcase_08 AC 138 ms
104,064 KB
testcase_09 AC 135 ms
104,440 KB
testcase_10 AC 70 ms
75,676 KB
testcase_11 AC 118 ms
86,272 KB
testcase_12 AC 73 ms
75,520 KB
testcase_13 AC 114 ms
85,700 KB
testcase_14 AC 73 ms
76,604 KB
testcase_15 AC 80 ms
78,096 KB
testcase_16 AC 74 ms
76,940 KB
testcase_17 AC 76 ms
78,096 KB
testcase_18 AC 139 ms
104,032 KB
testcase_19 AC 141 ms
104,424 KB
testcase_20 AC 142 ms
104,704 KB
testcase_21 AC 139 ms
104,576 KB
testcase_22 AC 144 ms
104,064 KB
testcase_23 AC 138 ms
104,320 KB
testcase_24 AC 141 ms
104,524 KB
権限があれば一括ダウンロードができます

ソースコード

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)
    used = set()
    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
                nk = k+(mul<<shift)
                nxt = (dp[k][0]+mx,dp[k][1]+my)
                dp[nk]=nxt
                if nxt == G: used.add(nk)
            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
                nk = k+(N<<shift)
                nxt = (dp[k][0]+mx,dp[k][1]+my)
                dp[nk]=nxt
                if nxt == G: used.add(nk)
    ans = 0 if G != (0,0) else 1
    for k in used:
        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