题目描述 description
给定一个N阶矩阵A,输出A的M次幂(M是非负整数)
例如:
A =
1 2
3 4
A的2次幂
7 10
15 22
输入 input
第一行是一个正整数N、M(1<=N<=30, 0<=M<=5),表示矩阵A的阶数和要求的幂数
接下来N行,每行N个绝对值不超过10的非负整数,描述矩阵A的值输出 output输出共N行,每行N个整数,表示A的M次幂所对应的矩阵。相邻的数之间用一个空格隔开
样例输入 sample input
2 2
1 2
3 4
样例输出 sample output
7 10
15 22
分析:
当矩阵A的列数与矩阵B的行数相等时才可进行矩阵乘法
矩阵相乘:矩阵A*矩阵B=矩阵C
则:
矩阵C的行数等于A的行数,矩阵C的列数等于B的列数
矩阵C的第m行第n列元素等于矩阵A的第m行元素与矩阵B的第n列元素的乘积之和
代码如下
#include <stdio.h>
#include <string.h>
#define N 100
int A[N][N],t[N][N],r[N][N];
int main()
{
int n,m,i,j,k;
scanf("%d%d",&n,&m);
for(i=0;i<n;i++) //输入矩阵A
for(j=0;j<n;j++)
scanf("%d",&A[i][j]);
for(i=0;i<n;i++) //初始化矩阵r (矩阵A的幂次方)
r[i][i]=1;
while(m--) //幂次方计算 (核心段)
{
memset(t,0,sizeof(t));
for(i=0;i<n;i++)
for(j=0;j<n;j++)
for(k=0;k<n;k++) //当一个矩阵列数等于行数时才可进行自乘,即幂次计算
t[i][j]+=r[i][k]*A[k][j]; // t[i][j]= t[i][j]+r[i][k]*A[k][j];
for(i=0;i<n;i++)
for(j=0;j<n;j++)
r[i][j]=t[i][j]; //将中间矩阵t的值赋给矩阵r
} //m次幂计算后输出矩阵r
for(i=0;i<n;i++)
{
for(j=0;j<n-1;j++)
printf("%d ",r[i][j]);
printf("%d\n",r[i][j]);
}
return 0;
}
运行结果
void *memset(void *s, int ch, size_t n);
函数解释:将s中当前位置后面的n个字节 (typedef unsigned int size_t )用 ch 替换并返回 s 。
memset:作用是在一段内存块中填充某个给定的值,它是对较大的结构体或数组进行清零操作的一种最快方法.
memset()函数原型是extern void *memset(void *buffer, int c, int count) buffer:为指针或是数组,c:是赋给buffer的值,count:是buffer的长度.