L2-018 多项式A除以B (25 分)
这仍然是一道关于A/B的题,只不过A和B都换成了多项式。你需要计算两个多项式相除的商Q和余R,其中R的阶数必须小于B的阶数。
输入格式:
输入分两行,每行给出一个非零多项式,先给出A,再给出B。每行的格式如下:
N e[1] c[1] ... e[N] c[N]
其中N
是该多项式非零项的个数,e[i]
是第i
个非零项的指数,c[i]
是第i
个非零项的系数。各项按照指数递减的顺序给出,保证所有指数是各不相同的非负整数,所有系数是非零整数,所有整数在整型范围内。
输出格式:
分两行先后输出商和余,输出格式与输入格式相同,输出的系数保留小数点后1位。同行数字间以1个空格分隔,行首尾不得有多余空格。注意:零多项式是一个特殊多项式,对应输出为0 0 0.0
。但非零多项式不能输出零系数(包括舍入后为0.0)的项。在样例中,余多项式其实有常数项-1/27
,但因其舍入后为0.0,故不输出。
输入样例:
4 4 1 2 -3 1 -1 0 -1
3 2 3 1 -2 0 1
输出样例:
3 2 0.3 1 0.2 0 -1.0
1 1 -3.1
//本题难点在于多项式除法的运算规则,运算法则详情百度
#include<bits/stdc++.h>
using namespace std;
double a[10005],b[10005];//每个数组保存一个多项式
int main (){
map<int,double> op;//记录结果
int n;
cin>>n;
int m;
int maxn1=0,maxn2=0;//maxn1保存多项式A的最大的指数,maxn2保存多项式B的最大的指数。
double l;
for(int i=0;i<n;++i){
scanf("%d %lf",&m,&l);
a[m]=l;
if(m>maxn1) maxn1=m;
}
int n1;
cin>>n1;
for(int i=0;i<n1;++i){
scanf("%d %lf",&m,&l);
b[m]=l;
if(m>maxn2) maxn2=m;
}
for(int i=maxn1;i>=0;--i){//根据运算法则计算结果
if(i<maxn2) break;//被除项最大指数小于除项最大指数结束运算
if(a[i]){
op[i-maxn2]=a[i]/b[maxn2];
for(int j=maxn2;j>=0;--j){
a[i-maxn2+j]-=op[i-maxn2]*b[j];//得到新的被除式
}
}
}
//计算过后,op中保存的是商,数组a中保存的是余
int cnt=0;
for(auto i=op.begin();i!=op.end();++i){
if(fabs(i->second)<0.05) op.erase(i);//小于0.05的数被当作0来看,所以要舍去。
}
printf("%d ",op.size());
if(op.size()==0) printf("0 0.0");
for(auto i=op.rbegin();i!=op.rend();++i){
if(cnt) printf(" ");
printf("%d %.1lf",i->first,i->second);
cnt++;
}
printf("\n");
int flag=0;
for(int i=maxn2-1;i>=0;i--){
if(fabs(a[i])>=0.05){
flag++;
}
}
cnt=0;
printf("%d ",flag);
if(!flag) printf("0 0.0");
else{
for(int i=maxn2-1;i>=0;--i){
if(fabs(a[i])>=0.05){
if(cnt) printf(" ");
printf("%d %.1lf",i,a[i]);
cnt++;
}
}
}
}