結果

問題 No.891 隣接3項間の漸化式
コンテスト
ユーザー 37zigen
提出日時 2019-09-21 03:39:36
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 860 bytes
コンパイル時間 646 ms
コンパイル使用メモリ 74,880 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-15 11:47:10
合計ジャッジ時間 1,868 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <string>
#include <vector>
#include <queue>
#include <cmath>

const long long MOD=(long long)1e9+7;

int main() {
  long long u[2][2], v[2][2], m[2][2];
  for (int i=0;i<2;++i) for (int j=0;j<2;++j) u[i][j]=v[i][j]=m[i][j]=0;
  long long a, b, n;
  std::cin>>a>>b>>n;
  m[0][0]=a;
  m[0][1]=b;
  m[1][0]=1;
  u[0][0]=u[1][1]=1;
  for (;n>0;n>>=1) {
    if (n%2==1) {
      for (int i=0;i<2;++i) for (int j=0;j<2;++j) {
	  v[i][j]=0;
	  for (int k=0;k<2;++k) v[i][j]=(v[i][j]+u[i][k]*m[k][j]%MOD)%MOD;
	}
      for (int i=0;i<2;++i) for (int j=0;j<2;++j) u[i][j]=v[i][j];
    }
    for (int i=0;i<2;++i) for (int j=0;j<2;++j) {
	v[i][j]=0;
	for (int k=0;k<2;++k) v[i][j]=(v[i][j]+m[i][k]*m[k][j])%MOD;
      }
    for (int i=0;i<2;++i) for (int j=0;j<2;++j) m[i][j]=v[i][j];

    
  }
  std::cout<<u[1][0]<<std::endl;
  return 0;
}
0