斐波那契数列的第N项

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>

using namespace std;
typedef long long ll;
#define MOD 1000000009

struct Mat   //矩阵
{
   ll mat[2][2];
}t;

Mat mat_x(Mat a,Mat b) //矩阵乘法
{
    Mat ta;
    memset(ta.mat,0,sizeof(ta.mat));
    for (int i=0;i<2;i++)
        for (int j=0;j<2;j++)
    {
           for (int k=0;k<2;k++)
         ta.mat[i][j]+=(a.mat[i][k]*b.mat[k][j])%MOD;
          ta.mat[i][j] = ta.mat[i][j]%MOD;
    }

         return ta;
}

Mat mat_ksm(ll w)
{
    Mat temp=t;
     if(w<0)
        return temp;
    while (w)
    {
        if (w&1)
           temp=mat_x(temp,t);
        t=mat_x(t,t);
        w=w>>1;
    }
    return temp;
}

void init()
{
    t.mat[0][0]=1;
    t.mat[0][1]=1;
    t.mat[1][0]=1;
    t.mat[1][1]=0;
}

int main()
{
    ll n;
    ios::sync_with_stdio(false);
    cin>>n;
    init();
    Mat tt=mat_ksm(n-2);
    cout << tt.mat[0][0] << endl;
    return 0;
}