結果

問題 No.1331 Moving Penguin
ユーザー H3PO4
提出日時 2022-03-19 09:12:21
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 761 ms / 1,500 ms
コード長 606 bytes
コンパイル時間 411 ms
コンパイル使用メモリ 81,920 KB
実行使用メモリ 88,064 KB
最終ジャッジ日時 2024-11-07 15:24:28
合計ジャッジ時間 30,356 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 49
権限があれば一括ダウンロードができます

ソースコード

diff # 

def isqrt(n):
    x = n
    y = (x + 1) // 2
    while y < x:
        x = y
        y = (x + n // x) // 2
    return x


N = int(input())
A = map(int, input().split())
sq = isqrt(N)
MOD = 10 ** 9 + 7

dp = [0] * N
dp[0] = 1

small_a = [[0] * i for i in range(sq)]

for i, a in enumerate(A):
    for b in range(1, sq):
        dp[i] += small_a[b][i % b]
    dp[i] %= MOD

    if i + 1 < N and a != 1:
        dp[i + 1] += dp[i]

    if sq <= a:
        for j in range(i + a, N, a):
            dp[j] += dp[i]
    else:
        small_a[a][i % a] += dp[i]
        small_a[a][i % a] %= MOD
print(dp[-1] % MOD)
0