結果

問題 No.1049 Zero (Exhaust)
ユーザー persimmon-persimmonpersimmon-persimmon
提出日時 2021-06-15 15:03:53
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 99 ms / 2,000 ms
コード長 743 bytes
コンパイル時間 1,061 ms
コンパイル使用メモリ 81,920 KB
実行使用メモリ 58,112 KB
最終ジャッジ日時 2024-06-08 01:53:14
合計ジャッジ時間 3,200 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 41 ms
51,968 KB
testcase_01 AC 41 ms
51,712 KB
testcase_02 AC 40 ms
51,456 KB
testcase_03 AC 84 ms
57,472 KB
testcase_04 AC 42 ms
51,712 KB
testcase_05 AC 72 ms
57,472 KB
testcase_06 AC 94 ms
57,728 KB
testcase_07 AC 48 ms
57,856 KB
testcase_08 AC 55 ms
57,728 KB
testcase_09 AC 80 ms
57,472 KB
testcase_10 AC 85 ms
58,112 KB
testcase_11 AC 83 ms
57,088 KB
testcase_12 AC 89 ms
57,216 KB
testcase_13 AC 53 ms
57,216 KB
testcase_14 AC 66 ms
57,728 KB
testcase_15 AC 83 ms
57,088 KB
testcase_16 AC 75 ms
57,216 KB
testcase_17 AC 77 ms
57,600 KB
testcase_18 AC 88 ms
57,472 KB
testcase_19 AC 78 ms
57,856 KB
testcase_20 AC 58 ms
57,472 KB
testcase_21 AC 79 ms
57,600 KB
testcase_22 AC 62 ms
57,088 KB
testcase_23 AC 97 ms
57,600 KB
testcase_24 AC 99 ms
58,112 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def main0(p,k):
	# 自明DP
	mod=10**9+7
	dp=[0]*p
	dp[0]=1
	for _ in range(k):
		ndp=[0]*p
		for j in range(p):
			# dp[j]からの遷移
			if dp[j]==0:continue
			for i in range(p):
				# (1)
				ndp[(j+i)%p]+=dp[j]
				ndp[(j+i)%p]%=mod
				# (2)
				ndp[(j*i)%p]+=dp[j]
				ndp[(j*i)%p]%=mod
		dp=ndp
	return dp[0]

def main1(p,k):
	mod=10**9+7
	# dp[i]:N=iとなる操作列の個数
	# 初期値dp[0]=0

	dp0=1 # dp[0]の値
	dp1=0 # sum(dp)-dp[0]
	for _ in range(k):
		ndp0=dp1+dp0 # (1)
		ndp0+=dp0*p%mod+dp1 # (2)
		ndp1=(dp0+dp1)*(p-1)%mod # (1)
		ndp1+=dp1*(p-1) # (2)
		dp0=ndp0%mod
		dp1=ndp1%mod

	return dp0

if __name__=='__main__':
	p,k=map(int,input().split())
	#ret0=main0(p,k)
	#print(ret0)
	ret1=main1(p,k)
	print(ret1)
0