一元多項式的表示及加減乘除運算

比如：怎樣實現用線性鏈表表示多項式的加法運算？

依據一元多項式相加的運算規則：對於兩個一元多項式中全部指數同樣的項。相應系數相加，若其和不為零，則構成“和多項式”中的一項。對於兩個一元多項式中全部指數不同樣的項，則分別復抄到“和多項式”中去。

#include <stdio.h>  
#include <stdlib.h>  
#include <malloc.h>  
  
typedef struct polyn  
{  
    float coef;  
    int expn;  
    struct polyn* next;  
}PolyNode,*PLinkList;  
  
PLinkList CreatePolyn();//創建一元多項式，使一元多項式呈指數遞減  
void OutPut(PLinkList head);//輸出一元多項式  
PLinkList Addition(PLinkList L1,PLinkList L2);//多項式的加法  
PLinkList Subtraction(PLinkList L1,PLinkList L2);//多項式的減法  
PLinkList Reverse(PLinkList head);//將生成的鏈表逆置。使一元多項式呈指數遞增形式  
PLinkList MultiplyPolyn(PLinkList L1,PLinkList L2);//多項式的乘法 





#include "test.h"  
  
  
PLinkList CreatePolyn()//創建一元多項式。使一元多項式呈指數遞減  
{  
    PolyNode *p,*q,*s;  
    PolyNode *head = NULL;  
    int expn2;  
    float coef2;  
    head = (PLinkList)malloc(sizeof(PolyNode));//動態生成頭結點  
    if(!head)  
    {  
        return NULL;  
    }  
    head->coef = 0.0;//初始化  
    head->expn = 0;  
    head->next = NULL;  
    do  
    {  
        printf("輸入系數coef(系數和指數都為0結束)");  
        scanf("%f",&coef2);  
        printf("輸入指數數exp(系數和指數都為0結束)");  
        scanf("%d",&expn2);  
        if((long)coef2 == 0 && expn2 == 0)  
        {  
            break;  
        }  
        s = (PLinkList)malloc(sizeof(PolyNode));  
        if(!s)  
        {  
            return NULL;  
        }  
        s->expn = expn2;  
        s->coef = coef2;  
        q = head->next ;  
        p = head;  
        while(q && expn2 < q->expn)  
        {  
            p = q;  
            q = q->next ;  
        }  
        if(q == NULL || expn2 > q->expn)  
        {  
            p->next = s;  
            s->next = q;  
        }  
        else  
        {  
            q->coef += coef2;  
        }  
    }while(1);  
    return head;  
}  
  
void OutPut(PLinkList head)//輸出一元多項式  
{  
    PolyNode *p = head->next ;  
    while(p)  
    {  
        printf("%1.1f",p->coef);  
        if(p->expn)  
        {  
            printf("*x^%d",p->expn);  
        }  
        if(p->next && p->next->coef > 0)  
        {  
            printf("+");  
        }  
        p = p->next ;  
    }  
}  
  
PolyNode *Addition(PLinkList L1,PLinkList L2)//多項式的加法  
{  
    PolyNode *pa,*pb,*pc,*u,*head;  
    head = (PLinkList)malloc(sizeof(PolyNode));  
    if(!head)  
    {  
        return NULL;  
    }  
    head->coef = 0.0;  
    head->expn = 0;  
    head->next = NULL;  
    pc = head;  
    L2 = Reverse(L2);  
    pa = L1->next ;  
    pb = L2->next ;  
    while(pa != NULL && pb != NULL)  
    {  
        if(pa->expn == pb->expn)  
        {  
            u = (PLinkList)malloc(sizeof(PolyNode));  
            if(!u)  
            {  
                return NULL;  
            }  
            u->coef = pa->coef + pb->coef ;  
            u->expn = pa->expn ;  
            pa = pa->next ;  
            pb = pb->next ;  
            u->next = pc->next ;  
            pc->next = u;  
            pc = u;  
        }  
        else if(pa->expn > pb->expn)  
        {  
            u = (PLinkList)malloc(sizeof(PolyNode));  
            if(!u)  
            {  
                return NULL;  
            }  
            u->coef = pa->coef ;  
            u->expn = pa->expn ;  
            pa = pa->next ;  
            u->next = pc->next ;  
            pc->next = u;  
            pc = u;  
        }  
        else  
        {  
            u = (PLinkList)malloc(sizeof(PolyNode));  
            if(!u)  
            {  
                return NULL;  
            }  
            u->coef = pb->coef ;  
            u->expn = pb->expn ;  
            pb = pb->next ;  
            u->next = pc->next ;  
            pc->next = u;  
            pc = u;  
        }  
    }  
    L2 = Reverse(L2);  
    return head;  
}  
  
PolyNode *Subtraction(PLinkList L1,PLinkList L2)//多項式的減法  
{  
    PolyNode *pa,*pb,*pc,*u,*head;  
    head = (PLinkList)malloc(sizeof(PolyNode));  
    if(!head)  
    {  
        return NULL;  
    }  
    head->coef = 0.0;  
    head->expn = 0;  
    head->next = NULL;  
    pc = head;  
    pa = L1->next ;  
    pb = L2->next ;  
    while(pa != NULL && pb != NULL)  
    {  
        if(pa->expn == pb->expn)  
        {  
            u = (PLinkList)malloc(sizeof(PolyNode));  
            if(!u)  
            {  
                return NULL;  
            }  
            u->coef = pa->coef - pb->coef ;  
            u->expn = pa->expn ;  
            pa = pa->next ;  
            pb = pb->next ;  
            u->next = pc->next ;  
            pc->next = u;  
            pc = u;  
        }  
        else if(pa->expn > pb->expn)  
        {  
            u = (PLinkList)malloc(sizeof(PolyNode));  
            if(!u)  
            {  
                return NULL;  
            }  
            u->coef = pa->coef ;  
            u->expn = pa->expn ;  
            pa = pa->next ;  
            u->next = pc->next ;  
            pc->next = u;  
            pc = u;  
        }  
        else  
        {  
            u = (PLinkList)malloc(sizeof(PolyNode));  
            if(!u)  
            {  
                return NULL;  
            }  
            u->coef = pb->coef ;  
            u->expn = pb->expn ;  
            pb = pb->next ;  
            u->next = pc->next ;  
            pc->next = u;  
            pc = u;  
        }  
    }  
    return head;  
}  
  
PolyNode *Reverse(PLinkList head)//將生成的鏈表逆置，使一元多項式呈指數遞增形式  
{  
    PolyNode *q,*r,*p = NULL;  
    q = head->next ;  
    while(q)  
    {  
        r = q->next ;  
        q->next = p;  
        p = q;  
        q = r;  
    }  
    head->next = p;  
    return head;  
}  
  
PolyNode *MultiplyPolyn(PLinkList L1,PLinkList L2)//多項式的乘法  
{  
    PolyNode *pa,*pb,*pc,*u,*head;  
    int k,maxExp;  
    float coef;  
    head = (PLinkList)malloc(sizeof(PolyNode));  
    if(!head)  
    {  
        return NULL;  
    }  
    head->coef = 0.0;  
    head->expn = 0;  
    head->next = NULL;  
    if(L1->next != NULL && L2->next != NULL)  
    {  
        maxExp = L1->next->expn +L2->next->expn ;  
    }  
    else  
    {  
        return head;  
    }  
    pc = head;  
    L2 = Reverse(L2);  
    for(k = maxExp;k >= 0;k--)  
    {  
        pa = L1->next ;  
        while(pa != NULL && pa->expn > k)  
        {  
            pa = pa->next ;  
        }  
        pb = L2->next ;  
        while(pb != NULL && pa != NULL && pa->expn+pb->expn < k)  
        {  
            pb= pb->next ;  
        }  
        coef = 0.0;  
        while(pa != NULL && pb != NULL)  
        {  
            if(pa->expn +pb->expn == k)  
            {  
                coef += pa->coef *pb->coef ;  
                pa = pa->next ;  
                pb = pb->next ;  
            }  
            else if(pa->expn +pb->expn > k)  
            {  
                pa = pa->next ;  
            }  
            else  
            {  
                pb = pb->next ;  
            }  
        }  
        if(coef != 0.0)  
        {  
            u = (PLinkList)malloc(sizeof(PolyNode));  
            u->coef = coef;  
            u->expn = k;  
            u->next = pc->next ;  
            pc->next = u;  
            pc = u;  
        }  
    }  
    L2 = Reverse(L2);  
    return head;  
}  





#include "test.h"  
  
int main(void)  
{  
    PLinkList A,B,C,D,E;  
    A = CreatePolyn();  
    printf("A(x) =");  
    OutPut(A);  
    printf("\n");  
    B = CreatePolyn();  
    printf("B(x) =");  
    OutPut(B);  
    printf("\n");  
    C = MultiplyPolyn(A,B);  
    printf("C(x) = A(x)*B(x) =");  
    OutPut(C);  
    printf("\n");  
    D = Addition(A,B);  
    printf("D(x) = A(x)+B(x) =");  
    OutPut(D);  
    printf("\n");  
    E = Subtraction(A,B);  
    printf("E(x) = A(x)-B(x) =");  
    OutPut(E);  
    printf("\n");  
    return 0;  
}