挑战题9 矩阵乘法

9 矩阵乘法

作者: Turbo时间限制: 1S章节: 基本练习（循环）

问题描述 :

给定一个N阶矩阵A，输出A的M次幂（M是非负整数）
　　例如：
　　A =
　　1 2
　　3 4
　　A的2次幂
　　7 10
　　15 22

输入说明 :

第一行是一个正整数N、M（1<=N<=30, 0<=M<=5），表示矩阵A的阶数和要求的幂数
　　接下来N行，每行N个绝对值不超过10的非负整数，描述矩阵A的值

输出说明 :

输出共N行，每行N个整数，表示A的M次幂所对应的矩阵。相邻的数之间用一个空格隔开

输入范例 :

2 2
1 2
3 4

输出范例 :

7 10
15 22

/**
矩阵乘法

第一行与每列相乘

第二行与每列相乘
*/

#include <stdio.h>

#define  max 100
int main(){


	int N, M;//N维  M幂 
	scanf("%d%d",&N,&M);

	int number[N][N];

	int n1[N][N];

	
	int c[N][N];
	int i, j, k;
	for (i = 0; i<N; i++)
	for (j = 0; j<N; j++)
	{
		scanf("%d", &number[i][j]);
		n1[i][j] = number[i][j];
		c[i][j]=number[i][j];
	}
    
	int count;
	for(count=2;count<=M;count++){
		
		for (i = 0; i<N; i++){


		for (j = 0; j<N; j++){

			int s = 0;
			for (k = 0; k<N; k++)
				s +=  n1[k][j]*number[i][k] ;
			c[i][j] = s;
		}
							}
	for(i=0;i<N;i++)
		for(j=0;j<N;j++)
		n1[i][j]=c[i][j];//不要用指针匹配 它会直接改变相应的值
		
	}
	
	if(M==0)
		{
			for(i=0;i<N;i++)
				for(j=0;j<N;j++)
				{
					if(i==j)
						c[i][j]=1;
					else c[i][j]=0;
				}
		}

	for (i = 0; i<N; i++){
		for (j = 0; j<N; j++){
			if(j==N-1)
			printf("%d", c[i][j]);
			else printf("%d ", c[i][j]);
		}
		printf("\n");
	}

	return 0;
}

注意点：
1.线性代数任一个（n阶矩阵）零次幂都是单位矩阵
2.指针可以直接操作二维数组元素值