若n阶矩阵A中的元满足下述性质(以数组表示)
A[i][j] = A[j][i] 0=<i,j<=n-1
则成为n阶对称矩阵。
对于n阶矩阵,我们可以为每一对对称元分配一个存储空间,可以将n*n个元压碎存储到n(n+1)/2个元的空间中。
不失一般性,我们可以以行序为主序存储其下三角(包括对角线)中的元。
假设以一位数组sa[n(n+1)/2]作为n阶对称矩阵A的存储结构,则sa[k]和矩阵数组A[i][j]之间存在着一一对应的关系:
k=i(i+1)/2+j 当i>=j时
k=j(j+1)/2+i 当i<j 时
#include"stdio.h"
#define n 4
int t[n*(n+1)/2];
void Store(int x,int i,int j)
{
if(x<0||j<0||i>=n||j>=n)
{
printf("数组出界!");
exit(1);
}
else if(i>=j)
t[i*(i+1)/2+j]=x;
}
void main()
{
int i,j;
int D[n][n]={2,4,6,0,4,1,9,5,6,9,4,7,0,5,7,0};
for(i=0;i<n;i++)
for(j=0;j<n;j++)
Store(D[i][j],i,j);
for(i=0;i<=n*n/2+1;i++)
printf("%d ",t[i]);
printf("\n");
}