/*
*Copyright (c) 2014,煙臺大學計算機學院
*All rights reserved.
*檔名稱:Annpion.cpp
*作者:王耀鵬
*完成日期:2015年12月16日
*版本號:v1.0
*
*問題描述:上三角矩陣的壓縮儲存及基本運算
*輸入描述:輸入第n行的階層-n+1個數
*輸出描述:輸出上三角矩陣
*/
#include <stdio.h>
#include <malloc.h>
#define N 4
void Init(int *&b)//為N階對稱矩陣初始化儲存資料的一維陣列B
{
    b= (int *)malloc(sizeof(int )*(N*(N+1)/2+1));
}
int Value(int b[], int i, int j)//返回儲存在b[M]中，對應二維陣列A[i][j]的值
{
    if(j<i)
        return 0;
    else
        return b[j*(j+1)/2+i];
}
void Assign(int b[], int e, int i, int j)//將e賦值給對應二維陣列元素A[i][j]，要儲存到B[M]中
{
    if(j<i)
        b[N*(N+1)/2]=e;
    else
        b[j*(j+1)/2+i]=e;
}
void Disp(int b[])//輸出壓縮儲存在b中的對稱矩陣
{
    for(int i=0;i<N;++i)
    {
        for(int j=0;j<N;++j)
            printf("%4d",Value(b,i,j));
        printf("\n");
    }
}
void Destroy(int b[])//銷燬儲存空間
{
    free(b);
}
int main()
{
    int *b1;  //指向整型的指標，待初始化
    int i, j;
    int v;
    Init(b1);
    printf("請輸入對稱矩陣(只需要輸入下三角部分即可)\n");
    for(i=0;i<N;i++)
    {
        printf("輸入第%d行的%d個數據元素: ", i+1, N-i);
        for(j=i; j<N; j++)
        {
             scanf("%d", &v);
             Assign(b1, v, i, j);
        }
    }
    Disp(b1);
    Destroy(b1);
    return 0;
}
 

執行結果：