C++类实现多项式——加减乘运算

一、创建类Polynml
        在多项式类Polynml的内部,分别由函数变量,构造函数,析构函数和成员函数组成;另外还有重载运算符:①重载加法(+),②重载减法(-),③重载乘法(*)。

        注意多项式的每一项是用二维数组来构成,每个数组的第一维表示系数,第二维表示指数。(例如 2x^4 ;2为系数,4为指数)


二、代码实现

#include<iostream>

using namespace std;

class Polynml{ //定义多项式类
    int n;
    double (*p)[2];
public:

    Polynml(double poly[][2],int num){ //构造函数,带两个参数
        n=num;
        p=new double[n][2];
        for(int i=0;i<n;i++){
            p[i][0]=poly[i][0];
            p[i][1]=poly[i][1];
        }
    }
    Polynml(int num){ //构造函数,带一个参数
        n=num;
        p=new double[num][2];
        for(int i=0;i<n;i++){
            p[i][0]=0;
            p[i][1]=0;
        }
    }

    Polynml(){} //析构函数

    friend Polynml operator+(Polynml &c1,Polynml &c2){ //成员函数,实现+运算
        int i,j,tn;
        Polynml t(c2.n+c1.n);
        tn=c1.n;
        for(i=0;i<c1.n;i++){
            t.p[i][0]=c1.p[i][0];
            t.p[i][1]=c1.p[i][1];
        }
        for(i=0;i<c2.n;i++){
            for(j=0;j<t.n;j++){
                if(c2.p[i][1]==t.p[j][1]){
                    t.p[j][0]+=c2.p[i][0];
                    break;
                }
                else if(j>=tn){
                    t.p[j][0]=c2.p[i][0];
                    t.p[j][1]=c2.p[i][1];
                    tn++;
                    break;
                }
                else if(c2.p[i][1]<t.p[j][1]){
                    for(int k=t.n-1;k>j;k--){
                        t.p[k][0]=t.p[k-1][0];
                        t.p[k][1]=t.p[k-1][1];
                    }
                    t.p[j][0]=c2.p[i][0];
                    t.p[j][1]=c2.p[i][1];
                    tn++;
                    break;
                }
                else if(c2.p[i][1]>t.p[j][1])
                    continue;
            }
        }
            t.n=tn;
            return t;
    }

    friend Polynml operator-(Polynml &c1,Polynml &c2){ //成员函数,实现-运算
        Polynml t(c2.n);
        for(int i=0;i<t.n;i++){
            t.p[i][0]=-c2.p[i][0];
            t.p[i][1]=c2.p[i][1];
        }
        return c1+t;
    }

    friend Polynml operator*(Polynml &c1,Polynml &c2){ //成员函数,实现*运算
        Polynml t(c1.n*c2.n);
        int i,j,m=0;
        for(i=0;i<c1.n;i++){
            for(j=0;j<c2.n;j++){
                t.p[m][0]=c1.p[i][0]*c2.p[j][0];
                t.p[m][1]=c1.p[i][1]+c2.p[j][1];
                m++;
            }
        }
        return t;
    }

    void display(){
        for(int i=0;i<n;i++)
            cout<<p[i][0]<<" "<<p[i][1]<<endl;
            cout<<endl;
    }

    ~Polynml(){
        delete []p;
    }
};

int main()
{
    double a[4][2]={{1,0},{2,2},{3,4},{4,8}},
            b[4][2]={{5,1},{8,4},{9,8},{-1,10}};
    Polynml A(a,4);
    Polynml B(b,4);
    Polynml C(4);
    C=A+B;
    C.display();
    C=A-B;
    C.display();
    C=A*B;
    C.display();

 return 0;

}

三、总结

        在实现 重载乘法 时仍有缺点,并没有合并同类项,故而导致多项式各自独立。

  • 11
    点赞
  • 89
    收藏
    觉得还不错? 一键收藏
  • 9
    评论
#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

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 9
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值