fbnq矩阵快速幂

#include <iostream>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define MOD 1000000007
typedef long long ll;
using namespace std;
struct mat
{
    ll a[2][2];
};
mat mat_mul(mat x,mat y)
{
    mat res;
    memset(res.a,0,sizeof(res.a));
    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++)
        for(int k=0;k<2;k++)
        res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;
    return res;
}
int  mat_pow(int n,int f1,int f2)//矩阵快速幂求解第n个fbnq数,f1=f2=1对mod取模
{
    mat c,res;
    c.a[0][0]=c.a[0][1]=c.a[1][0]=1;
    c.a[1][1]=0;
    memset(res.a,0,sizeof(res.a));
    res.a[0][0]=f2;res.a[0][1]=f2;res.a[1][0]=f1;res.a[1][1]=f1;
    while(n)
    {
        if(n&1) res=mat_mul(res,c);
        c=mat_mul(c,c);
        n=n>>1;
    }
    return res.a[0][1];
}
int main()
{
   long long int n;


       for(int i=0;i<100;i++)
       {int f3=mat_pow(i,1,1);
       printf("%d\n",f3);}
    return 0;
}