結果

問題 No.1846 Good Binary Matrix
コンテスト
ユーザー asumo0729
提出日時 2022-02-18 22:52:04
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,775 ms / 2,000 ms
コード長 1,753 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 304 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 253,440 KB
最終ジャッジ日時 2024-06-29 09:40:49
合計ジャッジ時間 30,375 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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