結果

問題 No.3038 フィボナッチ数列の周期
ユーザー chocoruskchocorusk
提出日時 2019-04-01 01:12:16
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,721 bytes
コンパイル時間 1,214 ms
コンパイル使用メモリ 82,688 KB
実行使用メモリ 77,276 KB
最終ジャッジ日時 2024-05-02 22:30:40
合計ジャッジ時間 7,024 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 53 ms
59,520 KB
testcase_01 AC 55 ms
62,976 KB
testcase_02 AC 48 ms
59,392 KB
testcase_03 AC 56 ms
62,464 KB
testcase_04 AC 49 ms
59,776 KB
testcase_05 AC 54 ms
61,696 KB
testcase_06 AC 68 ms
65,536 KB
testcase_07 AC 117 ms
76,648 KB
testcase_08 AC 124 ms
77,052 KB
testcase_09 AC 114 ms
76,416 KB
testcase_10 AC 117 ms
76,500 KB
testcase_11 AC 115 ms
76,160 KB
testcase_12 AC 124 ms
77,276 KB
testcase_13 AC 113 ms
76,800 KB
testcase_14 AC 123 ms
76,544 KB
testcase_15 AC 115 ms
76,556 KB
testcase_16 TLE -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
権限があれば一括ダウンロードができます

ソースコード

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