結果

問題 No.2857 Div Array
ユーザー tassei903tassei903
提出日時 2024-08-25 14:22:37
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 209 ms / 2,000 ms
コード長 1,541 bytes
コンパイル時間 296 ms
コンパイル使用メモリ 82,252 KB
実行使用メモリ 77,464 KB
最終ジャッジ日時 2024-08-25 14:22:41
合計ジャッジ時間 3,869 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

def matmul(a, b):
    n = len(a)
    m = len(b[0])
    l = len(b)
    c = [[0] * m for _ in range(n)]
    #print(a, b)
    for i in range(n):
        for j in range(m):
            for k in range(l):
                c[i][j] += a[i][k] * b[k][j]
                c[i][j] %= mod
    return c

def matmul2(v, a):
    n = len(v)
    m = len(a[0])
    c = [0] * m
    for i in range(n):
        for j in range(m):
            c[j] += v[i] * a[i][j]
            c[j] %= mod
    return c

def matpow(a, n):
    if n == 0:
        return [[1 if i == j else 0 for j in range(len(a))] for i in range(len(a))]
    elif n == 1:
        return a
    if n % 2 == 0:
        return matpow(matmul(a, a), n // 2)
    return matmul(a, matpow(a, n - 1))

mod = 998244353
n, m, k = na()

a = []
for i in range(1, m+1):
    if not a or a[-1][0] != m // i:
        a.append((m // i, 1))
    else:
        a[-1] = (m // i, a[-1][1] + 1)
#print(len(a))
A = [[0] * len(a) for _ in range(len(a))]

for i in range(len(a)):
    for j in range(len(a)):
        if abs(a[i][0] - a[j][0]) <= k:
            A[i][j] = a[j][1]

v = [a[i][1] for i in range(len(a))]


res = matmul2(v, matpow(A, n-1))
print(sum(res) % mod)
0