一个大型的稀疏矩阵要是采用三元存储，即数组RE，里面是矩阵的非零元素，数组ICOL，里面是非零元素的列号，这两个数组的个数是相同的，还有一个数组IRFR，里面是每一行第一个非零元素在ICOL的位置，

本人现在有一个大型的矩阵，要是想按照上述思想存储，怎么写程序啊，得到三个数组啊？