結果
| 問題 |
No.492 IOI数列
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-03-14 02:58:39 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 60 ms / 1,000 ms |
| コード長 | 737 bytes |
| コンパイル時間 | 358 ms |
| コンパイル使用メモリ | 82,376 KB |
| 実行使用メモリ | 65,664 KB |
| 最終ジャッジ日時 | 2024-10-15 15:01:12 |
| 合計ジャッジ時間 | 2,554 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 19 |
ソースコード
M1 = 1000000007
M2 = 101010101010101010101
A = [
[0, 1, 0],
[0, 100, 1],
[0, 0, 1]
]
def multiply(x, y, m):
z = [[0,0,0], [0,0,0], [0,0,0]]
for i in range(3):
for j in range(3):
for k in range(3):
z[i][j] += x[i][k] * y[k][j]
z[i][j] %= m
return z
def modpow(x, p, m):
r = [[1,0,0],[0,1,0],[0,0,1]]
base = x
while p > 0:
if p & 1 == 1:
r = multiply(r, base, m)
base = multiply(base, base, m)
p = p >> 1
return r
def main():
n = int(input())
x = modpow([[0,1,0],[0,100,1],[0,0,1]], n, M1)
y = modpow([[0,1,0],[0,100,1],[0,0,1]], n, M2)
r1 = (x[0][1] + x[0][2]) % M1
r2 = (y[0][1] + y[0][2]) % M2
print(r1)
print(r2)
if __name__ == '__main__':
main()