結果

問題 No.891 隣接3項間の漸化式
ユーザー QCFiumQCFium
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <stdio.h>
#include <stdint.h>
#include <assert.h>
#include <string.h>
int ri() {
int n;
scanf("%d", &n);
return n;
}
#define MOD 1000000007
int 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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0