結果

問題 No.2193 メガの下１桁
ユーザー shobonvipshobonvip
提出日時 2023-01-13 23:02:47
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,528 bytes
コンパイル時間 183 ms
コンパイル使用メモリ 82,584 KB
実行使用メモリ 61,880 KB
最終ジャッジ日時 2024-06-06 23:55:05
合計ジャッジ時間 2,980 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
60,932 KB
testcase_01 AC 46 ms
61,276 KB
testcase_02 AC 43 ms
61,236 KB
testcase_03 AC 42 ms
60,440 KB
testcase_04 AC 42 ms
61,120 KB
testcase_05 AC 39 ms
54,604 KB
testcase_06 AC 38 ms
54,840 KB
testcase_07 AC 41 ms
59,896 KB
testcase_08 AC 39 ms
54,540 KB
testcase_09 AC 39 ms
55,036 KB
testcase_10 AC 42 ms
60,844 KB
testcase_11 AC 42 ms
61,396 KB
testcase_12 AC 38 ms
54,236 KB
testcase_13 AC 43 ms
60,420 KB
testcase_14 AC 43 ms
61,500 KB
testcase_15 AC 38 ms
55,352 KB
testcase_16 AC 43 ms
60,432 KB
testcase_17 AC 42 ms
60,572 KB
testcase_18 WA -
testcase_19 AC 38 ms
54,560 KB
testcase_20 AC 39 ms
55,100 KB
testcase_21 WA -
testcase_22 AC 42 ms
61,332 KB
testcase_23 AC 38 ms
56,144 KB
testcase_24 AC 42 ms
60,612 KB
testcase_25 AC 40 ms
55,684 KB
testcase_26 AC 44 ms
60,832 KB
testcase_27 AC 44 ms
60,684 KB
testcase_28 AC 42 ms
60,552 KB
testcase_29 AC 42 ms
60,044 KB
testcase_30 AC 41 ms
60,644 KB
testcase_31 AC 42 ms
60,748 KB
testcase_32 AC 37 ms
55,844 KB
testcase_33 AC 40 ms
59,908 KB
testcase_34 AC 41 ms
59,980 KB
testcase_35 AC 42 ms
61,880 KB
testcase_36 AC 37 ms
54,452 KB
testcase_37 AC 43 ms
61,008 KB
testcase_38 AC 40 ms
56,204 KB
testcase_39 AC 42 ms
60,500 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 = random.randrange(0, 10)
	d = random.randrange(0, 10)
	n = random.randrange(0, 1000)
	b = random.randrange(2, 12)
	"""
	m = int(input())
	d = int(input())
	n = int(input())
	b = int(input())
	"""
	m = 9
	d = 3
	n = 82
	b = 9
	"""
	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 % nb]

		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]

	def naive(m,d,n,b):
		if n == 0:
			r = '0123456789ABCDEFG'
			return r[m % b]
		if m == 0:
			n -= 1
			if d == 0:
				m = 1
			else:
				if (m+d)%b == 0:
					m = 0
				else:
					m = pow(m+d, m, b)
		m %= b
		for i in range(n):
			if (m+d)%b != 0:
				m = pow(m+d, m, b)
			else:
				m = 0
		
		r = '0123456789ABCDEFG'

		return r[m]

	#if fast(m,d,n,b) != naive(m,d,n,b):
	#	print(m,d,n,b)
	#	print(fast(m,d,n,b), naive(m,d,n,b))
	#	break
	print(fast(m,d,n,b))
	break
0