結果
問題 |
No.344 ある無理数の累乗
|
ユーザー |
![]() |
提出日時 | 2019-06-08 21:10:07 |
言語 | C (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 1 ms / 2,000 ms |
コード長 | 831 bytes |
コンパイル時間 | 414 ms |
コンパイル使用メモリ | 30,336 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-10-05 23:53:31 |
合計ジャッジ時間 | 1,250 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 30 |
ソースコード
// yukicoder: 344 ある無理数の累乗 // 2019.6.8 bal4u #include <stdio.h> #define MOD 1000 void multiply(int F[2][2], int M[2][2]) { int x = (F[0][0] * M[0][0] + F[0][1] * M[1][0]) % MOD; int y = (F[0][0] * M[0][1] + F[0][1] * M[1][1]) % MOD; int z = (F[1][0] * M[0][0] + F[1][1] * M[1][0]) % MOD; int w = (F[1][0] * M[0][1] + F[1][1] * M[1][1]) % MOD; F[0][0] = x, F[0][1] = y, F[1][0] = z, F[1][1] = w; } void power(int F[2][2], int n) { int M[2][2] = { {1,3},{1,1} }; if (n == 0 || n == 1) return; power(F, n >> 1); multiply(F, F); if (n & 1) multiply(F, M); } int calc(int n) { int F[2][2] = { {1,3},{1,1} }; if (n == 0) return 0; power(F, n - 1); return F[0][0]; } int main() { int n, ans; scanf("%d", &n); n++; ans = (calc(n) << 1) % MOD; if (n & 1) ans--; printf("%d\n", ans); return 0; }