一元n次多项式的存储和处理
【前言】
这是我写的第一篇博客,如有错误请指正,共同进步。下面的程序代码实现了一元n次多项式的创建,输出,加法,减法,乘法,求函数值,如果没有读者想要实现的功能,可以节约时间离开了。当然如果想要直接借用代码,翻到【最终代码】查阅即可(平时自己在CSDN上查找资料的时候,效率总是很慢,笔者也希望用此提示告诉大家以下代码可以实现什么)。最后笔者作为一个小白,我相信读者稍花时间,一定能明白下面的代码。
【题目详细】
一元n次多项式定义如下:
其中Ai为实数,i为不小于0的整数。定义一元n次多项式操作包括:
(1)输入各个系数Ai和指数i,创建一个多项式;
(2) 输出多项式:将一元n次多项式输出为:F(X)=AnX^n+...+A1X+A0
(3) 两个多项式相加:输入两个多项式,求出并输出两个多项式的和;
(4) 两个多项式相减:输入两个多项式,求出并输出两个多项式的差;
(5) 两个多项式相乘:输入两个多项式,求出并输出两个多项式的积;
(6) 求函数值:输入X,求出并输出一元n次多项式的值。
输入格式:
有两个一元n次多项式,分别为:
f(X)=3X^2+ X+1
g(X)=−2X^2−X-1
其中系数为实数,指数取不小于0的整数,求x=10时f(x)的值。则输入分为3行,第1行为第一个一元n次多项式,第1个一元n次多项式按照第1项系数,指数 第2项系数,指数 .... 的格式输入,系数和指数以“,”分割,各项的系数和指数之间以空格分割,输入一元n次多项式不要求由高次项到低次项排列,最后以 0,0(即系数=0,指数=0)表示结束,输入多项式可以为空。第2行为第二个一元n次多项式,输入格式与第一个一元n次多项式相同。第3行为X的值。对上面的两个一元n次多项式:
输入样例:
3,2 1,1 1,0 0,0
-2,2 -1,1 -1,0 0,0
10
输出格式:
输出分为以下几行:第1行输出第1个一元n次多项式,第2行输出第2个一元n次多项式,第3行输出两个一元n次多项式的和,第4行输出两个一元n次多项式的差,第5行输出两个一元n次多项式的乘积,第6行输出第1个一元n次多项式带入X后求出的值。输出要求一元n次多项式的高次项在前,低次项在后,实数保留小数点后面1位数,一元多项式为空时输出为空,即:若g(x)输入为空,输为g(x)=。对上面2个一元n次多项式的输出为:
输出样例:
f(x)=3.0X^2+X+1.0
g(x)=-2.0X^2-X-1.0
f(x)+g(x)=X^2
f(x)-g(x)=5.0X^2+2.0X+2.0
f(x)*g(x)=-6.0X^4-5.0X^3-6.0X^2-2.0X-1.0
f(10.0)=311.0
【分步代码展示】
开头部分
#include<bits/stdc++.h>
using namespace std;
笔者的建议是直接万能头文件,省时好记,写代码前直接敲上去,管他三七二十一。
------------------------------------------------------我是分割线---------------------------------------------------------
main函数部分
int main()
{
poly pl1[1000];
poly pl2[1000];
int i = 0,j = 0;
double value;
int end_flag1 = 0,end_flag2 = 0;
do
{
i++;
scanf("%lf,%lf",&pl1[i].xishu,&pl1[i].zhishu);
end_flag1 = i;
}
while(pl1[i].xishu != 0 || pl1[i].zhishu != 0);
int a;
while ((a=getchar()) != '\n' && a != EOF);
do
{
j++;
scanf("%lf,%lf",&pl2[j].xishu,&pl2[j].zhishu);
end_flag2 = j;
}
while(pl2[j].xishu != 0.0 || pl2[j].zhishu != 0.0);
int b;
while ((b=getchar()) != '\n' && b != EOF);
cin >> value;
//标准化多项式,指数降序排列
sort(pl1+1,pl1+end_flag1+1,cmp);
sort(pl2+1,pl2+end_flag2+1,cmp);
}
首先定义了两个结构体数组(poly pl1,pl2)用于记录两多项式的系数、指数。结构体创建如下:
typedef struct feature
{
double zhishu;
double xishu;
} poly;
然后定义了两个循环用的万年常客i、j变量。一个双精度浮点数Double类型的用于读入输出函数值部分的自变量x的值。定义了两个记录结构体数组结束位置的变量(定义结构体数组的时候,我大手一挥,干了个1000,说白了这两个变量来收拾我的残局)。
后面便是循环输入两个多项式的系数,指数(遇见0,0停止)。scanf有个特殊的地方,他会把我本来转行输入的回车给读进去,导致输入后存储环节出问题。为了避免,笔者直接派a,b给读走了。
sort函数排序按照cmp规则,按照指数降序排序。cmp函数如下:
int cmp(poly a,poly b)
{
return a.zhishu > b.zhishu;
}
------------------------------------------------------我是分割线---------------------------------------------------------
函数部分
output输出
传入两个参数,记载着一个多项式系数、指数的结构体数组、该数组结束位置。
num:遍历结构体数组所记录的系数,如果为零,num++。如果num==end_flag说明改多项式系数都为零,即该多项式为0。
showpos:输出数值前的正负号。
noshowpos:不输出数值前的正负号。
fixed << setpresicion(x):控制输出x位小数。
这一块函数笔者可能处理的比较复杂,考虑到0,1在系数指数中的存在规则。笔者把第一项与其他项分开讨论,而在每一情况中又分开讨论了系数为0/1,不为0/1,指数为0/1,不为0/1。(系数为0这一项不用输出,系数为1要省去系数,指数为1要省去指数,指数为0要省去x)仔细想一想可以想明白应该。
void output(poly x[],int end_flag)
{
int i,num=0;
for(i=1; i<=end_flag; i++)
{
if(x[i].xishu == 0)
{
num++;
}
}
if(num == end_flag)
{
cout << "0";
}
else
{
if(x[1].xishu != 0.0)
{
if(x[1].xishu == 1.0)
{
if(x[1].zhishu == 1)
{
cout << "X" ;
}
else
{
cout << "X^" << noshowpos << fixed << setprecision(0) << x[1].zhishu;
}
}
else if(x[1].xishu > 0 && x[1].xishu != 1)
{
if(x[1].zhishu == 1)
{
cout << fixed << setprecision(1) << noshowpos << x[1].xishu << "X";
}
else
{
cout << fixed << setprecision(1) << noshowpos << x[1].xishu << "X^" ;
cout << noshowpos << fixed << setprecision(0) << x[1].zhishu;
}
}
else
{
if(x[1].zhishu == 1)
{
cout << fixed << setprecision(1) << showpos << x[1].xishu << "X";
}
else
{
cout << fixed << setprecision(1) << showpos << x[1].xishu << "X^" ;
cout << noshowpos << fixed << setprecision(0) << x[1].zhishu;
}
}
}
for(i = 2; i <=end_flag; i++)
{
if(x[i].xishu != 0.0)
{
if(x[i].xishu == 1.0)
{
if(x[i].zhishu == 1)
{
cout << "+" << "X" ;
}
else if(x[i].zhishu == 0)
{
cout << "+1.0";
}
else
{
cout << "+" << "X^" << noshowpos << x[i].zhishu;
}
}
else if(x[i].xishu == -1.0)
{
if(x[i].zhishu == 1)
{
cout << "-" << "X" ;
}
else if(x[i].zhishu == 0)
{
cout << "-1.0";
}
else
{
cout << "-" << "X^" << noshowpos << x[i].zhishu;
}
}
else
{
if(x[i].zhishu == 1)
{
cout << showpos << fixed << setprecision(1) << x[i].xishu << "X";
}
else if(x[i].zhishu == 0)
{
cout << showpos << fixed << setprecision(1) << x[i].xishu;
}
else
{
cout << showpos << fixed << setprecision(1) << x[i].xishu << "X^" ;
cout << noshowpos << fixed << setprecision(0) << x[i].zhishu;
}
}
}
}
}
}
加法
传入四个参数,两个记载着一个多项式系数、指数的结构体数组、两个该数组结束位置。
end_flag==1:说明多项式1=0,直接输出多项式2。
poly pl_add[1000]:定义一个新的结构体数组用来储存加法运算后的多项式系数、指数。首先把两个多项式接在一起,多项式2接在多项式1的后面。排序。然后合并同类项。
void add(poly x[],poly y[],int end_flag1,int end_flag2)
{
if(end_flag1 == 1)
{
output(y,end_flag2);
}
else if(end_flag2 == 1)
{
output(x,end_flag1);
}
else
{
int i,j,cnt;
cnt = end_flag1 + end_flag2;
poly pl_add[1000];
for(i=1; i<=end_flag1; i++)
{
pl_add[i].xishu = x[i].xishu;
pl_add[i].zhishu = x[i].zhishu;
}
for(j=end_flag1+1; j<=end_flag1+end_flag2; j++)
{
pl_add[j] .xishu = y[j-end_flag1].xishu;
pl_add[j].zhishu = y[j-end_flag1].zhishu;
}
sort(pl_add+1,pl_add+cnt+1,cmp);
for(i=1; i<=end_flag1+end_flag2; i++)
{
for(j=1; j<=end_flag1+end_flag2-1; j++)
{
if(pl_add[i].zhishu == pl_add[i+j].zhishu)
{
pl_add[i].xishu = pl_add[i].xishu + pl_add[i+j].xishu;
pl_add[i+j].xishu = 0;
}
}
}
output(pl_add,cnt);
}
}
减法
传入四个参数,两个记载着一个多项式系数、指数的结构体数组、两个该数组结束位置。
end_flag==1:说明多项式1=0,直接输出-多项式2。
poly pl_sub[1000]:定义一个新的结构体数组用来储存加法运算后的多项式系数、指数。
void sub(poly x[],poly y[],int end_flag1,int end_flag2)
{
if(end_flag1 == 1)
{
for(int k=1; k<=end_flag2; k++)
{
y[k].xishu = -1 * y[k].xishu;
}
output(y,end_flag2);
}
else if(end_flag2 == 1)
{
output(x,end_flag1);
}
else
{
int i,j,k,cnt;
cnt = end_flag1 + end_flag2;
poly pl_sub[1000];
for(i=1,j=1,k=1; k<=cnt;)
{
if(x[i].zhishu == y[j].zhishu)
{
pl_sub[k].zhishu = x[i].zhishu;
pl_sub[k].xishu = x[i].xishu - y[i].xishu;
i++;
j++;
k++;
}
else
{
if(x[i].zhishu > y[j].zhishu)
{
pl_sub[k].zhishu = x[i].zhishu;
pl_sub[k].xishu = x[i].xishu;
i++;
k++;
}
else
{
if(y[i].xishu < 0)
{
pl_sub[k].zhishu = y[i].zhishu;
pl_sub[k].xishu = -1 * y[i].xishu;
j++;
k++;
}
else
{
pl_sub[k].zhishu = y[i].zhishu;
pl_sub[k].xishu = y[i].xishu;
j++;
k++;
}
}
}
}
output(pl_sub,cnt);
}
}
乘法
传入四个参数,两个记载着一个多项式系数、指数的结构体数组、两个该数组结束位置。
poly pl_multi[1000]:定义一个新的结构体数组用来储存加法运算后的多项式系数、指数。
乘法与加法有着很大的相似性,乘法双层循环相乘后得出的数组,排序,然后合并同类项。
void multi(poly x[],poly y[],int end_flag1,int end_flag2)
{
int i,j,cnt=0;
poly pl_multi[1000];
poly pl_multi_[1000];
for(i=1; i<=end_flag1; i++)
{
for(j=1; j<=end_flag2; j++)
{
cnt++;
pl_multi[cnt].xishu = x[i].xishu * y[j].xishu;
pl_multi[cnt].zhishu = x[i].zhishu + y[j].zhishu;
}
}
for(i=1; i<=cnt; i++)
{
for(j=1; j<=cnt; j++)
{
if(pl_multi[i].zhishu == pl_multi[i+j].zhishu)
{
pl_multi[i].xishu = pl_multi[i].xishu + pl_multi[i+j].xishu;
pl_multi[i+j].xishu = 0;
}
}
}
output(pl_multi,cnt);
}
输出函数值
传入三个参数,记载着一个多项式系数、指数的结构体数组、自变量值、该数组结束位置。
很简单看看就懂。
void output_value(poly x[],double value,int endflag)
{
int ans = 0;
double temp;
for(int i = 1; i<endflag-1; i++)
{
temp = value;
for(int j=2; j<=x[i].zhishu; j++)
{
temp = temp * value;
}
ans = ans + x[i].xishu * temp;
}
if(ans ==0)
{
cout << "0";
}
else
{
cout << fixed << setprecision(1) << ans + x[endflag - 1].xishu;
}
}
调用函数
在main函数里调用写的函数,只不过要加入一些补充输出的成分。具体在最终代码部分展示。
------------------------------------------------------我是分割线---------------------------------------------------------
【最终代码展示】
#include<bits/stdc++.h>
using namespace std;
typedef struct feature
{
double zhishu;
double xishu;
} poly;
int cmp(poly a,poly b)
{
return a.zhishu > b.zhishu;
}
void output(poly x[],int end_flag)
{
int i,num=0;
for(i=1; i<=end_flag; i++)
{
if(x[i].xishu == 0)
{
num++;
}
}
if(num == end_flag)
{
cout << "0";
}
else
{
if(x[1].xishu != 0.0)
{
if(x[1].xishu == 1.0)
{
if(x[1].zhishu == 1)
{
cout << "X" ;
}
else
{
cout << "X^" << noshowpos << fixed << setprecision(0) << x[1].zhishu;
}
}
else if(x[1].xishu > 0 && x[1].xishu != 1)
{
if(x[1].zhishu == 1)
{
cout << fixed << setprecision(1) << noshowpos << x[1].xishu << "X";
}
else
{
cout << fixed << setprecision(1) << noshowpos << x[1].xishu << "X^" ;
cout << noshowpos << fixed << setprecision(0) << x[1].zhishu;
}
}
else
{
if(x[1].zhishu == 1)
{
cout << fixed << setprecision(1) << showpos << x[1].xishu << "X";
}
else
{
cout << fixed << setprecision(1) << showpos << x[1].xishu << "X^" ;
cout << noshowpos << fixed << setprecision(0) << x[1].zhishu;
}
}
}
for(i = 2; i <=end_flag; i++)
{
if(x[i].xishu != 0.0)
{
if(x[i].xishu == 1.0)
{
if(x[i].zhishu == 1)
{
cout << "+" << "X" ;
}
else if(x[i].zhishu == 0)
{
cout << "+1.0";
}
else
{
cout << "+" << "X^" << noshowpos << x[i].zhishu;
}
}
else if(x[i].xishu == -1.0)
{
if(x[i].zhishu == 1)
{
cout << "-" << "X" ;
}
else if(x[i].zhishu == 0)
{
cout << "-1.0";
}
else
{
cout << "-" << "X^" << noshowpos << x[i].zhishu;
}
}
else
{
if(x[i].zhishu == 1)
{
cout << showpos << fixed << setprecision(1) << x[i].xishu << "X";
}
else if(x[i].zhishu == 0)
{
cout << showpos << fixed << setprecision(1) << x[i].xishu;
}
else
{
cout << showpos << fixed << setprecision(1) << x[i].xishu << "X^" ;
cout << noshowpos << fixed << setprecision(0) << x[i].zhishu;
}
}
}
}
}
}
void add(poly x[],poly y[],int end_flag1,int end_flag2)
{
if(end_flag1 == 1)
{
output(y,end_flag2);
}
else if(end_flag2 == 1)
{
output(x,end_flag1);
}
else
{
int i,j,cnt;
cnt = end_flag1 + end_flag2;
poly pl_add[1000];
for(i=1; i<=end_flag1; i++)
{
pl_add[i].xishu = x[i].xishu;
pl_add[i].zhishu = x[i].zhishu;
}
for(j=end_flag1+1; j<=end_flag1+end_flag2; j++)
{
pl_add[j] .xishu = y[j-end_flag1].xishu;
pl_add[j].zhishu = y[j-end_flag1].zhishu;
}
sort(pl_add+1,pl_add+cnt+1,cmp);
for(i=1; i<=end_flag1+end_flag2; i++)
{
for(j=1; j<=end_flag1+end_flag2-1; j++)
{
if(pl_add[i].zhishu == pl_add[i+j].zhishu)
{
pl_add[i].xishu = pl_add[i].xishu + pl_add[i+j].xishu;
pl_add[i+j].xishu = 0;
}
}
}
output(pl_add,cnt);
}
}
void sub(poly x[],poly y[],int end_flag1,int end_flag2)
{
if(end_flag1 == 1)
{
for(int k=1; k<=end_flag2; k++)
{
y[k].xishu = -1 * y[k].xishu;
}
output(y,end_flag2);
}
else if(end_flag2 == 1)
{
output(x,end_flag1);
}
else
{
int i,j,k,cnt;
cnt = end_flag1 + end_flag2;
poly pl_sub[1000];
for(i=1,j=1,k=1; k<=cnt;)
{
if(x[i].zhishu == y[j].zhishu)
{
pl_sub[k].zhishu = x[i].zhishu;
pl_sub[k].xishu = x[i].xishu - y[i].xishu;
i++;
j++;
k++;
}
else
{
if(x[i].zhishu > y[j].zhishu)
{
pl_sub[k].zhishu = x[i].zhishu;
pl_sub[k].xishu = x[i].xishu;
i++;
k++;
}
else
{
if(y[i].xishu < 0)
{
pl_sub[k].zhishu = y[i].zhishu;
pl_sub[k].xishu = -1 * y[i].xishu;
j++;
k++;
}
else
{
pl_sub[k].zhishu = y[i].zhishu;
pl_sub[k].xishu = y[i].xishu;
j++;
k++;
}
}
}
}
output(pl_sub,cnt);
}
}
void multi(poly x[],poly y[],int end_flag1,int end_flag2)
{
int i,j,cnt=0;
poly pl_multi[1000];
poly pl_multi_[1000];
for(i=1; i<=end_flag1; i++)
{
for(j=1; j<=end_flag2; j++)
{
cnt++;
pl_multi[cnt].xishu = x[i].xishu * y[j].xishu;
pl_multi[cnt].zhishu = x[i].zhishu + y[j].zhishu;
}
}
for(i=1; i<=cnt; i++)
{
for(j=1; j<=cnt; j++)
{
if(pl_multi[i].zhishu == pl_multi[i+j].zhishu)
{
pl_multi[i].xishu = pl_multi[i].xishu + pl_multi[i+j].xishu;
pl_multi[i+j].xishu = 0;
}
}
}
output(pl_multi,cnt);
}
void output_value(poly x[],double value,int endflag)
{
int ans = 0;
double temp;
for(int i = 1; i<endflag-1; i++)
{
temp = value;
for(int j=2; j<=x[i].zhishu; j++)
{
temp = temp * value;
}
ans = ans + x[i].xishu * temp;
}
if(ans ==0)
{
cout << "0";
}
else
{
cout << fixed << setprecision(1) << ans + x[endflag - 1].xishu;
}
}
int main()
{
poly pl1[1000];
poly pl2[1000];
int i = 0,j = 0;
double value;
int end_flag1 = 0,end_flag2 = 0;
do
{
i++;
scanf("%lf,%lf",&pl1[i].xishu,&pl1[i].zhishu);
end_flag1 = i;
}
while(pl1[i].xishu != 0 || pl1[i].zhishu != 0);
int a;
while ((a=getchar()) != '\n' && a != EOF);
do
{
j++;
scanf("%lf,%lf",&pl2[j].xishu,&pl2[j].zhishu);
end_flag2 = j;
}
while(pl2[j].xishu != 0.0 || pl2[j].zhishu != 0.0);
int b;
while ((b=getchar()) != '\n' && b != EOF);
cin >> value;
//标准化多项式,指数降序排列
sort(pl1+1,pl1+end_flag1+1,cmp);
sort(pl2+1,pl2+end_flag2+1,cmp);
//下面输出第一个多项式
cout << "f(x)=";
if(end_flag1 != 1)
{
output(pl1,end_flag1);
}
cout << endl;
cout << "g(x)=";
if(end_flag2 != 1)
{
output(pl2,end_flag2);
}
cout << endl;
//下面输出两个多项式的加法
cout << "f(x)+g(x)=";
if(end_flag1 != 1 || end_flag2 !=1)
{
add(pl1,pl2,end_flag1,end_flag2);
}
cout << endl;
//下面输出两个多项式减法
cout << "f(x)-g(x)=";
if(end_flag1 != 1 || end_flag2 !=1)
{
sub(pl1,pl2,end_flag1,end_flag2);
}
cout << endl;
//下面输出两个多项式相乘
cout << "f(x)*g(x)=";
if(end_flag1 != 1 && end_flag2 !=1)
{
multi(pl1,pl2,end_flag1,end_flag2);
}
cout << endl;
//下面输出多项式的值
cout << "f(" << noshowpos << fixed << setprecision(1) << value << ")=";
output_value(pl1,value,end_flag1);
return 0;
}
博客代码编写不易,如有错误请指正,如赞同还请点赞关注。
------------------------------------------------------我是分割线---------------------------------------------------------
【效果展示】