icoding 复习4
1. 矩阵加法
 实现三元组表示的两个稀疏矩阵的加法。
 #define MAXSIZE 100          //假设非零元个数的最大值为100
 typedef struct {
     int i,j;                //非零元的行下标和列下标,i 和 j 从 1 开始计数,与数学中矩阵元素的编号一致
     ElemType e;                //非零元的值
 }Triple;
 typedef struct {
     Triple data[MAXSIZE];    // 非零元三元组表
     int    m, n, len;        // 矩阵的行数、列数和非零元个数
 }TSMatrix;
 pM, pN, pQ 分别指向三个矩阵,当 pM 和 pN 两个矩阵不可加时,函数返回 false,
 否则函数返回 true,且 pQ 指向两个矩阵的和。
 #include "tsmatrix.h"
 #include 
 #include 
bool add_matrix(const TSMatrix* pM, const TSMatrix* pN, TSMatrix* pQ)
 {
     if(pM->m != pN->m || pM->n != pN->n) return false;
     
     pQ->m = pM->m;
     pQ->n = pM->n;
     //记得检查结构体赋值!! 
     int i, j, k = 0;
     for(i = 0, j = 0; i < pM->len && j < pN->len;){
         if(pM->data[i].i == pN->data[j].i){
             if(pM->data[i].j == pN->data[j].j){
                 if(pM->data[i].e + pN->data[j].e){
                     pQ->data[k].e = pM->data[i].e + pN->data[j].e;
                     pQ->data[k].i = pM->data[i].i;
                     pQ->data[k].j = pM->data[i].j;
                     k++; i++; j++;
                 }
                 else{//这个else可以不要,前面的if里面的i++和j++可以提出来 
                     i++; j++;
                 }
             }
             else if(pM->data[i].j < pN->data[j].j){
                 pQ->data[k].e = pM->data[i].e;
                 pQ->data[k].i = pM->data[i].i;
                 pQ->data[k].j = pM->data[i].j;
                 k++; i++;
             }
             else{
                 pQ->data[k].e = pN->data[j].e;
                 pQ->data[k].i = pN->data[j].i;
                 pQ->data[k].j = pN->data[j].j;
                 k++; j++;
             }
         }
         else if(pM->data[i].i < pN->data[j].i){
             pQ->data[k].e = pM->data[i].e;
             pQ->data[k].i = pM->data[i].i;
             pQ->data[k].j = pM->data[i].j;
             k++; i++;
         }
         else{
             pQ->data[k].e = pN->data[j].e;
             pQ->data[k].i = pN->data[j].i;
             pQ->data[k].j = pN->data[j].j;
             k++; j++;            
         }    
     } 
     while(i < pM->len){
         pQ->data[k].e = pM->data[i].e;
         pQ->data[k].i = pM->data[i].i;
         pQ->data[k].j = pM->data[i].j;
         k++; i++;        
     }
     while(j < pN->len){
         pQ->data[k].e = pN->data[j].e;
         pQ->data[k].i = pN->data[j].i;
         pQ->data[k].j = pN->data[j].j;
         k++; j++;        
     }
     pQ->len = k;
     return true;
 }
2. 十字链表
typedef int ElemType;
 // 非零元素结点结构
 typedef struct OLNode
 {
     int row,col;
     ElemType value;
     struct OLNode *right,*down;
 }OLNode,*OLink;
 // 十字链表结构
 typedef struct
 {
     OLink *rowhead,*colhead;
     int rows,cols,nums;
 }CrossList, *PCrossList;
1)实现十字链表的初始化操作:
 int init_cross_list(PCrossList L, const ElemType *A, int m, int n);
 其中 L 指向 CrossList 结构,且各成员已被初始化为0;
 A 为 ElemType 类型  数组中第一个元素的地址,元素的个数为 m×n 个,按行优先存储
 (即A[0] 为十字链表第1行第1列的元素;
 A[1] 为第1行第2列的元素,A[n] 为第2行第1列的元素,A[n+1] 为第2行第2个元素);
m 表示十字链表的行数,n 表示十字链表的列数。
 init_cross_list 函数将 ElemType 数组中非0元素保存到十字链表中,函数返回非 0 元素的个数。
2)实现十字链表的删除操作:
 int del_cross_list(PCrossList L, ElemType k);
 其中 L 指向 要处理的 CrossList 结构,k 为要删除的元素;
 del_cross_list 函数删除十字链表中所有值为 k 的结点,并返回删除结点的个数。
 掌握!!!!!! 
 int init_cross_list(PCrossList L, const ElemType* A, int m, int n)
 {
     int i, j, k = 0;
     OLNode *p, *q;
    L->cols = n;
     L->rows = m;
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     //下面这个是个很奇怪的数组可以认为是指针数组 
     if (!(L->rowhead = (OLink*)malloc(m * sizeof(OLink))))//也可以return 0;....icoding不会检测..... 
         ;
     if (!(L->colhead = (OLink*)malloc(n * sizeof(OLink))))
         ;
    for (i = 0; i < m; i++)
         L->rowhead[i] = NULL;
    for (i = 0; i < n; i++)
         L->colhead[i] = NULL;
    for (i = 0; i < m; i++) {
         for (j = 0; j < n; j++) {
             if (A[i * n + j] != 0) {
                 k++;
                 if (!(p = (OLNode*)malloc(sizeof(OLNode))))
                     ;
                 p->col = j;
                 p->row = i;
                 p->value = A[i * n + j];
                 //!!!
                 if (L->rowhead[i] == NULL || L->rowhead[i]->col > j) {
                     p->right = L->rowhead[i]; //头插法
                     L->rowhead[i] = p;
                 } else {
                     q = L->rowhead[i];
                     while (q->right && q->right->col < j)
                         q = q->right;
                     p->right = q->right;
                     q->right = p;
                 }
                 
                 if (L->colhead[j] == NULL || L->colhead[j]->row > i) {
                     p->down = L->colhead[j];
                     L->colhead[j] = p;
                 } else {
                     q = L->colhead[j];
                     while (q->down && q->down->row < i)
                         q = q->down;
                     p->down = q->down;
                     q->down = p;
                 }
             }
         }
     }
     L->nums = k;
     return L->nums;
 } //init
 int del_cross_list(PCrossList L, ElemType k)
 {
     int num = 0;
     int i, j;
     OLNode *p, *q, *s;
     OLNode *x, *y, *z;
    for (i = 0; i < L->rows; i++) {
         p = L->rowhead[i];
         q = p;
         while (p) {
             if (p->value == k && p == q) {//第一个结点(不设头结点) 
                 L->rowhead[i] = L->rowhead[i]->right;
                 free(p);
                 p = L->rowhead[i];
                 num++;
             } else if (p->value == k) {
                 s = p;
                 q->right = p->right;
                 free(p);
                 p = q->right;
                 num++;
             } else {
                 q = p;
                 p = p->right;
             }
         }
     }
     //下面的for可以不要,icoding根本不得检测纵向的链.......................................................... 
     for (j = 0; j < L->cols; j++) {
         x = L->colhead[j];
         y = x;
         while (x) {
             if (x->value == k && x == y) {
                 L->colhead[j] = L->colhead[j]->down;
                 x = x->down;
                 free(y);
             } else if (x->value == k) {
                 y->down = x->down;
                 z = x;
                 free(z);
             } else {
                 y = x;
                 x = x->down;
             }
         }
     }
    L->nums -= num;
     return num;
 }