結果

問題 No.2487 Multiple of M
ユーザー Carpenters-CatCarpenters-Cat
提出日時 2023-09-29 22:30:20
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
RE  
実行時間 -
コード長 981 bytes
コンパイル時間 154 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 11,136 KB
最終ジャッジ日時 2024-07-23 06:54:33
合計ジャッジ時間 3,370 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 AC 34 ms
10,880 KB
testcase_02 RE -
testcase_03 RE -
testcase_04 AC 34 ms
10,752 KB
testcase_05 RE -
testcase_06 RE -
testcase_07 AC 32 ms
11,008 KB
testcase_08 AC 32 ms
10,880 KB
testcase_09 RE -
testcase_10 AC 31 ms
10,752 KB
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
testcase_38 RE -
testcase_39 RE -
testcase_40 RE -
testcase_41 RE -
testcase_42 RE -
testcase_43 RE -
testcase_44 RE -
testcase_45 RE -
testcase_46 RE -
testcase_47 RE -
testcase_48 RE -
testcase_49 RE -
testcase_50 RE -
testcase_51 RE -
testcase_52 RE -
testcase_53 RE -
testcase_54 RE -
testcase_55 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import *
N, M, K = map(int, input().split())
mod = 998244353
g = gcd(M, K)
if N == 1 :
	print(0)
	exit()
if (M - 1) % mod == 0 :
	ans = 0
	if N & 1 :
		ans = 1
	else :
		ans = mod - 1
	ans = (ans * (g - 1)) % mod
	print(ans)
else :
	X = pow(M - 1, N - 1, mod)
	iv = pow(M - 1, -1, mod)
	n = N - 1
	mat = [[mod - iv, (g * iv) % mod], 0, 1]
	ret = [[1, 0], [0, 1]]
	while n > 0 :
		if n & 1 :
			ret2 = [[0, 0], [0, 0]]
			for i in range(2) :
				for j in range(2) :
					for k in range(2) :
						ret2[i][j] += ret[i][k] * mat[k][j]
						ret2[i][j] %= mod
			for i in range(2) :
				for j in range(2) :
					ret[i][j] = ret2[i][j]
		n >>= 1
		ret2 = [[0, 0], [0, 0]]
		for i in range(2) :
			for j in range(2) :
				for k in range(2) :
					ret2[i][j] += mat[i][k] * mat[k][j]
					ret2[i][j] %= mod
		for i in range(2) :
			for j in range(2) :
				mat[i][j] = ret2[i][j]
	b = (g - 1) * iv
	b %= mod
	bn = b * ret[0][0] + ret[0][1]
	bn %= mod
	print((X + mod - bn) % mod)
0