結果

問題 No.1846 Good Binary Matrix
ユーザー asumo0729asumo0729
提出日時 2022-02-18 22:52:04
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,712 ms / 2,000 ms
コード長 1,753 bytes
コンパイル時間 406 ms
コンパイル使用メモリ 87,420 KB
実行使用メモリ 259,096 KB
最終ジャッジ日時 2023-09-11 19:55:06
合計ジャッジ時間 29,923 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 391 ms
258,064 KB
testcase_01 AC 384 ms
258,084 KB
testcase_02 AC 385 ms
257,496 KB
testcase_03 AC 385 ms
257,644 KB
testcase_04 AC 385 ms
258,048 KB
testcase_05 AC 390 ms
257,576 KB
testcase_06 AC 391 ms
257,724 KB
testcase_07 AC 389 ms
259,096 KB
testcase_08 AC 387 ms
257,944 KB
testcase_09 AC 386 ms
257,964 KB
testcase_10 AC 387 ms
257,920 KB
testcase_11 AC 390 ms
258,076 KB
testcase_12 AC 384 ms
257,712 KB
testcase_13 AC 388 ms
257,684 KB
testcase_14 AC 386 ms
257,580 KB
testcase_15 AC 387 ms
257,776 KB
testcase_16 AC 1,677 ms
257,896 KB
testcase_17 AC 1,712 ms
257,340 KB
testcase_18 AC 1,628 ms
257,680 KB
testcase_19 AC 1,652 ms
257,780 KB
testcase_20 AC 1,708 ms
257,684 KB
testcase_21 AC 852 ms
257,500 KB
testcase_22 AC 882 ms
257,620 KB
testcase_23 AC 387 ms
258,316 KB
testcase_24 AC 383 ms
257,748 KB
testcase_25 AC 1,328 ms
257,500 KB
testcase_26 AC 1,019 ms
257,452 KB
testcase_27 AC 428 ms
257,588 KB
testcase_28 AC 536 ms
259,096 KB
testcase_29 AC 802 ms
257,940 KB
testcase_30 AC 1,443 ms
257,492 KB
testcase_31 AC 1,325 ms
257,592 KB
testcase_32 AC 1,484 ms
257,608 KB
testcase_33 AC 1,634 ms
257,636 KB
testcase_34 AC 398 ms
258,192 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
from operator import itemgetter
from collections import defaultdict, deque
import heapq
import bisect
import math
import itertools
import copy

stdin=sys.stdin
sys.setrecursionlimit(10 ** 8)

ip=lambda: int(sp())
fp=lambda: float(sp())
lp=lambda:list(map(int,stdin.readline().split()))
sp=lambda:stdin.readline().rstrip()
Yp=lambda:print('Yes')
Np=lambda:print('No')
inf = 1 << 60
inf = float('inf')
mod = 10 ** 9 + 7
#mod = 998244353
eps = 1e-9
sortkey1 = itemgetter(0)
sortkey2 = lambda x: (x[0], x[1])

class Comb():
    def __init__(self, N, mod):
        self.N = N
        self.mod = mod
        self.fa = self.fa_fainv()[0]
        self.fainv = self.fa_fainv()[1]

    def fa_fainv(self):
        fa = [1]
        for i in range(1, self.N + 1):
            fa.append(fa[-1] * i % self.mod)

        fainv = [pow(fa[-1],self.mod - 2,self.mod)]
        for i in range(1, self.N + 1)[::-1]:
            fainv.append(fainv[-1] * i % self.mod)
        fainv = fainv[::-1]

        return fa, fainv

    def aCb(self, a, b):
        if b < 0 or a < b: 
            return 0
        return self.fa[a] * self.fainv[a - b] * self.fainv[b] % self.mod

    def aPb(self, a, b):
        if b < 0 or a < b: 
            return 0
        return self.fa[a] * self.fainv[a - b]

###############################################################

H, W = lp()
comb = Comb(10 ** 6 + 2, mod)
ans = 0
s = -1
inv_2 = pow(2, mod - 2, mod)
inv_2_w = pow(inv_2, W, mod)
HW_2 = pow(2, H * W, mod)
p = pow(2, H, mod)
for i in range(H + 1):
    s *= -1
    q = pow(p, mod - 2, mod)
    C = comb.aCb(H, i) * HW_2

    r = pow(1 - q, W, mod)
    z = s * C * r
    ans += z
    ans %= mod
    p *= inv_2
    p %= mod
    HW_2 *= inv_2_w
    HW_2 %= mod
    
print(ans)
0