結果

問題 No.1141 田グリッド
ユーザー persimmon-persimmonpersimmon-persimmon
提出日時 2021-06-17 11:31:55
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 356 ms / 2,000 ms
コード長 2,118 bytes
コンパイル時間 425 ms
コンパイル使用メモリ 82,708 KB
実行使用メモリ 108,628 KB
最終ジャッジ日時 2025-01-02 14:56:44
合計ジャッジ時間 7,550 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

def main0(h,w,a,q,rc):
  mod=10**9+7
  ans=[]
  for r,c in rc:
    r,c=r-1,c-1
    tmp=1
    for i in range(h):
      if i==r:continue
      for j in range(w):
        if j==c:continue
        tmp*=a[i][j]
        tmp%=mod
    ans.append(tmp)
  return ans

def main1(h,w,a,q,rc):
  mod=10**9+7
  mat=[]
  mat.append([1]*(w+1))
  ary0=[]
  ary0.append([0]*(w+1))
  for i in range(h):
    tmp=[1]
    tmp0=[0]
    for j in range(w):
      if a[i][j]:
        tmp.append(tmp[-1]*a[i][j]%mod)
        tmp0.append(tmp0[-1])
      else:
        tmp.append(tmp[-1])
        tmp0.append(tmp0[-1]+1)
    mat.append(tmp)
    ary0.append(tmp0)
  for j in range(w+1):
    for i in range(h):
      mat[i+1][j]*=mat[i][j]
      mat[i+1][j]%=mod
      ary0[i+1][j]+=ary0[i][j]
  ans=[]
  #for x in mat:print(x)
  for r,c in rc:
    r,c=r-1,c-1
    # 0含むかどうかチェック
    t=ary0[r][c]
    t+=ary0[-1][c]-ary0[r+1][c]
    t+=ary0[r][-1]-ary0[r][c+1]
    t+=ary0[-1][-1]+ary0[r+1][c+1]-ary0[r+1][-1]-ary0[-1][c+1]
    if t:
      ans.append(0)
      continue
    tmp0=mat[r][c] # 左上
    tmp1=mat[-1][c]*pow(mat[r+1][c],mod-2,mod)%mod # 左下
    tmp2=mat[r][-1]*pow(mat[r][c+1],mod-2,mod)%mod # 右上
    tmp3=mat[-1][-1]*mat[r+1][c+1]%mod # 右下
    tmp3=tmp3*pow(mat[r+1][-1],mod-2,mod)*pow(mat[-1][c+1],mod-2,mod)%mod # 右下
    ans.append(tmp0*tmp1*tmp2*tmp3%mod)
  return ans

if __name__=='__main__':
  h,w=map(int,input().split())
  a=[list(map(int,input().split())) for _ in range(h)]
  q=int(input())
  rc=[list(map(int,input().split())) for _ in range(q)]
  #ret0=main0(h,w,a,q,rc)
  #print(ret0)
  ret1=main1(h,w,a,q,rc)
  print(*ret1,sep="\n")

from random import randint
if __name__=='__main__1':
  for _ in range(100):
    h=randint(2,6)
    w=randint(2,6)
    a=[[randint(0,6) for _ in range(w)] for _ in range(h)]
    q=randint(2,6)
    rc=[[randint(1,h),randint(1,w)]for _ in range(q)]
    ret0=main0(h,w,a,q,rc)
    ret1=main1(h,w,a,q,rc)
    if ret0!=ret1:
      print(h,w)
      for x in a:print(*x)
      print(q)
      for x in rc:print(*x)
      print(ret0)
      print(ret1)
      break
0