結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 54 ms
55,424 KB
testcase_01 AC 46 ms
55,552 KB
testcase_02 AC 47 ms
55,040 KB
testcase_03 AC 47 ms
55,040 KB
testcase_04 AC 46 ms
55,040 KB
testcase_05 AC 46 ms
55,424 KB
testcase_06 AC 46 ms
55,296 KB
testcase_07 AC 46 ms
54,912 KB
testcase_08 AC 46 ms
55,168 KB
testcase_09 AC 48 ms
55,424 KB
testcase_10 AC 46 ms
55,040 KB
testcase_11 AC 46 ms
55,168 KB
testcase_12 AC 45 ms
55,808 KB
testcase_13 AC 48 ms
55,552 KB
testcase_14 AC 52 ms
55,040 KB
testcase_15 AC 48 ms
55,296 KB
testcase_16 AC 47 ms
54,912 KB
testcase_17 AC 47 ms
55,040 KB
testcase_18 AC 48 ms
55,296 KB
testcase_19 AC 47 ms
55,680 KB
testcase_20 AC 46 ms
54,912 KB
testcase_21 AC 46 ms
55,168 KB
testcase_22 AC 46 ms
54,912 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