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

chenhande

#include "stdio.h"
: H! g/ p# Q; i* C! b, j, T#include "string.h"
#define N 50
int jisuan(int a[N][N][N],int h,int z,int i,int j);
& B; k! d4 x: f8 A# U5 X7 Bvoid micifang(int a[N][N][N],int h,int t);
void buer(int a[N][N][N],int h);
void warshall(int a[N][N][N],int h);
( @5 Q2 \0 f" W% @main()
( R0 b1 ]  Y& n5 v# o- i( o{
int i,h,j,t;
int a[N][N][N];
printf("请输入邻接矩阵A的行数(或列数)：\n");
scanf("%d",&h);
. N6 i2 C0 d4 M9 y, o. jprintf("\n请输入邻接矩阵A：\n\n");
) o0 O0 ?3 i+ B0 e% w6 B& yfor(i=0;i<h;i++)
{
printf("请输入邻接矩阵A的第%d行.\n",i+1);
for(j=0;j<h;j++)
8 f$ \3 i* B% f5 D+ u; {! F" o/ F0 @scanf("%d",&a[0][i][j]);
1 I- [! g3 W& @+ k( {}
printf("\n邻接矩阵A为：\n");
for(i=0;i<h;i++)
& z7 x: S- |" I6 t6 \{
- Y, t$ ?; `# U# c; H6 g9 G, e' Sfor(j=0;j<h;j++)
printf(" %d ",a[0][i][j]);
$ e9 B) V; L$ y0 Pprintf("\n");
1 t/ d6 I) y1 l8 q/ I$ U% U6 E1 ]}
, O+ F" k5 ~  Q1 m6 t! S+ jprintf("\n\n\n算法一：\n");
1 P  B+ o1 s  U2 j* R1 x1 ^/ Oprintf("\n矩阵A的1次方为:\n");
! p/ D1 V4 K0 H! M1 Ofor(i=0;i<h;i++)
8 H6 Y2 `2 E( ~: ?# q" Y{
for(j=0;j<h;j++)
$ H, y$ c6 w: g  e/ lprintf(" %d ",a[0][i][j]);
printf("\n");
8 z6 V4 A) J0 D% q4 |}
micifang(a,h,2);
for(i=0;i<h;i++)
for(j=0;j<h;j++)
8 p$ D: b; B) r2 \a[h][i][j]=0; //全部赋值为0！
7 N7 l0 C9 A- A7 Y7 M) efor(i=0;i<h;i++)
for(j=0;j<h;j++)
for(t=0;t<h;t++)
a[h][i][j]=a[h][i][j]+a[t][i][j];
printf("\n\nB的值如下：\n");
  @6 J, g8 n0 O0 S8 Hfor(i=0;i<h;i++)
' P; C3 s  D1 a% E" @{
$ b: F3 N# `: ^# R! ~* W+ \for(j=0;j<h;j++)
printf(" %d ",a[h][i][j]);
printf("\n");
! C9 w4 g, i" g9 x- [2 y8 i2 r0 L: F4 l}
" y: _& y8 P% D0 x% `printf("\n\n\n算法二：\n");
printf("\n矩阵A的1次布尔幂次方为:\n");
for(i=0;i<h;i++)
{
5 A- J' K. E2 L; T" x5 }* Cfor(j=0;j<h;j++)
9 {$ l; x& y/ X: S{
if(a[0][i][j]>1)
a[0][i][j]=1;
% p! j1 `5 u/ Z  w2 [% rprintf(" %d ",a[0][i][j]);
- b5 L# r9 V% V0 g7 J% m}
printf("\n");
4 M# _% g1 a( f% e. A  n8 |}
buer(a,h,2);
/ A" x! K8 i4 iprintf("\n\nP的值如下：\n");
: n! |% S7 N: R. p7 O+ F% h+ m0 O! f4 Qfor(i=0;i<h;i++)
{
for(j=0;j<h;j++)
{
if(a[h][i][j]>1)
6 Z* H, E- [% E5 N. H/ l; @9 ba[h][i][j]=1;
printf(" %d ",a[h][i][j]);
}
/ E2 J* T$ p# I( E0 I7 V& dprintf("\n");
}
printf("\n\n\nwarshall算法:\n\n");

for(i=0;i<h;i++) //处理邻接矩阵A.转化成可达矩阵。
9 ~7 b3 t; J& @# x" \8 d' |# \for(j=0;j<h;j++)
if(a[0][i][j]!=0)
) y& Z) S$ {* t4 P, d1 ta[0][i][j]=1;
printf("整理后的可达矩阵为：\n");
: m- ^8 V% X" Zfor(i=0;i<h;i++)
{
for(j=0;j<h;j++)
printf(" %d ",a[0][i][j]);
printf("\n");
" w9 Q: F4 r. X4 C" \5 P}
warshall(a,h);
getch();
}
int jisuan(int a[N][N][N],int h,int z,int i,int j)
% L/ L0 g0 ?& w5 @2 V; j/ f{
" M7 a- J$ j  pint s,jg=0;
* r9 u9 e& v* _! [  N: Y; Gif(z==2)
for(s=0;s<h;s++)
+ K% k1 I, }  L: c. U1 N1 Ijg=jg+a[z-2][i][s]*a[z-2][s][j];
else
7 {( P2 x1 M6 V0 d+ f! A8 C8 ?# Gfor(s=0;s<h;s++)
jg=jg+a[z-2][i][s]*a[0][s][j];
) P9 Z) v, F6 D8 Dreturn jg;
}
- Y7 h9 i. ?6 H9 e( a& jvoid micifang(int a[N][N][N],int h,int t)
{
int i,j,s=0;
  q) R& f& P& s. Z- \3 z- L) Mint z=t;
7 u% g4 R+ ]/ _( W* ]' X$ X; W; kif(z<h) //矩阵的值不全为零。
" R! s4 H+ t+ ~1 P, }{
$ Y2 ^8 p" n7 g& Kfor(i=0;i<h;i++)
. y+ n# c2 M+ ]" ]for(j=0;j<h;j++)
0 e+ s+ b$ ]9 _9 z' F, o/ \a[z-1][i][j]=0; //全部赋值为0！
: v& s, W) p+ I& i4 H5 U( Hfor(i=0;i<h;i++)
4 ?; s  ?: _' w" D2 O" sfor(j=0;j<h;j++)
a[z-1][i][j]=jisuan(a,h,z,i,j);
printf("\n矩阵A的%d次方为:\n",z);
9 O9 O& N8 a4 b" n; afor(i=0;i<h;i++)
' u+ x9 L! l8 ?' y5 e{
3 G( [7 g: e5 q- e1 A6 l& Ofor(j=0;j<h;j++)
# u: R5 x3 C8 I+ O6 V" X3 uprintf(" %d ",a[z-1][i][j]);
  i5 y+ f, H2 ]/ q+ ]+ f7 K+ uprintf("\n");
; X$ j  V+ K7 C- F9 W3 e  `, b1 `+ a. }}
6 S4 |" l( Y% g7 S4 Xmicifang(a,h,++z);
}
' o8 L( {9 X2 J; i}
3 @1 e" W0 T- l. C8 j9 D$ gvoid buer(int a[N][N][N],int h,int t)
{
7 c* K( M/ b/ Q/ |int i,j,s=0;
9 u2 q. ]- ^% r. }int z=t;
if(z<h) //矩阵的值不全为零。
{
for(i=0;i<h;i++)
" F- \# s* M, ?for(j=0;j<h;j++)
a[z-1][i][j]=0; //全部赋值为0！
  e3 P+ {9 ]$ Pfor(i=0;i<h;i++)
for(j=0;j<h;j++)
a[z-1][i][j]=jisuan(a,h,z,i,j);
# j' q) w; b( p; f* i& |printf("\n矩阵A的%d次布尔幂次方为:\n",z);
for(i=0;i<h;i++)
{
* {, l, Z' z/ P8 B9 e2 [) ^2 |for(j=0;j<h;j++)
{
if(a[z-1][i][j]>1)
  c4 a2 v$ s; C) `7 P; |5 q+ f7 Ya[z-1][i][j]=1;
2 F+ @" \; w6 {. Eprintf(" %d ",a[z-1][i][j]);
}
4 L) s2 y& E* P3 k7 pprintf("\n");
$ O: ^$ F) k7 M8 f/ Y' G7 ?# p}
# p+ c1 z) L2 V7 K( p3 x9 z+ O/ mbuer(a,h,++z);
}
5 b$ m$ U3 t6 i% |" P! S}
void warshall(int a[N][N][N],int h)
0 q0 i1 C7 y1 h9 r2 G% }3 X{
+ T1 e" i' ^* C& {int i,j,p,q,r=0,s;
* E, Q* l$ ]  Rfor(j=0;j<h;j++)
6 D# }2 ?! v; I3 ?for(i=0;i<h;i++)
if(i!=j && a[0][i][j]==1)
{
8 W8 Z6 f  N& X0 i% C" d) qfor(s=0;s<h;s++)
{
a[0][i][s]=a[0][i][s]+a[0][j][s];
if(a[0][i][s]==2)
1 e$ c/ o  Z8 T: M8 Ma[0][i][s]=1;
}
& o. p$ C, ~1 `5 ?. M0 o2 Wprintf("\n第%d次运算的结果如下：\n",++r);
for(p=0;p<h;p++)
. n/ l4 {3 g4 [. I; h1 Z' t& r" W{
for(q=0;q<h;q++)
printf(" %d ",a[0][p][q]);
printf("\n");
/ [- ^+ g* N) Y7 w# k}
  N8 V6 f$ B5 }7 c- l; _9 c+ Q}
}

lizhaowei1990

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