bo2-7.cpp 多项式(存储结构由c2-6.h定义)的基本操作及算法

1.  // bo2-7.cpp 多项式(存储结构由c2-6.h定义)的基本操作及算法2.22、2.23等(8个)
2.  #include"c2-5.h"
3.  #include"bo2-6.cpp"
4.  typedef LinkList polynomial;
5.  #define DestroyPolyn DestroyList // 与bo2-6.cpp中的函数同义不同名
6.  #define PolynLength ListLength // 与bo2-6.cpp中的函数同义不同名
8.  void OrderInsertMerge(LinkList &L,ElemType e,int(* compare)(term,term))
9.  { // 按有序判定函数compare()的约定，将值为e的结点插入或合并到升序链表L的适当位置
10.    Position q,s;
11.    if(LocateElem(L,e,q,compare)) // L中存在该指数项
12.    {
13.      q->data.coef+=e.coef; // 改变当前结点系数的值
14.      if(!q->data.coef) // 系数为0
15.      { // 删除多项式L中当前结点
16.        s=PriorPos(L,q); // s为当前结点的前驱
17.        if(!s) // q无前驱
18.          s=L.head;
19.        DelFirst(L,s,q);
20.        FreeNode(q);
21.      }
22.    }
23.    else // 生成该指数项并插入链表
24.    {
25.      MakeNode(s,e); // 生成结点
26.      InsFirst(L,q,s);
27.    }
28.  }
30.  int cmp(term a,term b) // CreatPolyn()的实参
31.  { // 依a的指数值<、=或>b的指数值，分别返回-1、0或+1
32.    if(a.expn==b.expn)
33.      return 0;
34.    else
35.      return (a.expn-b.expn)/abs(a.expn-b.expn);
36.  }
38.  void CreatPolyn(polynomial &P,int m) // 算法2.22
39.  { // 输入m项的系数和指数，建立表示一元多项式的有序链表P
40.    Position q,s;
41.    term e;
42.    int i;
43.    InitList(P);
44.    printf("请依次输入%d个系数，指数：/n",m);
45.    for(i=1;i<=m;++i)
46.    { // 依次输入m个非零项(可按任意顺序)
47.      scanf("%f,%d",&e.coef,&e.expn);
48.      if(!LocateElem(P,e,q,cmp)) // 当前链表中不存在该指数项，cmp是实参
49.      {
50.        MakeNode(s,e); // 生成结点并插入链表
51.        InsFirst(P,q,s);
52.      }
53.    }
54.  }
56.  void PrintPolyn(polynomial P)
57.  { // 打印输出一元多项式P
58.    Link q;
59.    q=P.head->next; // q指向第1个结点
60.    printf("  系数    指数/n");
61.    while(q)
62.    {
63.      printf("%f  %d/n",q->data.coef,q->data.expn);
64.      q=q->next;
65.    }
66.  }
68.  void AddPolyn(polynomial &Pa,polynomial &Pb) // 算法2.23
69.  { // 多项式加法：Pa=Pa+Pb，并销毁一元多项式Pb
70.    Position ha,hb,qa,qb;
71.    term a,b;
72.    ha=GetHead(Pa);
73.    hb=GetHead(Pb); // ha和hb分别指向Pa和Pb的头结点
74.    qa=NextPos(ha);
75.    qb=NextPos(hb); // qa和qb分别指向Pa和Pb中当前结点(现为第1个结点)
76.    while(!ListEmpty(Pa)&&!ListEmpty(Pb)&&qa)
77.    { // Pa和Pb均非空且ha没指向尾结点(qa!=0)
78.      a=GetCurElem(qa);
79.      b=GetCurElem(qb); // a和b为两表中当前比较元素
80.      switch(cmp(a,b))
81.      {
82.        case -1:ha=qa; // 多项式Pa中当前结点的指数值小
83.                qa=NextPos(ha); // ha和qa均向后移1个结点
84.                break;
85.        case 0: qa->data.coef+=qb->data.coef; // 两者的指数值相等，修改Pa当前结点的系数值
86.                if(qa->data.coef==0) // 删除多项式Pa中当前结点
87.                {
88.                  DelFirst(Pa,ha,qa);
89.          FreeNode(qa);
90.                }
91.                else
92.                  ha=qa;
93.                DelFirst(Pb,hb,qb);
94.                FreeNode(qb);
95.                qb=NextPos(hb);
96.                qa=NextPos(ha);
97.                break;
98.        case 1: DelFirst(Pb,hb,qb); // 多项式Pb中当前结点的指数值小
99.                InsFirst(Pa,ha,qb);
100.                ha=ha->next;
101.                qb=NextPos(hb);
102.      }
103.    }
104.    if(!ListEmpty(Pb))
105.    {
106.      Pb.tail=hb;
107.      Append(Pa,qb); // 链接Pb中剩余结点
108.    }
109.    DestroyPolyn(Pb); // 销毁Pb
110.  }
112.  void AddPolyn1(polynomial &Pa,polynomial &Pb)
113.  { // 另一种多项式加法的算法：Pa=Pa+Pb，并销毁一元多项式Pb
114.    Position qb;
115.    term b;
116.    qb=GetHead(Pb); // qb指向Pb的头结点
117.    qb=qb->next; // qb指向Pb的第1个结点
118.    while(qb)
119.    {
120.      b=GetCurElem(qb);
121.      OrderInsertMerge(Pa,b,cmp);
122.      qb=qb->next;
123.    }
124.    DestroyPolyn(Pb); // 销毁Pb
125.  }
127.  void Opposite(polynomial Pa)
128.  { // 一元多项式Pa系数取反
129.    Position p;
130.    p=Pa.head;
131.    while(p->next)
132.    {
133.      p=p->next;
134.      p->data.coef*=-1;
135.    }
136.  }
138.  void SubtractPolyn(polynomial &Pa,polynomial &Pb)
139.  { // 多项式减法：Pa=Pa-Pb，并销毁一元多项式Pb
140.    Opposite(Pb);
141.    AddPolyn(Pa,Pb);
142.  }
144.  void MultiplyPolyn(polynomial &Pa,polynomial &Pb)
145.  { // 多项式乘法：Pa=Pa×Pb，并销毁一元多项式Pb
146.    polynomial Pc;
147.    Position qa,qb;
148.    term a,b,c;
149.    InitList(Pc);
150.    qa=GetHead(Pa);
151.    qa=qa->next;
152.    while(qa)
153.    {
154.      a=GetCurElem(qa);
155.      qb=GetHead(Pb);
156.      qb=qb->next;
157.      while(qb)
158.      {
159.        b=GetCurElem(qb);
160.        c.coef=a.coef*b.coef;
161.        c.expn=a.expn+b.expn;
162.        OrderInsertMerge(Pc,c,cmp);
163.        qb=qb->next;
164.      }
165.      qa=qa->next;
166.    }
167.    DestroyPolyn(Pb); // 销毁Pb
168.    ClearList(Pa); // 将Pa重置为空表
169.    Pa.head=Pc.head;
170.    Pa.tail=Pc.tail;
171.    Pa.len=Pc.len;
172.  }