求教高手啊，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);
4 e) b1 K7 `; g: u1 t9 j$ k1 z  evoid micifang(int a[N][N][N],int h,int t);
  \/ _3 ?; ?. H  Nvoid buer(int a[N][N][N],int h);
void warshall(int a[N][N][N],int h);
) z) S; v9 }0 ^/ y  `- b- Pmain()
{
int i,h,j,t;
2 g, O; h) o. Aint a[N][N][N];
printf("请输入邻接矩阵A的行数(或列数)：\n");
scanf("%d",&h);
printf("\n请输入邻接矩阵A：\n\n");
2 J( m' C3 x$ N, Vfor(i=0;i<h;i++)
{
printf("请输入邻接矩阵A的第%d行.\n",i+1);
& Q7 ~8 O' Y* O2 f1 G# `) a# X- Ffor(j=0;j<h;j++)
scanf("%d",&a[0][i][j]);
}
$ H) a3 W9 m( j0 N1 k6 g; U& Hprintf("\n邻接矩阵A为：\n");
for(i=0;i<h;i++)
9 q- w& N# g. A9 O{
7 G( F0 f. @, r1 F* hfor(j=0;j<h;j++)
& a# T, D( v( X9 ]. Y, I5 Z2 rprintf(" %d ",a[0][i][j]);
/ ?, v& Y/ l( D. y, vprintf("\n");
  i' _. f* H+ i+ l/ }}
printf("\n\n\n算法一：\n");
printf("\n矩阵A的1次方为:\n");
( }4 h1 j+ g5 B1 n, ]  s" ~for(i=0;i<h;i++)
* X; I$ y: g" v( l{
9 {9 I$ j6 r; o& C3 |for(j=0;j<h;j++)
printf(" %d ",a[0][i][j]);
printf("\n");
}
* p8 p# K$ O9 ~, Z; o8 N; amicifang(a,h,2);
9 u3 d* o+ c8 H) ^% N; ^for(i=0;i<h;i++)
for(j=0;j<h;j++)
a[h][i][j]=0; //全部赋值为0！
for(i=0;i<h;i++)
/ \* T1 x& g5 J: n: J+ mfor(j=0;j<h;j++)
3 |. N/ Q" F+ [9 n7 i  `for(t=0;t<h;t++)
a[h][i][j]=a[h][i][j]+a[t][i][j];
printf("\n\nB的值如下：\n");
for(i=0;i<h;i++)
{
for(j=0;j<h;j++)
printf(" %d ",a[h][i][j]);
! j$ Y5 ~  }) w, P. U6 sprintf("\n");
}
printf("\n\n\n算法二：\n");
( @) L: t: n9 Y5 L, b) K1 Fprintf("\n矩阵A的1次布尔幂次方为:\n");
for(i=0;i<h;i++)
{
for(j=0;j<h;j++)
{
if(a[0][i][j]>1)
a[0][i][j]=1;
printf(" %d ",a[0][i][j]);
}
printf("\n");
}
buer(a,h,2);
+ _# W9 f4 b) }5 Oprintf("\n\nP的值如下：\n");
for(i=0;i<h;i++)
{
for(j=0;j<h;j++)
{
if(a[h][i][j]>1)
a[h][i][j]=1;
printf(" %d ",a[h][i][j]);
}
7 u3 l8 S$ B3 b& |printf("\n");
( D; h& v% D' Z; A}
# o) ?$ J' V' S5 D( Uprintf("\n\n\nwarshall算法:\n\n");
/ A8 J! n6 v2 Z
for(i=0;i<h;i++) //处理邻接矩阵A.转化成可达矩阵。
. {' t0 H6 v* b7 K* rfor(j=0;j<h;j++)
* k1 G2 R0 k7 r% b1 u. h' xif(a[0][i][j]!=0)
8 w3 r) \7 _. f: U4 l# Na[0][i][j]=1;
printf("整理后的可达矩阵为：\n");
. O9 j& h! k+ r$ h) q/ Z. ?* s) Hfor(i=0;i<h;i++)
$ {# B  }2 i4 B7 l* T{
: c9 L" m* n* z1 g2 xfor(j=0;j<h;j++)
4 n2 O; E2 R0 F2 C! [" I; Mprintf(" %d ",a[0][i][j]);
0 m9 W7 @  g$ w- O1 I0 _printf("\n");
}
8 q$ d( C4 V7 q) {5 u" L$ Cwarshall(a,h);
getch();
}
int jisuan(int a[N][N][N],int h,int z,int i,int j)
, h* _0 P6 q' _. C( z3 r{
int s,jg=0;
if(z==2)
8 u9 y8 l" W+ m/ R2 {0 zfor(s=0;s<h;s++)
jg=jg+a[z-2][i][s]*a[z-2][s][j];
else
3 q& B1 S; Y; v# a, m0 w& ffor(s=0;s<h;s++)
$ X" M$ K9 T3 g; m8 Q; Ljg=jg+a[z-2][i][s]*a[0][s][j];
return jg;
}
void micifang(int a[N][N][N],int h,int t)
{
int i,j,s=0;
* B* S* i2 q! A. A/ x! J" Oint z=t;
* N( ]" G: k! y' u2 K  gif(z<h) //矩阵的值不全为零。
{
* \7 v# V$ ]9 O9 A' x3 R: C% Afor(i=0;i<h;i++)
for(j=0;j<h;j++)
a[z-1][i][j]=0; //全部赋值为0！
for(i=0;i<h;i++)
! p! h, P/ j  C. O, lfor(j=0;j<h;j++)
a[z-1][i][j]=jisuan(a,h,z,i,j);
9 Z$ A* ^6 t5 h2 B& V, e4 q4 E" \6 ~printf("\n矩阵A的%d次方为:\n",z);
4 m7 C  F% E- W3 K" N* _( w8 _for(i=0;i<h;i++)
{
for(j=0;j<h;j++)
printf(" %d ",a[z-1][i][j]);
; l  u, ]6 ?7 E1 ]3 W3 x: nprintf("\n");
8 c0 u- M* G& Z8 U; m}
; j8 E" I, q: Omicifang(a,h,++z);
}
}
void buer(int a[N][N][N],int h,int t)
/ p4 ^; N: a) J/ E{
% b4 j+ M/ h7 ]$ I3 N$ R$ s% pint i,j,s=0;
* t) y$ j# t, O# ]* zint z=t;
if(z<h) //矩阵的值不全为零。
0 y3 W* t8 w4 {7 T{
for(i=0;i<h;i++)
for(j=0;j<h;j++)
; x3 K! c8 k2 t; }) `a[z-1][i][j]=0; //全部赋值为0！
7 N& W6 j3 m  D  Qfor(i=0;i<h;i++)
for(j=0;j<h;j++)
( H4 \+ F6 ?: m3 N- y! a. h1 V5 |a[z-1][i][j]=jisuan(a,h,z,i,j);
printf("\n矩阵A的%d次布尔幂次方为:\n",z);
for(i=0;i<h;i++)
{
8 D$ z& s. W* D. S! s0 hfor(j=0;j<h;j++)
{
0 j2 [: I# I) V) bif(a[z-1][i][j]>1)
* p& ?1 P, w2 B& Ga[z-1][i][j]=1;
printf(" %d ",a[z-1][i][j]);
# w9 Q8 _1 F% u8 s}
printf("\n");
}
buer(a,h,++z);
9 D/ K2 t) ^+ S# H1 g}
( V' K/ @. r( N) t}
void warshall(int a[N][N][N],int h)
{
int i,j,p,q,r=0,s;
for(j=0;j<h;j++)
6 N' M2 a- t4 ]# \2 J: g4 Q  H, Ifor(i=0;i<h;i++)
" ~. S1 T$ z; \; rif(i!=j && a[0][i][j]==1)
{
2 z, u4 g& s0 Dfor(s=0;s<h;s++)
{
  S9 d7 D$ `* s, g+ m7 Ya[0][i][s]=a[0][i][s]+a[0][j][s];
if(a[0][i][s]==2)
( C% L. p9 o3 U/ W1 o, xa[0][i][s]=1;
" |& `2 D% L7 w' `$ E}
! c; o6 \5 Y- Y6 Mprintf("\n第%d次运算的结果如下：\n",++r);
8 ^' g$ z5 Y% Vfor(p=0;p<h;p++)
* A0 b  ^  R- N6 w{
for(q=0;q<h;q++)
printf(" %d ",a[0][p][q]);
$ d# [: O0 F" ~+ r* S6 lprintf("\n");
}
; I+ B6 W8 f. F1 w4 w}
}

lizhaowei1990

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