有没有一种 我输入N，程序就求出矩阵M的N次方的程序？？？？我自己写了一个，多一次方就多一个for循环，理论上何以求到N次，但是觉得不够简洁。。。。。。。

#include<stdio.h>
#include<conio.h>
void main(void){
int M[5][5],Mk[5][5],MKK[5][5],MKKK[5][5],i,j,p;
M[0][0]=0;M[0][1]=0;M[0][2]=0;M[0][3]=0;M[0][4]=1;
M[1][0]=1;M[1][1]=0;M[1][2]=1;M[1][3]=0;M[1][4]=0;
M[2][0]=0;M[2][1]=0;M[2][2]=0;M[2][3]=0;M[2][4]=1;
M[3][0]=1;M[3][1]=0;M[3][2]=1;M[3][3]=0;M[3][4]=0;
M[4][0]=0;M[4][1]=1;M[4][2]=0;M[4][3]=1;M[4][4]=0;
for(i=0;i<5;i++){
  for(j=0;j<5;j++){Mk[i][j]=0;
    for( p=0;p<5;p++)
        Mk[i][j]=Mk[i][j]+M[i][p]*M[p][j];
  }
}
for(i=0;i<5;i++){
  for(j=0;j<5;j++){MKK[i][j]=0;
    for( p=0;p<5;p++)
        MKK[i][j]=MKK[i][j]+Mk[i][p]*M[p][j];
  }
}
for(i=0;i<5;i++){
  for(j=0;j<5;j++){MKKK[i][j]=0;
    for( p=0;p<5;p++)
        MKKK[i][j]=MKKK[i][j]+MKK[i][p]*M[p][j];
  }
}
for(i=0;i<5;i++){
  for(j=0;j<5;j++)
     printf("%8d",MKKK[i][j]);
     printf("\n");}
getch();
}

矩阵N次幂已经忘了怎么算了，不过就从你的程序上来看，前面三个for循环的规则都是一样的，可以编写成一个函数，根据输入的N值，循环调用函数而已

谢谢，我在根据这个思路想想

简单说说：
可以先编写anew[N][N]=a[N]*b[N]的int anew[][N] f1(int a[][N],int b[][N])的函数
然后编写anew[N][N]=a[N][N] ^P次方的程序    int anew[][N] f2(int a[][N],int p)
其中f2(int a[][N],int p)  
{
  其中若p=2  return f1(a[][N],a[][N])
  否则
  return f1(a[][N],f2(a[][N],p-1))
}

亲！虽然我没弄过矩阵的乘法，但是，个人觉得啊！
你前面那一团的东西干嘛不直接用一个数组    int M[5][5]={0,0,0,0,1,0,0,0,0,1,1,0,1,0,0,0,1,0,1,0}; 这样赋值啊！
你那个看着都有点乱乱的感觉。
而且，我也也觉得下面那三重循环 机构大致一样，何不用重载函数搞呢！ 看着都清晰点！

一个个赋值，便于查找修改，等程序写的成功了，就该用定义时赋值，如果写函数的话，矩阵的幂不知道怎么表示？有的人说用循环，有的人说递归，搞不懂了。。。

经过几天思考，现在来结贴，以下是关于矩阵N次幂的求法：

#include<stdio.h>
#include<conio.h>
#define A  5
#define B  5//A,B的值可根据具体的矩阵来确定；
void main(void){
    int M[A][B]={0,0,0,0,1,1,0,1,0,0,0,0,0,0,1,1,0,1,0,0,0,1,0,1,0},MT[A][B],
       ME[A][B]={0,0,0,0,1,1,0,1,0,0,0,0,0,0,1,1,0,1,0,0,0,1,0,1,0},i,j,p,c,k;//M[A][B]为初始矩阵，MT[A][A]为中转矩阵；ME[A][B]为结果矩阵
printf("M的()次幂:");
scanf("%d",&k);
for(c=0;c<k-1;c++){
    for(i=0;i<5;i++){
      for(j=0;j<5;j++){
         MT[i][j]=0;
         for( p=0;p<5;p++)
             MT[i][j]=MT[i][j]+ME[i][p]*M[p][j];//迭代
          }
          } 
for(i=0;i<5;i++){
      for(j=0;j<5;j++){
        ME[i][j]=*(MT[i]+j);}
      
    }
}
for(i=0;i<5;i++){
      for(j=0;j<5;j++){
        printf("%4d",ME[i][j]);}
      printf("\n");}
getch();
}

[ 本帖最后由 追寻XT 于 2014-9-17 14:34 编辑 ]