$求两个对称矩阵之和与乘积

#include <stdio.h>
#define X 3
#define Y 6

static int value(int a[], int i, int j)
{
if(i >= j)
        return a[(i * (i + 1)) / 2 + j];  //下三角 

    else
        return a[(j * (j + 1)) / 2 + i];  //上三角 
}

static void add(int a[], int b[], int c[][X])
{
    int i, j;

    for(i = 0; i < X; i++)
    {
        for(j = 0; j < X; j++)
        {
            c[i][j] = value(a, i, j) + value(b, i, j);
        }
    }
}

static void mul(int a[], int b[], int c[][X])
{
    int i, j, k, sum;

    for(i = 0; i < X; i++)
    {
        for(j = 0; j < X; j++)
        {
            sum = 0;
            for(k = 0; k < X; k++)
            {
                sum = sum + value(a, i, k) * value(b, k, j);
            }
            c[i][j] = sum;
        }
    }
}

static void dis1(int a[])
{
    int i, j;

    for(i = 0; i < X; i++)
    {
        for(j = 0; j < X; j++)
        {
            printf("%4d", value(a, i, j));
        }
        printf("\n");
    }
}

static void dis2(int c[][X])
{
    int i, j;

    for(i = 0; i < X; i++)
    {
        for(j = 0; j < X; j++)
        {
            printf("%4d", c[i][j]);
        }
        printf("\n");
    }
}


int main(int argc, char *argv[])
{
    int a[Y] = {1, 2, 3, 4, 5, 6}; 
    int b[Y] = {1, 1, 1, 1, 1, 1};
    int c1[X][X], c2[X][X];

    printf("A矩阵:\n");
    dis1(a);
    printf("B矩阵:\n");
    
    dis1(b);

    add(a, b, c1);
    printf("A+B:\n");
    dis2(c1);

    mul(a, b, c2);
    printf("A*B:\n");
    dis2(c2);

    int c[Y],d[Y],ad[X][X],s[X][X];
    printf("请输入C矩阵（例如A矩阵输入时应为1 2 3 4 5 6）:\n");
    for(int i=0; i<Y; i++)
        scanf("%d",c+i);        
    printf("C矩阵:\n");
    dis1(c);
    
    for(int i=0; i<Y; i++)
        scanf("%d",d+i);
    printf("D矩阵:\n");
    dis1(d);
    
    add(c, d, ad);
    printf("C+D:\n");
    dis2(ad);

    mul(c, d, s);
    printf("C*D:\n");
    dis2(s);

    return 0;
}