mult array .cpp

改程序可以用于 实现数据结构中的对称矩阵相乘

//对称矩阵相乘的程序代码
#include<stdio.h>
#include<malloc.h>
//数组结构类型的定义.h
const int n=3;
const int size=n*(n+1)/2;
typedef int datatype;
typedef struct{
	datatype A[size],B[size],C[n][n];
}array;

void input(datatype[]);
void output(datatype[][n]);
void mult(array*);

void main()
{
	array*pa;
	pa=(array*)malloc(sizeof(array));
	printf("以行为主序输入矩阵A的下三角:\n");
    input(pa->A);//以行为主序输入矩阵A的下三角
	printf("以行为主序输入矩阵B的下三角:\n");
	input(pa->B);//以行为主序输入矩阵B的下三角
	mult(pa);
	output(pa->C);//输出矩阵C
}

//对称矩阵的输入
void input(datatype x[])
{
	for(int i=0;i<size;i++)
		scanf("%d",&x[i]);
}

//矩阵的输出
void output(datatype x[][n])
{
	for(int i=0;i<n;i++)
		{for(int j=0;j<n;j++)
			printf("%5d",x[i][j]); 
		printf("\n");
	}
}

//添加对称矩阵相乘算法
void mult(array* pa)
{	
	int a[n][n],b[n][n];
	int i=0,j=0,k=0;
    for(i=0;i<n;i++)     
		for(j=0;j<n;j++)
		{	a[i][j]=0;
            b[i][j]=0;
			pa->C[i][j]=0;                   
		}
	for(i=0;i<n;i++)      
		for(j=0;j<=i;j++)
		{	a[i][j]=pa->A[i*(i+1)/2+j];
			b[i][j]=pa->B[i*(i+1)/2+j];      
         }
	
		
    for(j=1;j<n;j++)
	  for(i=0;i<j;i++)
	  {a[i][j]=a[j][i];
		b[i][j]=b[j][i];
	  }
	pa->C[0][0]=0;       
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
		for(k=0;k<n;k++)
		{pa->C[i][j]+=a[i][k]*b[k][j];
		}
			
			
}