算法笔记-矩阵求幂

1.矩阵求幂

#include<stdio.h>
//如果把Pow写在main函数下面的话，需要加一句函数调用
void Pow(int result[20][20],int a[20][20],int n){//注意函数的写法
    int i,j,k;
    int c[20][20];//用于保存中间结果
    for(i=0;i<n;i++)
        for(j=0;j<n;j++){
            c[i][j]=0;
        }
    for(i=0;i<n;i++){  //核心代码，求矩阵与矩阵相乘
        for(j=0;j<n;j++){
            for(k=0;k<n;k++)
                c[i][j]=c[i][j]+result[i][k]*a[k][j];
        }
    }
    for(i=0;i<n;i++){
        for(j=0;j<n;j++){
            result[i][j]=c[i][j];
        }
    }
}
int main(){
    int n,k;
    int i,j;
    int a[20][20],result[20][20];
    while(scanf("%d %d",&n ,&k)!=EOF){//注意输入俩个数时的写法
        for(i=0;i<n;i++){
            for(j=0;j<n;j++){
                scanf("%d",&a[i][j]);
                result[i][j]=a[i][j];
            }
        }
        for(i=1;i<k;i++)   
            Pow(result,a,n);
        for(i=0;i<n;i++){
            for(j=0;j<n;j++){
                printf("%d",result[i][j]);
                if(j==n-1)
                    printf("\n");
                else
                    printf(" ");
            }
        }
    }
    return 0;
}

 题目分析：

（1）先确定两矩阵相乘的求法，若要求c[i][j]，则为第一个矩阵第i行的各个元素与第二个矩阵第j列的各个元素分别相乘再相加。用代码表示出来即为：

for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			for(k=0;k<n;k++){
				c[i][j]=c[i][j]+result[i][k]*a[k][j];
			}
		}
	}

（2）求k次幂，把上述求两个矩阵相乘的方法写成函数，则k次幂需要调用k-1次，因为调用1次是两矩阵相乘，调用2次是矩阵的三次幂，所以为k-1次

注意点：

（1）虽然求的是矩阵的幂，但求时用的是result[i][j]*a[i][j]，因为经过一次运算后，求二次幂时两矩阵已不相同；

（2）题目中对输出格式有规定，所以当j=n-1时要换行，这里要有一个条件判断；

（3）如果把Pow写在main函数下面的话，需要加一句函数调用。