实验七 对称矩3

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实验七对称矩阵

#include

int C[100];

int D[100][100];

void show(int a[],int n)

{

int i,j;

for(i=0;i

{

for(j=0;j

{

if(i>=j)

printf("%d\t",a[i*(i+1)/2+j]);

else

printf("%d\t",a[j*(j+1)/2+i]);

}

printf("\n");

}

printf("\n");

}

void addA_B(int a[],int b[],int n) {

int i;

for(i=0;i

{

C[i]=a[i]+b[i]; }

}

void AchengB(int a[],int b[],int n) {

int k,l,i,sum=0;

for(k=0;k

{

for(l=0;l

{

for(i=0;i

{

if(k>=i&&i>=l) sum=sum+a[k*(k+1)/2+i]*b[i*(i+ 1)/2+l];

if(k=l)

sum=sum+a[i*(i+1)/2+k]*b[i*(i+1 )/2+l];

if(k>=i&&i

sum=sum+a[k*(k+1)/2+i]*b[l*(l+ 1)/2+i];

if(k

sum=sum+a[i*(i+1)/2+k]*b[l*(l+1 )/2+i];

}

D[k][l]=sum;

sum=0;

}

}

void main()

{

int k,l;

int

A[10]={1,2,3,4,5,6,7,8,9,10};

int

B[10]={1,2,3,4,5,6,7,8,9,10};

printf("A:\n");

show(A,4);

printf("B:\n");

show(B,4);

addA_B(A,B,4);

printf("A+B:\n");

show(C,4);

AchengB(A,B,4);

printf("A*B:\n");

for(k=0;k<4;k++)

{

for(l=0;l<4;l++)

{

printf("%d\t",D[k][l]);

}

printf("\n");

}

}

相关文档
最新文档