結果

問題 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
コンパイル時間 342 ms
コンパイル使用メモリ 82,200 KB
実行使用メモリ 61,264 KB
最終ジャッジ日時 2024-05-03 16:08:53
合計ジャッジ時間 3,284 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
54,824 KB
testcase_01 AC 38 ms
55,004 KB
testcase_02 AC 37 ms
54,976 KB
testcase_03 AC 36 ms
56,032 KB
testcase_04 AC 39 ms
54,948 KB
testcase_05 AC 39 ms
54,384 KB
testcase_06 AC 40 ms
55,008 KB
testcase_07 AC 39 ms
55,972 KB
testcase_08 AC 40 ms
54,924 KB
testcase_09 AC 39 ms
55,108 KB
testcase_10 AC 39 ms
54,380 KB
testcase_11 AC 39 ms
55,472 KB
testcase_12 AC 39 ms
55,056 KB
testcase_13 AC 42 ms
60,136 KB
testcase_14 AC 43 ms
54,488 KB
testcase_15 AC 38 ms
54,740 KB
testcase_16 AC 37 ms
55,268 KB
testcase_17 AC 40 ms
61,064 KB
testcase_18 AC 39 ms
54,944 KB
testcase_19 AC 36 ms
55,332 KB
testcase_20 AC 35 ms
54,384 KB
testcase_21 AC 35 ms
55,380 KB
testcase_22 AC 35 ms
55,564 KB
testcase_23 AC 37 ms
54,968 KB
testcase_24 AC 37 ms
55,000 KB
testcase_25 AC 40 ms
60,520 KB
testcase_26 AC 40 ms
60,736 KB
testcase_27 AC 38 ms
54,912 KB
testcase_28 AC 38 ms
54,904 KB
testcase_29 AC 38 ms
54,512 KB
testcase_30 AC 37 ms
54,304 KB
testcase_31 AC 36 ms
55,032 KB
testcase_32 AC 35 ms
54,980 KB
testcase_33 AC 37 ms
55,740 KB
testcase_34 AC 40 ms
60,688 KB
testcase_35 AC 37 ms
54,600 KB
testcase_36 AC 38 ms
55,308 KB
testcase_37 AC 39 ms
54,380 KB
testcase_38 AC 42 ms
55,552 KB
testcase_39 AC 41 ms
54,860 KB
testcase_40 AC 43 ms
61,264 KB
testcase_41 AC 39 ms
56,292 KB
testcase_42 AC 41 ms
54,880 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