結果

問題 No.793 うし数列 2
ユーザー Navier_BoltzmannNavier_Boltzmann
提出日時 2024-04-21 20:47:06
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 65 ms / 2,000 ms
コード長 1,268 bytes
コンパイル時間 326 ms
コンパイル使用メモリ 81,848 KB
実行使用メモリ 56,064 KB
最終ジャッジ日時 2024-04-21 20:47:09
合計ジャッジ時間 2,474 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 48 ms
55,424 KB
testcase_01 AC 48 ms
55,680 KB
testcase_02 AC 48 ms
55,424 KB
testcase_03 AC 48 ms
55,296 KB
testcase_04 AC 49 ms
55,360 KB
testcase_05 AC 47 ms
55,296 KB
testcase_06 AC 48 ms
55,296 KB
testcase_07 AC 54 ms
55,680 KB
testcase_08 AC 48 ms
55,680 KB
testcase_09 AC 49 ms
56,064 KB
testcase_10 AC 48 ms
55,296 KB
testcase_11 AC 48 ms
55,552 KB
testcase_12 AC 48 ms
55,552 KB
testcase_13 AC 48 ms
55,680 KB
testcase_14 AC 48 ms
55,168 KB
testcase_15 AC 48 ms
55,552 KB
testcase_16 AC 48 ms
55,680 KB
testcase_17 AC 65 ms
55,300 KB
testcase_18 AC 48 ms
56,064 KB
testcase_19 AC 49 ms
55,296 KB
testcase_20 AC 49 ms
55,168 KB
testcase_21 AC 49 ms
55,808 KB
testcase_22 AC 49 ms
55,296 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import *
from itertools import *
from functools import *
from heapq import *
import sys,math
input = sys.stdin.readline


N = int(input())
mod = 10**9 + 7
def bostan_mori(C,A,N,mod):
    #A = [a1,a2,a3,,,,](初項)
    #An = -C[1]*A_(n-1) - C[2]*A_(n-2)...
    N -= 1
    n = len(A)
    P = [0]*n
    Q = C
    for i in range(n):
        ai = A[i]
        for j in range(n-i):
            cj = C[j]
            P[i+j] = (P[i+j] + ai*cj)%mod
    
    while N:
        U = [0]*(2*n)
        for i in range(n):
            pi = P[i]
            for j in range(n):
                qj = Q[j]
                if j%2:
                    qj *= -1
                U[i+j] = (U[i+j] + pi*qj)%mod
        if N%2==0:
            P = [p for i,p in enumerate(U) if i%2==0]
        else:
            P = [p for i,p in enumerate(U) if i%2==1]
        Q2 = [0]*(2*n)
        for i in range(n):
            qi = Q[i]
            for j in range(n):
                qj = Q[j]
                if j%2:
                    qj *= -1
                Q2[i+j] = (Q2[i+j] + qi*qj)%mod
        Q = [q for i,q in enumerate(Q2) if i%2==0]
        N //= 2
    return P[0]

if N==1:
    print(13)
    exit()

inv3 = pow(3,mod-2,mod)
ans = (13+inv3)*pow(10,N-1,mod)-inv3
print(ans%mod)
0