結果

問題 No.8038 フィボナッチ数列の周期
コンテスト
ユーザー chocorusk
提出日時 2019-04-01 01:12:16
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,721 bytes
コンパイル時間 494 ms
コンパイル使用メモリ 82,080 KB
実行使用メモリ 166,672 KB
最終ジャッジ日時 2024-11-23 17:35:41
合計ジャッジ時間 38,708 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 16 TLE * 9
権限があれば一括ダウンロードができます

ソースコード

diff # 

import math
n=int(input())
f=[[1, 1], [1, 0]]
MOD=10**9+7
def mul(a, b):
    ret=[[0, 0], [0, 0]]
    for i in range(2):
        for j in range(2):
            for k in range(2):
                ret[i][j]+=a[i][k]*b[k][j]
    return ret
def powmod(a, k, mod):
    ret=[[1, 0], [0, 1]]
    ap=a
    while k>0:
        if k&1:
            ret=mul(ret, ap)
            for i in range(2):
                for j in range(2):
                    ret[i][j]%=mod
        ap=mul(ap, ap)
        for i in range(2):
            for j in range(2):
                ap[i][j]%=mod
        k>>=1
    return ret
ans=1
def lcm(a, b):
    return a//math.gcd(a, b)*b
for i in range(n):
    p, k=map(int, input().split())
    mod=p**k
    e=(p**(2*k-2))*(p*p-1)
    x=p-1
    fac=[]
    for j in range(2, x+1):
        if j*j>x:
            break
        if x%j==0:
            t=0
            while x%j==0:
                x//=j
                t+=1
            fac.append([j, t])
    if x>1:
        fac.append([x, 1])
    x=p+1
    for j in range(2, x+1):
        if j*j>x:
            break
        if x%j==0:
            t=0
            while x%j==0:
                x//=j
                t+=1
            if j==2:
                fac[0][1]+=t
            else:
                fac.append([j, t])
    if x>1:
        fac.append([x, 1])
    for j in range(2*k-2):
        e//=p
        fp=powmod(f, e, mod)
        if fp[0][0]==1 and fp[1][0]==0:
            continue
        e*=p
        break
    for q in fac:
        for j in range(q[1]):
            e//=q[0]
            fp=powmod(f, e, mod)
            if fp[0][0]==1 and fp[1][0]==0:
                continue
            e*=q[0]
            break
    ans=lcm(ans, e)
print(ans%MOD)
0