結果

問題 No.64 XORフィボナッチ数列
コンテスト
ユーザー nenuon
提出日時 2017-07-23 16:08:49
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,188 bytes
コンパイル時間 850 ms
コンパイル使用メモリ 94,880 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-09 10:57:28
合計ジャッジ時間 1,496 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 11
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <cstdio>
#include <iostream>
#include <map>
#include <cmath>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <vector>
#include <stdlib.h>
#include <stdio.h>
#include <bitset>
using namespace std;
#define FOR(I,A,B) for(int I = (A); I < (B); ++I)
typedef long long ll;
typedef vector<ll> vec;
typedef vector<vec> mat;

// 無理やり行列累乗を使って

mat mul(mat &a, mat &b) {
  mat ret(a.size(), vec(b[0].size()));
  FOR(i,0,a.size()) {
    FOR(j,0,b[0].size()) {
      FOR(k,0,b.size()) {
        ret[i][j] ^= a[i][k] * b[k][j];
      }
    }
  }
  return ret;
}

mat fastpow(mat &a, ll n) {
  mat ret(a.size(), vec(a.size()));
  FOR(i,0,a.size()) FOR(j,0,a.size()) ret[i][i] = 1;
  while(n) {
    if(n & 1) ret = mul(ret, a);
    a = mul(a, a);
    n >>= 1;
  }
  return ret;
}

int main()
{
  ll f1,f2,N;
  cin>>f1>>f2>>N;
  if(N==0) {
    cout << f1 << endl;
    return 0;
  }
  mat m(2, vec(2));
  m[0][0] = 1;m[0][1] = 1;
  m[1][0] = 1;m[1][1] = 0;
  N--;
  mat mm = fastpow(m, N);
  mat f(2,vec(1));
  f[0][0] = f2;f[1][0] = f1;
  mat ans = mul(mm, f);
  cout << ans[0][0] << endl;
  return 0;
}
0