結果

問題 No.2193 メガの下１桁
ユーザー shobonvipshobonvip
提出日時 2023-01-13 23:03:39
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 916 bytes
コンパイル時間 209 ms
コンパイル使用メモリ 82,408 KB
実行使用メモリ 62,120 KB
最終ジャッジ日時 2024-06-06 23:56:12
合計ジャッジ時間 2,865 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 44 ms
60,736 KB
testcase_01 AC 43 ms
61,272 KB
testcase_02 AC 41 ms
61,472 KB
testcase_03 AC 41 ms
60,972 KB
testcase_04 AC 41 ms
61,012 KB
testcase_05 AC 39 ms
54,764 KB
testcase_06 AC 38 ms
54,460 KB
testcase_07 AC 42 ms
60,388 KB
testcase_08 AC 39 ms
54,780 KB
testcase_09 AC 38 ms
54,720 KB
testcase_10 AC 42 ms
60,556 KB
testcase_11 AC 42 ms
60,864 KB
testcase_12 AC 40 ms
55,592 KB
testcase_13 AC 41 ms
60,004 KB
testcase_14 AC 43 ms
61,396 KB
testcase_15 AC 38 ms
54,656 KB
testcase_16 AC 42 ms
60,628 KB
testcase_17 AC 42 ms
59,716 KB
testcase_18 WA -
testcase_19 AC 37 ms
54,708 KB
testcase_20 AC 37 ms
54,980 KB
testcase_21 WA -
testcase_22 AC 41 ms
61,400 KB
testcase_23 AC 38 ms
54,492 KB
testcase_24 AC 41 ms
60,648 KB
testcase_25 AC 37 ms
54,620 KB
testcase_26 AC 40 ms
60,160 KB
testcase_27 AC 43 ms
60,984 KB
testcase_28 AC 42 ms
60,516 KB
testcase_29 AC 42 ms
59,816 KB
testcase_30 AC 42 ms
62,120 KB
testcase_31 AC 42 ms
60,156 KB
testcase_32 AC 40 ms
55,288 KB
testcase_33 AC 40 ms
59,644 KB
testcase_34 AC 40 ms
60,040 KB
testcase_35 AC 41 ms
60,780 KB
testcase_36 AC 38 ms
54,592 KB
testcase_37 AC 39 ms
60,400 KB
testcase_38 AC 36 ms
54,568 KB
testcase_39 AC 40 ms
60,632 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import random
from math import gcd

def pfact(m):
	pf = {}
	for i in range(2,int(m**0.5)+1):
		while m%i == 0:
			pf[i] = pf.get(i,0) + 1
			m //= i
	if m>1 : pf[m]=1
	return pf


def euler(m):
	#pfactのdefが必須
	tmp = pfact(m)
	ans = 1
	for i, j in tmp.items():
		ans = ans * (i**(j-1)) * (i-1)
	return ans


while True:
	m = int(input())
	d = int(input())
	n = int(input())
	b = int(input())
	def fast(m,d,n,b):
		nb = b * euler(b) // gcd(b, euler(b))
		v = [[0] * nb for i in range(61)]
		for i in range(nb):
			if (i+d)%nb != 0:
				v[0][i] = pow(i+d, i, nb)
			else:
				v[0][i] = 0

		for i in range(60):
			for j in range(nb):
				v[i+1][j] = v[i][v[i][j]]

		if n == 0:
			r = '0123456789ABCDEFG'
			return r[m % b]

		if m == 0 and d == 0:
			n -= 1
			m = 1

		m %= nb
		for i in range(61):
			if (n >> i) & 1:
				m = v[i][m]

		r = '0123456789ABCDEFG'
		return r[m % b]

	print(fast(m,d,n,b))
	break
0