这仍然是一道关于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
思路
自从输入重新使用scanf就老是忘记写&然后浪费大量时间找半天bug,还是要多加练习
多项式除法A/B
求商A最高项除以B最高项所得系数若>=0.05则为一项商q[i],
多项式A=A-B*q[i],(A的最高指数每次循环–)然后继续A/B直到A的最高指数小于B的最高指数
此时A为余数
代码
#include <iostream>
#include <map>
#include <cmath>
#include <cstdio>
using namespace std;
typedef struct stu{
int e;//指数
double c;//系数
}mm;
int main(){
int n,m,a,ma;
cin>>n;
double b;
map<int,double> xx;//不能设置成int数组,因为不知道会有多大的指数
for(int i=0;i<n;i++){
cin>>a>>b;
xx[a]=b;
if(i==0) ma=a;//ma被除数最大指数
}
cin>>m;
mm xx1[m],q[m];//xx1[m]除数,q[m]是商
for(int i=0;i<m;i++){
cin>>xx1[i].e>>xx1[i].c;//xx1[0]处,除数最大指数
}
int p=0;
while(ma>=xx1[0].e){//被除数指数大于除数指数
int dif=ma-xx1[0].e;//指数比
double xs=xx[ma]/xx1[0].c;//系数比
if(fabs(xs)>=0.05){//系数不可舍去
q[p].e=dif;
q[p++].c=xs;
for(int i=0;i<m;i++){
xx[xx1[i].e+dif]-=xx1[i].c*xs;
}
}else ma--;//系数不可舍去
while(ma>=xx1[0].e&&fabs(xx[ma])<0.05) ma--;
}
cout<<p;
if(!p) printf(" 0 0.0\n");
else{
for(int i=0;i<p;i++){
printf(" %d %.1f",q[i].e,q[i].c);
}
printf("\n");
}
p=0;
while(ma>=0){
if(fabs(xx[ma])>=0.05){
q[p].e=ma;
q[p++].c=xx[ma];
}
ma--;
}
cout<<p;
if(!p) printf(" 0 0.0\n");
for(int i=0;i<p;i++){
printf(" %d %.1f",q[i].e,q[i].c);
}
return 0;
}