結果

問題 No.980 Fibonacci Convolution Hard
ユーザー Ryuhei Mori
提出日時 2020-02-01 00:19:10
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 420 ms / 2,000 ms
コード長 1,000 bytes
コンパイル時間 408 ms
コンパイル使用メモリ 29,312 KB
最終ジャッジ日時 2025-01-08 21:28:27
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 17
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:48:8: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   48 |   scanf("%d%d", &p, &q);
      |   ~~~~~^~~~~~~~~~~~~~~~
main.cpp:51:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   51 |     scanf("%d", &x);
      |     ~~~~~^~~~~~~~~~

ソースコード

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