2.9 Fiboncaai 数列

#include <iostream>
#include <stdlib.h>
using namespace std;
class Matrix{
public:
    Matrix(int mm, int nn)
    {
        m=mm;n=nn;
        int i,j;
        data=new int *[m];
        for(i=0;i<m;i++)
        {
            data[i]=new int [n];
        }
        for(i=0;i<m;i++)
            for(j=0;j<n;j++)
                data[i][j]=0;
    }
    Matrix(const Matrix & src)
    {
        m=src.m;
        n=src.n;
        int i,j;
        data=new int *[m];
        for(i=0;i<m;i++)
        {
            data[i]=new int[n];
        }
        for(i=0;i<m;i++)
            for(j=0;j<n;j++)
                data[i][j]=src.data[i][j];
    }
    ~Matrix()
    {
        for(int i=0;i<m;i++)
            delete data[i];
        delete data;
    }
    Matrix & operator =(const Matrix &src);
    Matrix  operator *(const Matrix &s);
    void init();
    void input();
    void display();
private:
    int ** data;
    int m,n;
};
Matrix & Matrix::operator =(const Matrix &src)
{
    m=src.m;
    n=src.n;
    int i,j;
    data=new int *[m];
    for(i=0;i<m;i++)
    {
        data[i]=new int[n];
    }
    for(i=0;i<m;i++)
        for(j=0;j<n;j++)
            data[i][j]=src.data[i][j];
    return *this;
}
Matrix  Matrix::operator *(const Matrix &s)
{
    if(this->n!=s.m)
    {
        cout<<"can not multiply"<<endl;
        exit(0);
    }
    Matrix result(this->m,s.n);
    int i,j,k;
    for(k=0;k<m;k++)
        for(i=0;i<m;i++)
            for(j=0;j<s.n;j++)
                result.data[k][i]+=this->data[k][j]*s.data[j][i];
    return result;

}
void Matrix::input()
{
    for(int i=0;i<m;i++)
        for(int j=0;j<n;j++)
            cin>>data[i][j];
}
void Matrix::display()
{
    for(int i=0;i<m;i++)
    {
        for(int j=0;j<n;j++)
        {
            cout<<data[i][j]<<" ";
        }
        cout<<endl;
    }
}
void Matrix::init()
{
    if(m!=n)
    {
        cout<<"can not init"<<endl;
        exit(0);
    }
    for(int i=0;i<m;i++)
        data[i][i]=1;

}
int main()
{
    int n;
    Matrix i(2,2);
    i.init();//初始化单位矩阵
    Matrix a(2,2);
    cout<<"please input A"<<endl;
    a.input();
    cout<<"A:"<<endl;
    a.display();
    cout<<"n:"<<endl;
    cin>>n;
    for(;n;n>>=1)//快速幂算法 时间复杂度 O(log2^n)
    {
        if(n&1==1)
            i=i*a;
        a=a*a;

    }
    cout<<"I:"<<endl;
    i.display();
    return 1;
}