多项式操作的c++实现

这篇博客介绍了如何使用C++实现多项式操作,包括结构体定义、链表操作、多项式节点插入、查找、初始化、删除等基本操作,以及多项式相加、相减和相乘的算法实现。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

#include"define.h"

typedef struct
{
 float coef;//系数
 int expn;//指数
}term,ElemType;
typedef struct LNode
{
    ElemType data;
    struct LNode *next;
}*Link;
typedef struct
{
    Link head,tail;
    int len;
}LinkList;

typedef LinkList polynomial;

status MakeNode(Link &p,ElemType e)
{
    //分配由p指向的值为e的节点
    p=(Link)malloc(sizeof(LNode));
    if(p)
    {
        p->data.coef=e.coef;
  p->data.expn=e.expn;
        p->next=NULL;
        return OK;
    }
    else
        return 0;
}

Link GetHead(LinkList L)
{
    return L.head;
}

status SetCurElem(Link &p,ElemType e)
{
    //用e更新p所指节点
 if(!p)
 {
  printf("p is NULL!");
  return ERROR;
 }
    p->data.coef=e.coef;
 p->data.expn=e.expn;
 return OK;
}

Link LocateElem(LinkList L,ElemType e,status(* compare)(ElemType,ElemType))
{
    //返回线性链表中第一个与e满足compare()判定关系的位置;没有则返回NULL
    Link p;
    p=L.head->next;
    while(p!=NULL&&!(*compare)(p->data,e))
        p=p->next;
    if(p!=NULL)
        return p;
    else
        return NULL;
}

status InitList(LinkList &L)
{
    L.head=L.tail=(Link)malloc(sizeof(LNode));
    if(!L.head)
        exit(-2);
    L.len=0;//长度不算头结点
    L.head->next=NULL;
    return OK;
}

status ecqual(ElemType a,ElemType b)
{
 if(a.expn==b.expn)
  return OK;
 return ERROR;
}

void InsFirst(Link h,Link s)
{
    //插入节点(头插)
    s->next=h->next;
    h->next=s;
}

void GetCurElem(Link p,ElemType &a)
{
    //p所直节点的值赋给a
    a.coef=p->data.coef;
 a.expn=p->data.expn;
}

Link NextPos(LinkList L,Link p)
{
    if(!p->next)
  //cout<<"The End!"<<endl;
        return NULL;
    return p->next;
}

status cmp(ElemType a,ElemType b)
{
    if(a.expn>b.expn)
        return 1;
    else if(a.expn==b.expn)
        return 0;
    else
        return -1;
}

status ListEmpty(LinkList L)
{
    //空,返回true
    if(L.head->next==NULL)
        return true;
    else
        return false;
}

status DelFirst(Link h,Link &p)
{
    //删除表中第一个节点,并以p返回
    if(h->next==NULL)
    {
        cout<<"this is empty!"<<endl;
        p=h->next;
        return 0;
    }
    p=h->next;
    h->next=p->next;
    p->next=NULL;
    return OK;
}

void FreeNode(Link &p)
{
    free(p);
}

status ClearList(LinkList &L)
{
    Link p,h;
 h=L.head;
    p=h->next;
    while(p!=NULL)
    {
        h=p->next;
        FreeNode(p);
        p=h;
    }
 
    if(!p)
    {
  L.head->next=NULL;
        L.tail=L.head;
        L.len=0;
        return OK;
    }
    else
        return ERROR;
}

void Append(LinkList &L,Link s)
{
    //将s所指的一串节点链接到L的表尾
    int i=0;
 while(L.tail->next)
 {
  L.tail=L.tail->next;
 }
 L.tail->next=s;
    while(L.tail->next)
    {
        L.tail=L.tail->next;
        i++;
    }
    L.len+=i;
}

void DestroyPolyn(polynomial &p)
{
 //销毁多项式
 if(p.head)
 {
  FreeNode(p.head);
  p.head=p.head->next;
 }
}

void PrintPolyn(polynomial p)
{
    //打印
 Link h;
 h=GetHead(p);
 h=h->next;
 if(h->data.expn==0)
  printf("%f",h->data.coef);
  printf("%fx^%d",h->data.coef,h->data.expn);
 h=h->next;
 while(h)
 {
  if(h->data.expn==0)
   printf("%f",h->data.coef);
  else
  {
   if(h->data.coef>0)
    printf("+");
   printf("%fx^%d",h->data.coef,h->data.expn);
  }
  h=h->next;
 }
    printf("/n");
}

int PolynLength(polynomial p)
{
 //返回多项式的项数
 int i;
 i=0;
 while(p.head)
 {
  i++;
  p.head=p.head->next;
 }
 return i;
}

void CreatePolyn(polynomial &p,int m)
{
    //输入m项的指数和系数,建立表示一元多项式的有序链表p
 Link h,s;
 term e;
 int i;
 InitList(p);
    h=GetHead(p);
    e.coef=0.0;
 e.expn=-1;//设置头结点
 SetCurElem(h,e);
 //cout<<h->data.coef<<h->data.expn;
 cout<<"Please input coef and expn"<<"("<<m<<",Descending)"<<":"<<endl;
 for(i=1;i<=m;i++)
    {
  //cin>>e.coef>>e.expn;
  scanf("%f",&e.coef);
  scanf("%d",&e.expn);
  if(!LocateElem(p,e,ecqual))
   if(MakeNode(s,e))
    InsFirst(h,s);
 }
}

void AddPolyn(polynomial &pa,polynomial &pb,status(* cmp)(ElemType,ElemType))
{
 Link ha,hb,qa,qb;
 ElemType a,b;
 //int T;
 float sum;
 ha=GetHead(pa);
 hb=GetHead(pb);
 qa=NextPos(pa,ha);
 qb=NextPos(pb,hb);
 while(qa&&qb)
 {
        GetCurElem(qa,a);
  GetCurElem(qb,b);
  //T=(*cmp)(a,b);
  switch((*cmp)(a,b))
  {
  case -1://多项式pa中当前节点的值较小
   ha=qa;
   qa=NextPos(pa,qa);
   break;
  case 0:
   sum=a.coef+b.coef;
   if(sum!=0)
   {
    qa->data.coef=sum;//修改系数
    ha=qa;
   }
   else
   {
    DelFirst(ha,qa);
    FreeNode(qa);
   }
   DelFirst(hb,qb);
   FreeNode(qb);
   qb=NextPos(pb,hb);
   qa=NextPos(pa,ha);
   break;
  case 1:
   DelFirst(hb,qb);
   InsFirst(ha,qb);
   qb=NextPos(pb,hb);
   ha=NextPos(pa,ha);
   break;
  }
 }
 if(!ListEmpty(pb))
  Append(pa,qb);
 FreeNode(hb);
}

void SubtractPolyn(polynomial &pa,polynomial &pb)
{
 //多项式相减,并销毁pb
 Link ha,hb,qa,qb;
 ElemType a,b;
 float sub;
 ha=GetHead(pa);
 hb=GetHead(pb);
 qa=NextPos(pa,ha);
 qb=NextPos(pb,hb);
    while(qa&&qb)
 {
  GetCurElem(qa,a);
  GetCurElem(qb,b);
  switch((*cmp)(a,b))
  {
  case -1:
   ha=qa;
   qa=NextPos(pa,qa);
   break;
  case 0:
   sub=a.coef-b.coef;
   if(sub!=0.0)
   {
    qa->data.coef=sub;
    ha=qa;
   }
   else
   {
                DelFirst(ha,qa);
    FreeNode(qa);
   }
   DelFirst(hb,qb);
   FreeNode(qb);
   qb=NextPos(pb,hb);
   qa=NextPos(pa,ha);
   break;
   
  case 1:
   DelFirst(hb,qb);
   qb->data.coef*=-1;
   InsFirst(ha,qb);
   qb=NextPos(pb,hb);
   ha=NextPos(pa,ha);
   break;
  }
 }
 if(!ListEmpty(pb))
  Append(pa,qb);
 FreeNode(hb);
}

/*void  contrary(polynomial &q)
{
 //逆序
 Link p,c,n;
    p=NULL;//前一个节点
    c=q.head;//当前节点
    n=q.head->next;//下一个节点
    while(n){//如果没有到链表尾
        c->next=p;//逆序当前节点
        p=c;//节点前进
        c=n;//节点前进
        n=n->next;//节点前进
    }
    c->next=p;
    q.head=c;
}*/

void MultiplyPolyn(polynomial &pa,polynomial &pb)
{
 //多项式相乘
    polynomial sum,T;//sum保存相加结果,T保存每步相乘结果
 Link hs,hT,ha,hb,qa,qb,s,t;
 term e;
 // ElemType a,b;
 // int m,n
 /*int la,lb;
 la=PolynLength(pa);
 lb=PolynLength(pb);*/
    InitList(sum);
 InitList(T);
 cout<<1<<endl;
 e.coef=0.0;
 e.expn=-1;//设置头结点
 hT=GetHead(T);
 hs=GetHead(sum);
 SetCurElem(hs,e);
 SetCurElem(hT,e);
    ha=GetHead(pa);
 if(ha->next==NULL)
  cout<<"re"<<endl;
 qa=NextPos(pa,ha);
 cout<<2<<endl;
 hb=GetHead(pb);
    while(qa)
 {
  qb=NextPos(pb,hb);
  cout<<3<<endl;
  while(qb)
  {
   
   e.coef=qa->data.coef*qb->data.coef;
   e.expn=qa->data.expn+qb->data.expn;
   if(MakeNode(s,e))
   {
    if(hT->next==NULL)
    {
     hT->next=s;
     t=s;
                }
    else
    {
     t->next=s;
     s->next=NULL;
     t=s;
    }
   }
   qb=qb->next;
  }
  qa=qa->next;
        AddPolyn(sum,T,cmp);
  PrintPolyn(sum);
        InitList(T);
  hT=GetHead(T);
  cout<<4<<endl;
 }
    DestroyPolyn(pb);
 hs=GetHead(sum);
 hs=NextPos(sum,hs);
 ClearList(pa);
    Append(pa,hs);
}

void main()
{
 polynomial a,b;
 CreatePolyn(a,2);
 CreatePolyn(b,2);
 PrintPolyn(a);
 PrintPolyn(b);
 //AddPolyn(a,b,cmp);
 //MultiplyPolyn(a,b);
 SubtractPolyn(a,b);
 PrintPolyn(a);
}

#ifndef FUNCTION_H_ #define FUNCTION_H_ #include #include #include "polyfit.h" #include using namespace std; dxs::dxs() { ifstream fin("多项式拟合.txt"); fin>>n; x=new float[n]; y=new float[n]; for(int i=0;i>x[i]; } for(i=0;i>y[i]; } cout<>nn; m=nn+1; u=new float*[m]; for(i=0;i<m;i++) { u[i]=new float[m+1]; }//创建m行,m+1列数组 } void dxs::dfine() { for(int i=0;i<m;i++) { for(int j=0;j<m+1;j++) { u[i][j]=0; } } for(i=0;i<m;i++) { for(int j=0;j<m;j++) { for(int k=0;k<n;k++) { u[i][j]=u[i][j]+pow(x[k],j+i); } } } for(i=0;i<m;i++) { for(int k=0;k<n;k++) { u[i][m]=u[i][m]+pow(x[k],i)*y[k]; } } } void dxs::show() { for(int i=0;i<m;i++) { for(int j=0;j<m+1;j++) { cout<<u[i][j]<<" ";//<<endl; } cout<<endl; } ////显示具有m行m+1列u数组的各元素值 } void dxs::select_main(int k,float **p,int m) { double d; d=*(*(p+k)+k); //cout<<d; int l=k; int i=k+1; for(;i fabs(d)) { d=*(*(p+i)+k); l=i; } else continue; } if(d==0) cout<<"错误"; else { if(k!=l) { for(int j=k;j<m+1;j++) { double t; t=*(*(p+l)+j); *(*(p+l)+j)=*(*(p+k)+j); *(*(p+k)+j)=t; } } } } void dxs::gaosi() { for(int k=0;k<m;k++) { select_main(k,u,m);//调用列主元函数 for(int i=1+k;i<m;i++) { // *(*(p+i)+k)=(float) *(*(p+i)+k) / *(*(p+k)+k); u[i][k]=(float) u[i][k] / u[k][k]; } for(i=k+1;i<m;i++) { for(int j=k+1;j=0;i--) { float a=0; for(int j=i+1;j<m;j++) { //a=a + (*(*(p+i)+j) * *(*(p+j)+m)); a=a+u[i][j] * u[j][m]; } //*(*(p+i)+n-1)= (*(*(p+i)+n-1) - a) / *(*(p+i)+i); u[i][m]= (u[i][m] -a) / u[i][i]; } cout<<"方程组的解为:"<<endl; for(i=0;i<m;i++) { cout<<"a"<<i+1<<"="; cout<<u[i][m]<<endl; // l[i]=*(*(p+i)+n-1); } cout<<"y="<<u[0][m]; for(i=1;i<m;i++) { cout<<showpos<<u[i][m]<<"x"; if(i!=1)cout<<"^"<<noshowpos<<i; } cout<<endl; } dxs::~dxs() { delete[]x,y; delete []*u; } #endif
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值