求教高手啊，C语言改matlab。帮帮忙，感激不尽

chenhande

#include "stdio.h"
" y* J5 G; i8 n+ {. v#include "string.h"
#define N 50
int jisuan(int a[N][N][N],int h,int z,int i,int j);
/ z  D* P! D) P$ w5 d( q: Wvoid micifang(int a[N][N][N],int h,int t);
# K- P, ]# ?/ i7 A: Cvoid buer(int a[N][N][N],int h);
void warshall(int a[N][N][N],int h);
9 h2 L5 w1 r1 p: l4 n, C- W0 Fmain()
- X' ^0 C3 N- O{
int i,h,j,t;
int a[N][N][N];
printf("请输入邻接矩阵A的行数(或列数)：\n");
scanf("%d",&h);
% J7 K" A# Z! }3 l7 t7 Wprintf("\n请输入邻接矩阵A：\n\n");
for(i=0;i<h;i++)
+ \1 ~& R% S* n" t  R: h* q{
printf("请输入邻接矩阵A的第%d行.\n",i+1);
for(j=0;j<h;j++)
6 m% X) P4 r' W+ ^8 K8 oscanf("%d",&a[0][i][j]);
}
printf("\n邻接矩阵A为：\n");
: U0 T; q& C* `& ?" y: N) Hfor(i=0;i<h;i++)
  |  ?' ~6 Q8 j+ o; C4 [{
) c3 Y9 `2 L2 W( y: H0 E+ C* Dfor(j=0;j<h;j++)
printf(" %d ",a[0][i][j]);
printf("\n");
}
! U8 t4 Q6 I0 h! @5 [! T8 jprintf("\n\n\n算法一：\n");
printf("\n矩阵A的1次方为:\n");
& q1 M$ ], W6 c, y1 O, T/ rfor(i=0;i<h;i++)
+ k1 D* d; o! ~4 _( O% J{
. r3 S7 e: O% x$ ~4 [: Bfor(j=0;j<h;j++)
9 X6 @2 h; p9 gprintf(" %d ",a[0][i][j]);
! R- o; _( g3 J& J7 uprintf("\n");
, w3 H) `- p+ p" [; I}
: r# ~, t, u' v# x  r9 c  mmicifang(a,h,2);
' A. E9 Q& J( g1 u9 R$ |* Z; Z2 Qfor(i=0;i<h;i++)
for(j=0;j<h;j++)
& u" i, w6 L7 z( k# za[h][i][j]=0; //全部赋值为0！
+ Z, q; |/ F# t9 rfor(i=0;i<h;i++)
for(j=0;j<h;j++)
for(t=0;t<h;t++)
) o! A$ q5 Z1 ]* N& A. j9 ka[h][i][j]=a[h][i][j]+a[t][i][j];
: P5 }' h; A; X+ o+ I: zprintf("\n\nB的值如下：\n");
for(i=0;i<h;i++)
{
for(j=0;j<h;j++)
printf(" %d ",a[h][i][j]);
printf("\n");
3 k# S9 N& E5 F) s; {: f8 m}
* `: [  [; W% A: T; x$ Zprintf("\n\n\n算法二：\n");
printf("\n矩阵A的1次布尔幂次方为:\n");
for(i=0;i<h;i++)
3 {% [* L7 ^0 Q: Q* L+ I{
! p: L7 g* O& G' Pfor(j=0;j<h;j++)
6 B& U0 t/ B$ B7 H) x{
& z. [" f% m3 k% g2 b9 o% Rif(a[0][i][j]>1)
a[0][i][j]=1;
printf(" %d ",a[0][i][j]);
}
printf("\n");
}
buer(a,h,2);
printf("\n\nP的值如下：\n");
5 V+ D! l0 R: b; _% G% ^1 L+ afor(i=0;i<h;i++)
  }/ ^) p/ V& O, t{
for(j=0;j<h;j++)
9 W, Q9 {7 J& x6 U( L; ~/ l& n{
if(a[h][i][j]>1)
1 Q  R$ A, }: g  N( G7 s8 p' r0 F5 r% V, Ha[h][i][j]=1;
  x% V& N- _" I9 wprintf(" %d ",a[h][i][j]);
}
printf("\n");
+ r8 F2 E3 d2 c- v* `+ z: h0 K}
printf("\n\n\nwarshall算法:\n\n");
' a8 j: D7 M( A% j4 Q
) L2 M8 J6 ?0 Rfor(i=0;i<h;i++) //处理邻接矩阵A.转化成可达矩阵。
for(j=0;j<h;j++)
if(a[0][i][j]!=0)
a[0][i][j]=1;
4 {  j8 t  n% g# k3 H0 iprintf("整理后的可达矩阵为：\n");
* A; O6 I6 x% ?2 v9 u6 }; ifor(i=0;i<h;i++)
{
8 g2 ^, s( B2 o( u1 i9 pfor(j=0;j<h;j++)
. }- b4 p3 i* y3 q8 kprintf(" %d ",a[0][i][j]);
printf("\n");
}
. N; N' J. n2 M6 e! N- owarshall(a,h);
0 }$ {4 e; l+ d" ugetch();
}
6 X! g2 o. H, D6 uint jisuan(int a[N][N][N],int h,int z,int i,int j)
! g3 ?& _( \  K* R5 q% D2 s/ W! v6 o: D{
int s,jg=0;
# q* H3 J2 p' a' q& uif(z==2)
% `% q' U4 q, m* @for(s=0;s<h;s++)
+ j! z5 D0 m$ D( l9 E* B5 Z8 Ejg=jg+a[z-2][i][s]*a[z-2][s][j];
else
for(s=0;s<h;s++)
/ b2 M- ]+ Z% a, w% ojg=jg+a[z-2][i][s]*a[0][s][j];
return jg;
}
  I" D0 c% P. q- j! p6 ivoid micifang(int a[N][N][N],int h,int t)
{
$ g) u8 L! M+ w$ |, x4 Lint i,j,s=0;
int z=t;
if(z<h) //矩阵的值不全为零。
6 \8 f3 H* [7 }/ W( ^. h5 T{
' ~, y* V# l  x) Y: }6 ^% h. Zfor(i=0;i<h;i++)
for(j=0;j<h;j++)
& l/ Y, }! K7 W3 Y8 E7 L* @a[z-1][i][j]=0; //全部赋值为0！
9 a1 W1 f! b2 g- U. K, ^! P2 m$ ffor(i=0;i<h;i++)
for(j=0;j<h;j++)
a[z-1][i][j]=jisuan(a,h,z,i,j);
1 _$ i, j" _+ Dprintf("\n矩阵A的%d次方为:\n",z);
7 ^8 y! r9 U6 D2 ^6 afor(i=0;i<h;i++)
{
  K& c+ p3 i5 z" D* @2 Qfor(j=0;j<h;j++)
  T1 H+ `" S; |printf(" %d ",a[z-1][i][j]);
+ ~3 R, G, q4 B( U+ gprintf("\n");
}
micifang(a,h,++z);
" x% x) ]2 K+ ?. x! \  J}
}
void buer(int a[N][N][N],int h,int t)
% k3 l" g5 y3 z{
int i,j,s=0;
$ k6 J& R0 `$ rint z=t;
- a0 r  ?+ S7 G4 l+ z$ T" C' l: bif(z<h) //矩阵的值不全为零。
{
for(i=0;i<h;i++)
) ?5 B; l% D. l  h# Hfor(j=0;j<h;j++)
- b# |3 T1 ^# m) @& Ga[z-1][i][j]=0; //全部赋值为0！
for(i=0;i<h;i++)
+ L6 s. V8 ^( i& ^for(j=0;j<h;j++)
# ]- O4 g+ F3 o6 v# j; N0 Xa[z-1][i][j]=jisuan(a,h,z,i,j);
printf("\n矩阵A的%d次布尔幂次方为:\n",z);
for(i=0;i<h;i++)
% {$ D$ v3 m. F{
for(j=0;j<h;j++)
0 E8 y/ h$ A7 B+ P{
if(a[z-1][i][j]>1)
a[z-1][i][j]=1;
% q6 k3 d$ j  N- b  T0 |( \printf(" %d ",a[z-1][i][j]);
9 ]- |& {6 v0 F$ @0 J) u}
printf("\n");
. }* x' ~1 Z& h/ Q5 X+ b}
buer(a,h,++z);
  `8 s+ l' `& ?. k3 d; A  L, E( @}
" }7 l& Z5 r  d& Z}
+ W$ X9 N5 F8 H8 tvoid warshall(int a[N][N][N],int h)
{
int i,j,p,q,r=0,s;
for(j=0;j<h;j++)
& h( n0 ]0 m2 C3 s2 hfor(i=0;i<h;i++)
& p2 v0 [( `+ W7 J' e. Qif(i!=j && a[0][i][j]==1)
{
2 L( K+ H* _1 ]% g  Efor(s=0;s<h;s++)
7 X8 g' |' Y( b5 F* R  X: ~$ Q  z{
1 h) h3 u0 C& a5 H) g' E7 ?1 ga[0][i][s]=a[0][i][s]+a[0][j][s];
( o1 K/ V. d/ `if(a[0][i][s]==2)
1 r( K' t+ L% ya[0][i][s]=1;
}
8 @( E# L% ~, w" ^8 g" Vprintf("\n第%d次运算的结果如下：\n",++r);
0 ^; \& E7 n8 d! [( w3 g/ lfor(p=0;p<h;p++)
  S: x" f& K: f' T+ G; ]{
* }1 [: w- w1 ?5 |for(q=0;q<h;q++)
printf(" %d ",a[0][p][q]);
' q4 \) q/ ~# W. Q* mprintf("\n");
}
}
5 P1 Q4 E+ k: u}

lizhaowei1990

先点一点这个是要干什么？