多项式A除以B
题目要求:
这仍然是一道关于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
这道题对我来说,还是挺有难度的。。。
对于多项式除法貌似还是第一次做。
被除数/除数 = 商……余
(图片来源:https://www.bbsmax.com/A/MyJxevQM5n/)
仔细一想的话,其实还可以,不算很难理解。
首先,每一步运算它的被除数都是在变的。被除数 = (被除数-(被除数最高项/除数最高项)*除数),直到被除数最高项的指数小于除数。剩下的被除数就是商。这样一想还是可以的,但实现起来,真是。。。依旧是看别人的优秀代码。。
代码如下:
#include<iostream>
#include<cmath>
using namespace std;
int main(){
int a,b;int ma,mb;int zhi,xi;double c[100000];
cin>>a;
double* data1 = new double[1000000]();
for(int i=0;i<a;i++){//数组的索引为多项式的指数,数组的内容为数组的系数
cin>>zhi;
if(i == 0)ma = zhi;
cin>>xi;
data1[zhi] = xi;
}
cin>>b;
double* data2 = new double[1000000]();
for(int i=0;i<b;i++){
cin>>zhi;
if(i == 0)mb = zhi;
cin>>xi;
data2[zhi] = xi;
}
//这一步比较关键,这个指数容易搞混,但理解之后就清楚了
//data1[i]/data[mb]就相当于被除数(指数大于除数最高指数)的某一项的系数除以除数最高项的系数,然后的商:指数为i-mb,系数为data[i]/data[mb]
for(int i=ma;i>=mb;i--){
c[i-mb] = data1[i]/data2[mb];//i-mb为商的指数
for(int j=mb;j>=0;j--){
data1[i+j-mb] -= data2[j]*c[i-mb];//除数每一项都乘以商
//这里为什么是i+j-mb?ennnn我只能说(i+j-mb) = (j)+(i-mb)
}
}
//下面就好理解了
int num1 = 0,num2 = 0;
//商的个数
for(int i=ma-mb;i>=0;i--){
if(abs(c[i]) >= 0.1)num1++;
}
//余数的个数 余数的最高指数比除数的最高指数小
for(int i=mb-1;i>=0;i--){
if(abs(data1[i]) >= 0.1)num2++;
}
if(num1 == 0){
printf("0 0 0.0\n");
}else{
printf("%d",num1);
for(int i=ma-mb;i>=0;i--){
if(abs(c[i]) >=0.1)printf(" %d %.1lf",i,c[i]);
}
printf("\n");
}
if(num2 == 0){
printf("0 0 0.0\n");
}else{
printf("%d",num2);
for(int i=mb-1;i>=0;i--){
if(abs(data1[i])>0.1)printf(" %d %.1lf",i,data1[i]);
}
}
delete data1;
delete data2;
return 0;
}
就这样。
一个集坚强与自信于一身的菇凉。