这仍然是一道关于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
思路:
模拟二项式除法,数据存在数组A,B里面**,下标作为指数,内容作为系数,**每次寻找指数最高位进行匹配,确保每次可以消掉A的最高位。
坑点:
1.注意是绝对值大于0.1,绝对值不大于0.1的值作为0舍去。
2.不止是余数可能为0,商也可能为0。
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int MAX = 100000;
int aNum, bNum;
float A[MAX];
float B[MAX];
float out[MAX];
int aMax;
int bMax;//这个是固定的,判断循环是否继续
void chuyi() {//A除以B一次,并将商放在out中
int k = aMax - bMax;
out[k] = A[aMax] / B[bMax];//放进去一位
for (int i = aMax; i>=aMax-bMax ; i--) {//A减去out[k]倍的B
A[i] =A[i]- B[i-k]* out[k];
}
return;
}
int main() {
cin >> aNum;
for (int i = 0; i < aNum; i++) {
int index;//指数
int temp;//系数
cin >> index >> temp;
A[index] = temp;
if (i == 0)
aMax = index;
}
cin >> bNum;
for (int i = 0; i < bNum; i++) {
int index;//指数
int temp;//系数
cin >> index >> temp;
B[index] = temp;
if (i == 0)
bMax = index;
}
//数据读入完成
while (aMax >= bMax) {
chuyi();//除一次
int i;
bool flag = false;
for (i = aMax; i >= 0; i--) {
if (abs(A[i]) >= 0.1) {
aMax = i;
flag = true;
break;
}
}
if (!flag) {//被除尽了
break;
}
if (aMax == -1)//数据处理完成
break;
}
//输出
aNum = 0;
int outNum = 0;
for (int i = 0; i < MAX; i++) {
if (abs(A[i]) >= 0.1)
aNum++;
if (abs(out[i]) >= 0.1)
outNum++;
}
cout << outNum;
for (int i = MAX - 1; i >= 0; i--) {
if(abs(out[i])>=0.1)
printf(" %d %.1f", i, out[i]);
}
if (outNum == 0) {
printf(" 0 0.0");
}
cout << endl;
cout << aNum;
for (int i = MAX-1; i >=0; i--) {
if (abs(A[i]) >= 0.1)
printf(" %d %.1f", i, A[i]);
}
if (aNum == 0) {
printf(" 0 0.0");
}
return 0;
}