結果

問題 No.181 A↑↑N mod M
ユーザー vwxyzvwxyz
提出日時 2024-05-03 16:08:48
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 43 ms / 5,000 ms
コード長 1,510 bytes
コンパイル時間 176 ms
コンパイル使用メモリ 82,124 KB
実行使用メモリ 61,052 KB
最終ジャッジ日時 2024-11-24 18:29:35
合計ジャッジ時間 3,087 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
54,980 KB
testcase_01 AC 40 ms
54,952 KB
testcase_02 AC 40 ms
54,424 KB
testcase_03 AC 40 ms
55,420 KB
testcase_04 AC 40 ms
55,364 KB
testcase_05 AC 40 ms
54,636 KB
testcase_06 AC 39 ms
54,680 KB
testcase_07 AC 39 ms
55,044 KB
testcase_08 AC 40 ms
55,160 KB
testcase_09 AC 40 ms
54,264 KB
testcase_10 AC 40 ms
55,320 KB
testcase_11 AC 39 ms
55,136 KB
testcase_12 AC 39 ms
55,052 KB
testcase_13 AC 42 ms
59,716 KB
testcase_14 AC 40 ms
55,168 KB
testcase_15 AC 40 ms
54,484 KB
testcase_16 AC 39 ms
54,132 KB
testcase_17 AC 42 ms
60,680 KB
testcase_18 AC 40 ms
54,700 KB
testcase_19 AC 38 ms
54,256 KB
testcase_20 AC 39 ms
55,300 KB
testcase_21 AC 38 ms
55,608 KB
testcase_22 AC 39 ms
54,648 KB
testcase_23 AC 39 ms
55,232 KB
testcase_24 AC 40 ms
54,984 KB
testcase_25 AC 41 ms
60,056 KB
testcase_26 AC 43 ms
60,520 KB
testcase_27 AC 40 ms
55,284 KB
testcase_28 AC 39 ms
54,980 KB
testcase_29 AC 38 ms
54,924 KB
testcase_30 AC 40 ms
54,444 KB
testcase_31 AC 39 ms
55,748 KB
testcase_32 AC 39 ms
54,628 KB
testcase_33 AC 39 ms
55,048 KB
testcase_34 AC 42 ms
61,052 KB
testcase_35 AC 40 ms
55,188 KB
testcase_36 AC 40 ms
55,436 KB
testcase_37 AC 41 ms
55,456 KB
testcase_38 AC 40 ms
54,836 KB
testcase_39 AC 40 ms
54,580 KB
testcase_40 AC 42 ms
60,460 KB
testcase_41 AC 39 ms
55,204 KB
testcase_42 AC 39 ms
54,308 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
readline=sys.stdin.readline
from collections import defaultdict
import math

def LCM(n,m):
    if n or m:
        return abs(n)*abs(m)//math.gcd(n,m)
    return 0

def Totient(N,mod):
    retu=N
    primes=list(Factorize(mod).keys())
    l=len(primes)
    dp_popcount=[0]*(1<<l)
    dp_count=[N]*(1<<l)
    for bit in range(1,1<<l):
        i=(bit&-bit).bit_length()-1
        dp_popcount[bit]=dp_popcount[bit^1<<i]+1
        dp_count[bit]=dp_count[bit^1<<i]//primes[i]
        if dp_popcount[bit]%2==0:
            retu+=dp_count[bit]
        else:
            retu-=dp_count[bit]
    return retu

def Factorize(N):
    assert N>=1
    factors=defaultdict(int)
    for p in range(2,N):
        if p**2>N:
            break
        while N%p==0:
            factors[p]+=1
            N//=p
    if N!=1:
        factors[N]+=1
    return factors

def Permutation_Product(x,n,f):
    lst=[]
    dct={}
    idx=0
    cur=x
    while True:
        if cur in dct.keys():
            loop_begin=dct[cur]
            loop_end=idx
            break
        else:
            lst.append(cur)
            dct[cur]=idx
            cur=f(cur)
            idx+=1
    if loop_end-1>=n:
        return lst[n]
    n-=loop_begin
    loop_length=loop_end-loop_begin
    lst=lst[loop_begin:]
    return lst[n%(loop_length)]

A,N,M=map(int,readline().split())
mod=1
MM=M*M
while MM>=2:
    mod=LCM(mod,MM)
    for p in Factorize(MM):
        MM=MM*(p-1)//p
ans=Permutation_Product(1,N,lambda x:pow(A,x,mod))%M
print(ans)
0