結果
問題 |
No.978 Fibonacci Convolution Easy
|
ユーザー |
![]() |
提出日時 | 2020-08-05 10:23:33 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 93 ms / 2,000 ms |
コード長 | 1,042 bytes |
コンパイル時間 | 163 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 75,264 KB |
最終ジャッジ日時 | 2024-09-16 00:55:57 |
合計ジャッジ時間 | 2,375 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 21 |
ソースコード
import sys def input(): return sys.stdin.readline().strip() def list2d(a, b, c): return [[c] * b for i in range(a)] def list3d(a, b, c, d): return [[[d] * c for j in range(b)] for i in range(a)] def list4d(a, b, c, d, e): return [[[[e] * d for j in range(c)] for j in range(b)] for i in range(a)] def ceil(x, y=1): return int(-(-x // y)) def INT(): return int(input()) def MAP(): return map(int, input().split()) def LIST(N=None): return list(MAP()) if N is None else [INT() for i in range(N)] def Yes(): print('Yes') def No(): print('No') def YES(): print('YES') def NO(): print('NO') sys.setrecursionlimit(10 ** 9) INF = 10 ** 19 MOD = 10 ** 9 + 7 EPS = 10 ** -10 def fermat(x, y, MOD): return x * pow(y, MOD-2, MOD) % MOD N, p = MAP() if N == 1: print(0) exit() A = [0] * N A[1] = 1 for i in range(2, N): A[i] = (A[i-1]*p + A[i-2]) % MOD allsm = pow(sum(A), 2, MOD) centersm = 0 for i in range(N): centersm += pow(A[i], 2) centersm %= MOD ans = (fermat(allsm, 2, MOD) + fermat(centersm, 2, MOD)) % MOD print(ans)