问题描述
实现一元稀疏多项式的表示及基本操作(建立、销毁、输出、加法、减法等操作)
算法设计
实现要定义的一元多项式应采用链式存储结构。根据一元多项式相加的运算法则,对于两个多项式中所有指数相同的项,对应系数相加,若其和不为零,则构成新多项式的一项;对于两个多项式中所有指数不同的项,分别复制到新多项式中。新多项式不必另外生成,而是在原来的两个多项式中摘取结点即可。
采用顺序存储结构存储多项式系数A,使多项式的某些运算变得更简洁。但在实际使用时,多项式的阶数可能很高,不同的多项式阶数可能相差很大,这使顺序存储结构的最大长度难以确定。而且,最高次幂与最低次幂项之间缺项很多时,采用顺序存储结构显然十分浪费存储空间,因此,一般情况下多采用链式存储结构来存储高阶多项式。
在用线性链表来存储一个多项式时,多项式中的每个非零项系数对应一一个结点,结点由数据元素项和指针组成。数据元素项中包含系数和指数值,设计中先定义定义数据元素中的数据,其中有指数、系数级指针next等。并要先构造一元多项式。在定义输出、删除、合并等模块。假设指针qa和qb分别指向多项式A和B中当前进行操作的某个结点,比较这个结点的指数项,可能有三种情况:
①指针qa->exp<qb->exp,则摘取qa指针所指结点插入到和多项式中;
②指针qa->exp<qb->exp,则摘取qb指针所指结点插入到和多项式中;
③指针qa->exp=qb->exp,则将系数相加,若和数不为零,则修改qa->coef的值,同时释放qb所指结点,反之,从多项式A的链表中删除相应的结点,并释放指针Pa
算法实现
#define _CRT_SECURE_NO_WARNINGS
#include<stdio.h>
#include<stdlib.h>
/**********************一元多项式储存结构表示******************************/
typedef struct Polynomial {
int coefficient; //系数
int index; //指数
struct Polynomial *next;
}Plm;
#define SIZE sizeof(Plm)
/**********************模块函数******************************/
Plm *CreatPolyn(int m) //输入m项的系数和指数,建立一元多项式
{
Plm *p,*p1, *p2;
p1 = (Plm *)malloc(SIZE);
scanf("%d %d", &p1->coefficient, &p1->index);
m--;
p = p1;
while (m)
{
p2 = (Plm*)malloc(SIZE);
p1->next = p2;
scanf("%d %d", &p2->coefficient, &p2->index);
p1 = p2;
m--;
}
p1->next = NULL;
return p;
}
void DestroyPolyn(Plm *p) //一元多项式p存在,销毁p
{
Plm *q;
if (p == NULL)
printf("error");
else
while (p != NULL)
{
q = p;
p = p->next;
free(q);
}
}
void PrintPolyn(Plm *p) //输出一元多项式
{
while (p!=NULL)
{
printf("%dX^%d ", p->coefficient, p->index);
if (p->next != NULL && p->next->coefficient>=0) printf("+");
p = p->next;
}
printf("\n");
}
int PolynLength(Plm *p) //返回一元多项式p的项数
{
int i = 0;
while (p != NULL)
{
i++;
p = p->next;
}
return i;
}
Plm *AddPolyn(Plm *pa,Plm *pb) //完成多项式相加并销毁pb
{
Plm *p,*q, *q1;
q = (Plm*)malloc(SIZE);
q1 = q;
while (pa&&pb != NULL)
{
if (pa->index == pb->index)
{
q->next = pa;
q = q->next;
pa->coefficient = pa->coefficient + pb->coefficient;
p = pb;
pb = pb->next;
free(p);
pa = pa->next;
}
else if (pa->index < pb->index)
{
q->next = pa;
q = q->next;
pa = pa->next;
}
else
{
q->next = pb;
q = q->next;
pb = pb->next;
}
}
if (pa != NULL) q->next = pa;
else q->next = pb;
q1 = q1->next;
return q1;
}
Plm *SubtractPolyn(Plm *pa,Plm *pb) //完成多项式相减,pa = pa - pb;并销毁pb
{
Plm *p, *q, *q1;
q = (Plm*)malloc(SIZE);
q1 = q;
while (pa&&pb != NULL)
{
if (pa->index == pb->index)
{
q->next = pa;
q = q->next;
pa->coefficient = pa->coefficient - pb->coefficient;
p = pb;
pb = pb->next;
free(p);
pa = pa->next;
}
else if (pa->index < pb->index)
{
q->next = pa;
q = q->next;
pa = pa->next;
}
else
{
pb->coefficient = -pb->coefficient;
q->next = pb;
q = q->next;
pb = pb->next;
}
}
if (pa != NULL) q->next = pa;
else
{
q->next = pb;
while (pb != NULL)
{
pb->coefficient = -pb->coefficient;
pb = pb->next;
}
}
q1 = q1->next;
return q1;
}
int main()
{
Plm *head_A, *head_B;
int a = 0, b = 0;
//加
printf("输入A的项数:");
scanf("%d", &a);
printf("输入多项式A:");
head_A = CreatPolyn(a);//如何通过creatpolyn改变head的值?把head作为参数传进去而不是让函数返回头指针
printf("输入B的项数:");
scanf("%d", &b);
printf("输入多项式B:");
head_B = CreatPolyn(b);
head_A = AddPolyn(head_A, head_B);
PrintPolyn(head_A);
DestroyPolyn(head_A);
//减
printf("输入A的项数:");
scanf("%d", &a);
printf("输入多项式A:");
head_A = CreatPolyn(a);//如何通过creatpolyn改变head的值?把head作为参数传进去而不是让函数返回头指针
printf("输入B的项数:");
scanf("%d", &b);
printf("输入多项式B:");
head_B = CreatPolyn(b);
head_A = SubtractPolyn(head_A,head_B);
PrintPolyn(head_A);
DestroyPolyn(head_A);
return 0;
}