一个简单的矩阵乘法（快速幂）

　　//求FIB的第n项
 1 #include <iostream>
#include <vector>
#include <cstdio>

using namespace std;
#define LL long long

typedef vector<int> vec;
typedef vector<vec> mat;

mat mul(mat& A, mat& B)
{
    mat C ( A.size(),vec(B[0].size()));

    for(int i=0;i<A.size();i++)
    {
        for(int k = 0;k<B.size();k++)
        {
            for(int j = 0;j<B[0].size();j++)
            {
                C[i][j] = C[i][j] + A[i][k]*B[k][j] ;
            }
        }
    }
    return C;
}

mat pow(mat A, LL n)
{
    mat B( A.size(),vec(A.size()));

    for(int i=0;i<A.size();i++)
    {
        B[i][i] =1;
    }
    while( n>0)
    {
        if( n&1) B = mul(B,A);
        A = mul(A,A);
        n>>=1;
    }
    return B;


}

int main()
{
    LL n ;
    cin>>n;
    mat A ( 2, vec(2));
    A[0][0] =1; A[0][1] = 1;
    A[1][0] = 1; A[1][1] = 0;
    A = pow( A, n);
    printf("%d
",A[1][0]);
    return 0;
}