矩阵的压缩存储--特殊矩阵--三角矩阵

矩阵的压缩存储–特殊矩阵–三角矩阵

三角矩阵大体分为三类：对称矩阵、下三角矩阵和上三角矩阵。

#include <stdio.h>
#include <stdlib.h>
#define N 100
int a[N][N], sa[2*N];
int main(void)
{
    int i, j, k;
    int n, m;
    printf("请输入矩阵的行数和列数： ");
    scanf("%d%d", &n, &m);
    printf("请输入一个对称矩阵(a[i][j]=a[j][i])：\n");
    for(i = 1; i <= n; i++)
        for(j = 1; j <= m; j++)
            scanf("%d", &a[i][j]);
    for(i = 1; i <= n; i++){
        for(j = 1; j <= m; j++){
            if(i >= j)
                sa[i*(i-1)/2+j-1] = a[i][j];
        }
    }
    printf("当i>=j时，k=i(i-1)/2+j-1\n");
    printf("当i<j时，k=j(j-1)/2+i-1\n");
    printf("sa[]中的数据为： \n");
    for(i = 0; i <= n*(n-1)/2+m-1; i++)
        printf("sa[%d] = %d\n", i, sa[i]);
    printf("其对应关系为：\n");
    for(i = 1; i <= n; i++)
        for(j = 1; j <= m; j++)
            if(i >= j)
                printf("a[%d][%d] = sa[%d] = %d\n", i, j, i*(i-1)/2+j-1, sa[i*(i-1)/2+j-1]);
            else
                printf("a[%d][%d] = sa[%d] = %d\n", i, j, j*(j-1)/2+i-1, sa[j*(j-1)/2+i-1]);
    return 0;
}

示例运行结果