結果

問題 No.980 Fibonacci Convolution Hard
ユーザー Ryuhei MoriRyuhei Mori
提出日時 2020-02-01 00:19:10
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 437 ms / 2,000 ms
コード長 1,000 bytes
コンパイル時間 181 ms
コンパイル使用メモリ 33,152 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-17 12:16:58
合計ジャッジ時間 10,068 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 437 ms
6,816 KB
testcase_01 AC 426 ms
6,940 KB
testcase_02 AC 422 ms
6,944 KB
testcase_03 AC 419 ms
6,940 KB
testcase_04 AC 413 ms
6,944 KB
testcase_05 AC 421 ms
6,940 KB
testcase_06 AC 418 ms
6,944 KB
testcase_07 AC 417 ms
6,940 KB
testcase_08 AC 423 ms
6,940 KB
testcase_09 AC 433 ms
6,940 KB
testcase_10 AC 436 ms
6,940 KB
testcase_11 AC 426 ms
6,940 KB
testcase_12 AC 424 ms
6,944 KB
testcase_13 AC 428 ms
6,940 KB
testcase_14 AC 418 ms
6,940 KB
testcase_15 AC 419 ms
6,944 KB
testcase_16 AC 358 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>

using u32 = unsigned;
using u64 = unsigned long long;

constexpr u32 mod = 1000000007;
u32 p;


/*
[x^{q-2}] [x/(1-p*x-x^2)]^2
*/

void mult(const u32 *P, const u32 *Q, u32 *R){
  for(int i = 0; i < 9; i++) R[i] = 0;
  for(int i = 0; i < 5; i++){
    for(int j = 0; j < 5; j++){
      R[i+j] = (R[i+j] + (u64) P[i] * Q[j]) % mod;
    }
  }
}

int qth(int x){
  u32 P[5]={0,1,0,0,0};
  u32 Q[5]={1,mod-p,mod-1,0,0};
  u32 R[9];
  mult(P, P, R);
  for(int i = 0; i < 5; i++) P[i] = R[i];
  mult(Q, Q, R);
  for(int i = 0; i < 5; i++) Q[i] = R[i];
  for(; x; x >>= 1){
    u32 V[5];
    for(int i = 0; i < 5; i++){
      if(i&1) V[i] = mod - Q[i];
      else V[i] = Q[i];
    }
    mult(P, V, R);
    for(int i = 0; i < 5; i++) P[i] = R[2*i+(x&1)];
    mult(Q, V, R);
    for(int i = 0; i < 5; i++) Q[i] = R[2*i];
  } 
  return P[0] % mod;
}


int main(){
  int q;
  scanf("%d%d", &p, &q);
  while(q--){
    int x;
    scanf("%d", &x);
    printf("%d\n", qth(x-2));
  }
  return 0;
}
0