結果
問題 | No.1202 お菓子の食べ方 |
ユーザー |
![]() |
提出日時 | 2020-08-30 02:55:28 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 354 ms / 2,000 ms |
コード長 | 1,745 bytes |
コンパイル時間 | 643 ms |
コンパイル使用メモリ | 82,372 KB |
実行使用メモリ | 116,988 KB |
最終ジャッジ日時 | 2024-11-14 17:51:10 |
合計ジャッジ時間 | 13,230 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 45 |
ソースコード
import sysread = sys.stdin.buffer.readreadline = sys.stdin.buffer.readlinereadlines = sys.stdin.buffer.readlinesn,m = map(int,readline().split())s = [list(map(int,i.split())) for i in readlines()]mod = 1_000_000_007## nCkのmodを求める関数# テーブルを作る(前処理)max_n = 10**6 + 10**4fac, finv, inv = [0]*max_n, [0]*max_n, [0]*max_ndef comInit(max_n):fac[0] = fac[1] = 1finv[0] = finv[1] = 1inv[1] = 1for i in range(2,max_n):fac[i] = fac[i-1]* i% modinv[i] = mod - inv[mod%i] * (mod // i) % modfinv[i] = finv[i-1] * inv[i] % modcomInit(max_n)# 二項係数の計算def com(n,k):if(n < k):return 0if( (n<0) | (k < 0)):return 0return fac[n] * (finv[k] * finv[n-k] % mod) % modfor i in range(n):s[i].append(0)s.append([0]*(m+1))ans = 0for i in range(n+1):for j in range(m+1):if(i==n)&(j==m):continueif(s[i][j] > 0):move = i+j+s[i][j]-1tmp = com(move,i) * com(move-i,j)tmp %= modans += tmp#飛び降りif(i!=0):if(s[i-1][j] - s[i][j] > 0):tmp = com(i-1+j+s[i-1][j], s[i-1][j]-1)tmp -= com(i-1+j+s[i][j], s[i][j]-1)tmp %= modtmp *= com(i-1+j,j)tmp %= modans += tmpif(j!=0):if(s[i][j-1] - s[i][j] > 0):tmp = com(i+j-1+s[i][j-1], s[i][j-1]-1)tmp -= com(i+j-1+s[i][j], s[i][j]-1)tmp %= modtmp *= com(i+j-1,i)tmp %= modans += tmpans %= mod# print(i,j,ans)print(ans)