結果
問題 | No.891 隣接3項間の漸化式 |
ユーザー |
![]() |
提出日時 | 2019-09-15 15:09:42 |
言語 | C (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 1 ms / 2,000 ms |
コード長 | 795 bytes |
コンパイル時間 | 170 ms |
コンパイル使用メモリ | 31,104 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-06 22:45:40 |
合計ジャッジ時間 | 1,256 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
#include <stdio.h>#include <stdint.h>#include <assert.h>#include <string.h>int ri() {int n;scanf("%d", &n);return n;}#define MOD 1000000007int main() {int a, b, n;scanf("%d%d%d", &a, &b, &n);assert(0 <= a && a < MOD);assert(0 <= b && b < MOD);assert(0 <= n && n < MOD);int x[2][2] = {{0, 1}, {b, a}};int v[2] = {0, 1};for (; n; n >>= 1) {if (n & 1) {int t1 = ((int64_t) v[0] * x[0][0] + (int64_t) v[1] * x[0][1]) % MOD;int t2 = ((int64_t) v[0] * x[1][0] + (int64_t) v[1] * x[1][1]) % MOD;v[0] = t1, v[1] = t2;}int new_x[2][2];for (int i = 0; i < 2; i++)for (int j = 0; j < 2; j++) new_x[i][j] = ((int64_t) x[i][0] * x[0][j] + (int64_t) x[i][1] * x[1][j]) % MOD;memcpy(x, new_x, sizeof(new_x));}printf("%d\n", v[0]);return 0;}