結果

問題 No.980 Fibonacci Convolution Hard
コンテスト
ユーザー Ryuhei Mori
提出日時 2020-02-01 00:19:10
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 1,000 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 447 ms
コンパイル使用メモリ 43,068 KB
実行使用メモリ 7,852 KB
最終ジャッジ日時 2026-01-04 12:21:06
合計ジャッジ時間 10,702 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 13 WA * 4
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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