L2-018 多项式A除以B (25 分)

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++;
        }
      }
    }
}

猜你喜欢

转载自blog.csdn.net/q852092907/article/details/83625642