斐波那契优化 快速幂+矩阵乘法

#include <iostream>
using namespace std;

//矩阵 
struct Mta
{
    int a[2][2];
};

//矩阵相乘操作，a*b 
Mta mul(Mta a,Mta b)
{
    Mta c;
    for(int i=0;i<2;i++)
    {
        for(int j=0;j<2;j++)
        {
            c.a[i][j] = 0;
            for(int k=0;k<2;k++)
            {
                c.a[i][j]+= (a.a[i][k]*b.a[k][j])%32768;
            }    
        }
    }
    return c;
}

int main()
{
    int n;
    Mta res;
    Mta base;
    //求第n个斐波那契额数 
    while(cin>>n)
    {
    
    res.a[0][0] = 1;
    res.a[0][1] = 0;
    res.a[1][0] = 0;
    res.a[1][1] = 1;
    base.a[0][0] = 1;
    base.a[0][1] = 1;
    base.a[1][0] = 1;
    base.a[1][1] = 0;
    
    
    
    //矩阵快速幂 
    while(n!=0)
    {
        if(n&1)
        {
            res = mul(res,base);
        }
        base = mul(base,base);
        n /= 2;
    }
    
    cout<<res.a[0][1]<<endl;
}
    return 0;
}