多项式加法

1. /*多项式(Polyomial)加法
2.     code by zzr 
3. http://blog.youkuaiyun.com/ZhengZhiRen/   */
4. #include <stdio.h>
5. #include <malloc.h>
6. #include <stdlib.h>
7. #define LEN sizeof(node)
10. typedef   struct  node
11. {
12.      float  coef;      /*系数(coefficient)*/
13.      int  exp;         /*指数(exponent)*/
14.      struct  node *next;
15. }node;
18. node * create_poly()     /*建立多项式*/
19. {
20.     node *head;
21.     node *p;
22.     node *pt;
24.     head=(node *)malloc(LEN);    /*头结点*/
25.     p=head;
26.     pt=(node *)malloc(LEN);
27.     scanf( "%f%d" ,&pt->coef,&pt->exp);
28.      while (pt->coef!=0)
29.     {
30.         p->next=pt;
31.         p=pt;
32.         pt=(node *)malloc(LEN);
33.         scanf( "%f%d" ,&pt->coef,&pt->exp);
34.     }
35.     p->next=NULL;
37.      return  head;
38. }
41. int  cmp(node *pa,node *pb)
42. {
43. /*比较两个项的指数大小*/
44.      if (!pa||!pb)
45.          return  -2;   /*返回错误*/
46.      if (pa->exp>pb->exp)
47.          return  1;
48.      else   if (pa->exp==pb->exp)
49.          return  0;
50.      else
51.          return  -1;
52. }
55. void  add_poly(node *ha,node *hb)         /*两个多项式相加*/
56. {
57.     node *pa;
58.     node *pb;
59.     node *pt;
61.     pa=ha->next;
62.     pb=hb->next;
63.      while (pa&&pb)
64.     {
65.          switch (cmp(pa,pb))
66.         {
67.          case  1:      /*pa指数大*/
68.             pt=pb;
69.             hb->next=pt->next;
70.             pt->next=pa;
71.             ha->next=pt;
72.              /*指针后移*/
73.             ha=ha->next;
74.             pa=ha->next;
75.             pb=hb->next;
76.              break ;
77.          case  0:      /*指数相等*/
78.             pa->coef+=pb->coef;
79.             pt=pb;
80.              /*指针后移*/
81.             ha=ha->next;
82.             pa=ha->next;
83.             hb->next=pb->next;
84.             pb=hb->next;
85.             free(pt);
86.              break ;
87.          case  -1:     /*pb指数大*/
88.             ha=pa;
89.             pa=ha->next;
90.              break ;
91.         }
92.     }
93.      if (NULL==pa)
94.         ha->next=pb;
95. }
98. void  print_poly(node *poly)      /*输出多项式*/
99. {
100.      int  i;   /*控制是否打印加号*/
102.      for (i=0,poly=poly->next;poly;poly=poly->next,i++)
103.     {
104.          if (i>0)
105.             printf( "+" );
106.         printf( "%gX%d" ,poly->coef,poly->exp);
107.     }
108. }
111. void  main()
112. {
113.     node *poly_a;
114.     node *poly_b;
116.      /*创建多项式A、B*/
117.     printf( "请输入多项式A项的系数和指数/n" );
118.     poly_a=create_poly();
119.     printf( "请输入多项式B项的系数和指数/n" );
120.     poly_b=create_poly();
121.      /*输出多项式A、B*/
122.     printf( "A=" );
123.     print_poly(poly_a);
124.     printf( "/n" );
125.     printf( "B=" );
126.     print_poly(poly_b);
127.     printf( "/n" );
128.      /*A、B相加*/
129.     add_poly(poly_a,poly_b);
130.      /*输出结果*/
131.     printf( "A+B=" );
132.     print_poly(poly_a);
133.     printf( "/n" );
134. }