結果

問題 No.1049 Zero (Exhaust)
ユーザー persimmon-persimmonpersimmon-persimmon
提出日時 2021-06-15 15:03:53
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 104 ms / 2,000 ms
コード長 743 bytes
コンパイル時間 441 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 57,856 KB
最終ジャッジ日時 2024-12-26 16:12:55
合計ジャッジ時間 3,356 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
51,712 KB
testcase_01 AC 42 ms
51,584 KB
testcase_02 AC 41 ms
51,584 KB
testcase_03 AC 81 ms
56,960 KB
testcase_04 AC 45 ms
51,456 KB
testcase_05 AC 77 ms
57,344 KB
testcase_06 AC 104 ms
57,088 KB
testcase_07 AC 56 ms
56,960 KB
testcase_08 AC 62 ms
57,728 KB
testcase_09 AC 86 ms
57,088 KB
testcase_10 AC 89 ms
57,344 KB
testcase_11 AC 92 ms
57,600 KB
testcase_12 AC 99 ms
57,856 KB
testcase_13 AC 59 ms
57,472 KB
testcase_14 AC 68 ms
57,728 KB
testcase_15 AC 83 ms
57,216 KB
testcase_16 AC 76 ms
57,600 KB
testcase_17 AC 78 ms
57,088 KB
testcase_18 AC 88 ms
57,088 KB
testcase_19 AC 82 ms
57,728 KB
testcase_20 AC 60 ms
57,216 KB
testcase_21 AC 87 ms
57,600 KB
testcase_22 AC 69 ms
57,088 KB
testcase_23 AC 102 ms
57,344 KB
testcase_24 AC 101 ms
57,472 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