求：简单又困难的 “斐波那契数列”代码

下面给出上篇博客的代码解释具体的我也在注释里面写清楚了。

至于矩阵构造嘛。。还是要看个人悟性（也有运气），显然这个我还是不行的，这个矩阵初始化我复制的。

这里我没有弄取模之类的东西。其实弄取模也很简单的。在快速幂那个地方加上去弄就行了

#include <cstdio>
#include <cstring>
const int M = 1E9 + 7;
struct Matrix {
  long long a[2][2];
  Matrix() { memset(a, 0, sizeof(a)); } //构造矩阵
  //定义矩阵乘法
  Matrix operator * (const Matrix y) {  
    Matrix ans;
    for (int i = 0; i <= 1; ++ i)
      for (int j = 0; j <= 1; ++ j)
        for (int k = 0; k <= 1; ++ k)
          ans.a[i][j] += a[i][k] * y.a[k][j];
    for (int i = 0; i <= 1; ++i)
      for (int j = 0; j <= 1; ++ j)
        ans.a[i][j] %= M;
    return ans;
  }
  //重载等于号
  void operator = (const Matrix b) {
    for (int i = 0; i <= 1; ++i)
      for (int j = 0; j <= 1; ++j)
        a[i][j] = b.a[i][j];
  }
};
//矩阵快速幂
int fast_pow(long long x) {
  Matrix ans, trs;
  ans.a[0][0] = ans.a[1][1] = 1;  //单位矩阵
  trs.a[0][0] = trs.a[1][0] = trs.a[0][1] = 1; //神奇的"斐波那契矩阵"
  //同数字的快速幂一样
  while (x) {
  //若指数为奇数，则用单位矩阵来纪录
    if (x & 1)  
      ans = ans * trs;
    trs = trs * trs;
    x >>= 1;
  }
  //因为最后指数一定会变为1，且矩阵遵守结合律
  //所以开始的单位矩阵在最后就是结果
  return ans.a[0][0];
}
int main() {
  int n;
  scanf("%d", &n);
  //结果矩阵左上角是斐波那契数列的第n+1位，所以要算第n-1位
  printf("%d\n", fast_pow(n - 1));
  return 0;
}

题目原链接——斐波那契