結果

問題 No.1706 Many Bus Stops (hard)
ユーザー tassei903tassei903
提出日時 2021-10-08 23:44:16
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 52 ms / 2,000 ms
コード長 1,183 bytes
コンパイル時間 698 ms
コンパイル使用メモリ 82,340 KB
実行使用メモリ 64,956 KB
最終ジャッジ日時 2024-07-23 08:10:07
合計ジャッジ時間 3,671 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
52,952 KB
testcase_01 AC 36 ms
52,992 KB
testcase_02 AC 50 ms
64,764 KB
testcase_03 AC 50 ms
64,112 KB
testcase_04 AC 49 ms
63,404 KB
testcase_05 AC 49 ms
64,840 KB
testcase_06 AC 48 ms
64,244 KB
testcase_07 AC 51 ms
64,516 KB
testcase_08 AC 48 ms
63,928 KB
testcase_09 AC 51 ms
63,856 KB
testcase_10 AC 51 ms
64,392 KB
testcase_11 AC 49 ms
63,348 KB
testcase_12 AC 50 ms
63,852 KB
testcase_13 AC 49 ms
63,564 KB
testcase_14 AC 51 ms
64,556 KB
testcase_15 AC 52 ms
63,332 KB
testcase_16 AC 51 ms
63,484 KB
testcase_17 AC 49 ms
63,828 KB
testcase_18 AC 50 ms
63,668 KB
testcase_19 AC 49 ms
64,680 KB
testcase_20 AC 50 ms
64,580 KB
testcase_21 AC 49 ms
63,556 KB
testcase_22 AC 49 ms
64,276 KB
testcase_23 AC 48 ms
64,048 KB
testcase_24 AC 49 ms
63,728 KB
testcase_25 AC 51 ms
64,104 KB
testcase_26 AC 50 ms
64,956 KB
testcase_27 AC 50 ms
63,780 KB
testcase_28 AC 49 ms
64,280 KB
testcase_29 AC 49 ms
63,920 KB
testcase_30 AC 49 ms
63,032 KB
testcase_31 AC 50 ms
63,844 KB
testcase_32 AC 50 ms
64,832 KB
testcase_33 AC 50 ms
64,756 KB
testcase_34 AC 49 ms
63,564 KB
testcase_35 AC 49 ms
63,040 KB
testcase_36 AC 51 ms
63,716 KB
testcase_37 AC 51 ms
64,220 KB
testcase_38 AC 50 ms
63,340 KB
testcase_39 AC 52 ms
64,656 KB
testcase_40 AC 51 ms
63,780 KB
testcase_41 AC 51 ms
63,712 KB
testcase_42 AC 52 ms
63,260 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = lambda :sys.stdin.readline()[:-1]
ni = lambda :int(input())
na = lambda :list(map(int,input().split()))
sys.setrecursionlimit(10**7)
yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")
no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")
#######################################################################



#a*bの計算
def mat_mul(a,b):
    n,m,l=len(a),len(b[0]),len(b)
    c=[[0]*m for _ in range(n)]
    for i in range(n):
        for k in range(l):
            for  j in range(m):
                c[i][j]=(c[i][j]+a[i][k]*b[k][j])%mod
    return c
 
#a^nの計算
def mat_pow(a,n):
    m=len(a)
    b=[[0]*m for _ in range(m)]
    for i in range(m):
        b[i][i]=1
    while n>0:
        if n&1:
            b=mat_mul(b,a)
        a=mat_mul(a,a)
        n>>=1
    return b



c,n,m = na()
p = mod = 10**9+7
z = pow(c,p-2,p)
s = list(zip([z,0,1,0]))
b = [[z,0,0,z],[0,z,z*(c-1),z*(c-2)],[1,0,0,0],[0,1,0,0]]


if n==1:
    print((1-pow(1-z, m, mod))%mod)
else:
    p = mat_mul(mat_pow(b,n-1),s)
    #print(mat_mul(mat_pow(b,n-1),s))
    print((1-pow(1-mat_mul(mat_pow(b,n-1),s)[0][0],m,mod))%mod)
0