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

chenhande

#include "stdio.h"
#include "string.h"
#define N 50
int jisuan(int a[N][N][N],int h,int z,int i,int j);
void micifang(int a[N][N][N],int h,int t);
4 p6 q* U. |( o" i* m& l3 Kvoid buer(int a[N][N][N],int h);
1 x9 f* `, K3 Avoid warshall(int a[N][N][N],int h);
main()
( h+ T+ G6 b' M2 f/ Z, q. a7 M& @) ~/ x{
int i,h,j,t;
& m8 |3 X% M8 |$ z. v# q; Bint a[N][N][N];
0 {& W' \% e6 ?5 ~0 q; Tprintf("请输入邻接矩阵A的行数(或列数)：\n");
scanf("%d",&h);
printf("\n请输入邻接矩阵A：\n\n");
: B: ]0 K$ Z  C4 T6 T" W$ U: Ufor(i=0;i<h;i++)
1 c2 ~/ c: d* K4 w$ c- O$ `{
% }0 R. W- ?3 k' V6 C7 Yprintf("请输入邻接矩阵A的第%d行.\n",i+1);
6 a" `! Z; M8 I5 Lfor(j=0;j<h;j++)
scanf("%d",&a[0][i][j]);
}
  k2 k+ Q4 `+ d9 |5 R% `printf("\n邻接矩阵A为：\n");
for(i=0;i<h;i++)
{
for(j=0;j<h;j++)
printf(" %d ",a[0][i][j]);
printf("\n");
}
printf("\n\n\n算法一：\n");
printf("\n矩阵A的1次方为:\n");
5 A3 `: x+ w+ O3 S/ ]for(i=0;i<h;i++)
* f8 B% m. ~, H* H{
for(j=0;j<h;j++)
printf(" %d ",a[0][i][j]);
printf("\n");
}
micifang(a,h,2);
6 p0 I% o: P) N# V1 _for(i=0;i<h;i++)
for(j=0;j<h;j++)
6 {' h! \% k; U  |+ x3 Wa[h][i][j]=0; //全部赋值为0！
for(i=0;i<h;i++)
for(j=0;j<h;j++)
" @5 W- x$ C4 xfor(t=0;t<h;t++)
a[h][i][j]=a[h][i][j]+a[t][i][j];
2 ?3 t1 g3 z& B0 K  wprintf("\n\nB的值如下：\n");
for(i=0;i<h;i++)
1 s& e0 @6 O8 ~; e/ {  J- P0 `  e: }{
, q* T0 ?! a; [, p. D3 t/ Q0 Qfor(j=0;j<h;j++)
printf(" %d ",a[h][i][j]);
9 T/ I; U  T+ D" H) Iprintf("\n");
}
printf("\n\n\n算法二：\n");
printf("\n矩阵A的1次布尔幂次方为:\n");
% |" s( o0 f9 @* Q3 ufor(i=0;i<h;i++)
{
6 E) P  J7 b( f# A. pfor(j=0;j<h;j++)
& i+ _2 Z% B) \{
; A1 P1 I0 e5 c/ k) g* aif(a[0][i][j]>1)
a[0][i][j]=1;
% f8 D9 C4 y6 M9 eprintf(" %d ",a[0][i][j]);
}
printf("\n");
}
buer(a,h,2);
+ W% d; S1 C* S3 q! {+ d' uprintf("\n\nP的值如下：\n");
8 u  D! j, S7 m( A* Jfor(i=0;i<h;i++)
; D7 S& }8 G' a( W& b) |{
for(j=0;j<h;j++)
{
: g+ z$ m( e- R8 \8 Y- ]: X- fif(a[h][i][j]>1)
1 s6 t% O. v9 |a[h][i][j]=1;
printf(" %d ",a[h][i][j]);
8 u7 Z, [1 U& }( @  S}
printf("\n");
* J% g& N) N  d- p7 V}
printf("\n\n\nwarshall算法:\n\n");

for(i=0;i<h;i++) //处理邻接矩阵A.转化成可达矩阵。
for(j=0;j<h;j++)
if(a[0][i][j]!=0)
7 e. l# z1 u: {8 Ta[0][i][j]=1;
printf("整理后的可达矩阵为：\n");
for(i=0;i<h;i++)
* s' R' o6 }: ]+ ~; P& p{
: x" o! m" L3 @7 r$ B0 @for(j=0;j<h;j++)
printf(" %d ",a[0][i][j]);
printf("\n");
}
+ p5 E' U% W7 J: o. o/ v/ n% Zwarshall(a,h);
getch();
}
9 b6 \9 p/ `* `7 C4 C. Wint jisuan(int a[N][N][N],int h,int z,int i,int j)
{
- y) g2 _* N, p6 uint s,jg=0;
if(z==2)
9 F8 S% O% n9 n0 m: D7 X' c, M# W4 afor(s=0;s<h;s++)
jg=jg+a[z-2][i][s]*a[z-2][s][j];
else
for(s=0;s<h;s++)
jg=jg+a[z-2][i][s]*a[0][s][j];
, h% ]4 v5 b( Kreturn jg;
# |) Z, a. ]- [& r8 m, Z}
void micifang(int a[N][N][N],int h,int t)
: ?# M! v* k. g1 z/ Z{
int i,j,s=0;
int z=t;
( e1 B- M& K  d0 x4 s, D  K: hif(z<h) //矩阵的值不全为零。
{
for(i=0;i<h;i++)
for(j=0;j<h;j++)
a[z-1][i][j]=0; //全部赋值为0！
- @  H+ Y. H4 h/ U/ Lfor(i=0;i<h;i++)
4 Z- @( i) |% Afor(j=0;j<h;j++)
  a/ t4 K  o8 X* Ta[z-1][i][j]=jisuan(a,h,z,i,j);
6 l- a' H9 \' G) Vprintf("\n矩阵A的%d次方为:\n",z);
for(i=0;i<h;i++)
{
& R8 {3 }, T) n/ I& H$ \' Hfor(j=0;j<h;j++)
printf(" %d ",a[z-1][i][j]);
printf("\n");
. i+ E7 N; \8 s0 w. J8 f}
micifang(a,h,++z);
) d7 f# L7 q7 V# W! |}
. w' \* J# Y$ I  d/ e& f. P1 ?; D}
/ V1 t# T; I$ m# C" F: A9 Z: tvoid buer(int a[N][N][N],int h,int t)
{
int i,j,s=0;
int z=t;
if(z<h) //矩阵的值不全为零。
{
for(i=0;i<h;i++)
8 q& l0 @/ G( M; H; N$ {- L& ]for(j=0;j<h;j++)
% x6 [. _) Q3 z( F( `0 p- K4 v/ R' aa[z-1][i][j]=0; //全部赋值为0！
* n9 |: n7 Y( H& g. H- Zfor(i=0;i<h;i++)
# N, P. T7 a9 Q  z( sfor(j=0;j<h;j++)
+ _" j! v* a$ Na[z-1][i][j]=jisuan(a,h,z,i,j);
# E* k  E6 \) |$ O$ e) qprintf("\n矩阵A的%d次布尔幂次方为:\n",z);
for(i=0;i<h;i++)
{
. Y0 y  Z5 c8 D! e$ d, Pfor(j=0;j<h;j++)
{
, H  ]! f4 P( \' h$ sif(a[z-1][i][j]>1)
a[z-1][i][j]=1;
1 a5 d, [0 w: ~printf(" %d ",a[z-1][i][j]);
3 r0 C$ k! Z6 S  F* b}
0 N/ X2 D! f$ R9 t! cprintf("\n");
}
buer(a,h,++z);
$ L) Y8 t$ ^+ c- a/ q}
}
void warshall(int a[N][N][N],int h)
. U" ]2 K# V5 y+ ^: q; Z8 @{
int i,j,p,q,r=0,s;
- T' E' F/ ^( ^$ p' b& Q1 [for(j=0;j<h;j++)
for(i=0;i<h;i++)
if(i!=j && a[0][i][j]==1)
{
% z+ R( Y' o! V' L* q# B: pfor(s=0;s<h;s++)
{
a[0][i][s]=a[0][i][s]+a[0][j][s];
if(a[0][i][s]==2)
' g' @0 _+ k. }+ q6 O$ c' z* U$ Oa[0][i][s]=1;
/ D4 o: r& t' v/ r+ {}
* Q9 k& N8 y" h% L1 D, J/ J$ |, s! yprintf("\n第%d次运算的结果如下：\n",++r);
for(p=0;p<h;p++)
. M2 T) u' Z# y, M$ E, L( {" L) T$ l9 \{
. M- W$ H9 h* q; M+ mfor(q=0;q<h;q++)
printf(" %d ",a[0][p][q]);
printf("\n");
6 U) W) c. x# ?}
: P, D; o; _* k# z}
}

lizhaowei1990

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