以下为展示“矩阵上三角为0值、下三角非0值”的简单示例:
1. 用c语言实现的版本
#include
#include
#define n 4
int t[n*n/2];
/* 把x存为L(i,j) */
void Store(int x, int i, int j)
{
if (i<0||j<0||i>=n||j>=n)
{
printf("数组出界!");
exit(1);
}
/* 当且仅当i ≥ j 时,(i,j)位于下三角*/
if (i >= j)
{
t[i*(i+1)/2+j] = x;
}
else if (x != 0)
{
printf("(上三角)非对角线上元素值必须为零");
exit(1);
}
}
void main()
{
int i,j;
int D[n][n]={2,0,0,0,5,1,0,0,0,3,1,0,4,2,7,0};
for(i=0;i
{
for(j=0;j
{
Store(D[i][j],i,j); /* 把x存为L(i,j) */
}
}
printf("输出原始矩阵:\n");
for(i=0; i
{
for(j=0; j
{
printf("%d ",D[i][j]);
}
printf("\n");
}
printf("诸行输出矩阵下三角的值:\n");
for(i=0;i<=n*n/2+1;i++)
{
printf("%d ",t[i]);
}
printf("\n");
}
运行结果如下图所示: