多项式A除以B（C++）--- 模拟法解题

这仍然是一道关于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

 

 

---

 

 

因为输入格式中保证所有指数是各不相同的非负整数，然后又说输出格式与输入格式相同，就说明不会出现如：x^2 / x^3 这样的情况，因为这样输出的指数为负数了。

然后多项式相除，也就是长除，可以使用竖式进行演算，如对题给样例进行竖式除法演示，演示如下图所示：

 

 

代码如下

#include <iostream>
#include <cmath>
#include <cstdio>
using namespace std;
const int maxn = 1e6 + 10;
int b[maxn]; //b[]用来存储除数各个项的指数 
double a[maxn], c[maxn], d[maxn]; //a[i]=x表示被除数指数为i的项的系数为x，后面还可以用来表示除完后的余 
//c[]用来存储除数各个项的系数  d[i]=x表示 商中指数为i的项 并且系数为x 

int main()
{
	int n, m, limit = -1, cnt1 = 0, cnt2 = 0, exp;
	cin >> n;
	for(int i = 0; i < n; i++) {
		cin >> exp;
		cin >> a[exp];
		if (limit==-1) limit = exp;
	}
	cin >> m;
	for(int i = 0; i < m; i++) cin >> b[i] >> c[i];
	
	int exp_temp;
	for(int i = limit; i >= b[0]; i--) { //i>=b[0]，这一步确保余R的阶数小于B的阶数
		exp_temp = i - b[0];
		d[exp_temp] = a[i] / c[0];
		for(int j = 0; j < m; j++) a[exp_temp+b[j]] = a[exp_temp+b[j]] - c[j] * d[exp_temp];
	}
	
	/*分别统计系数不为0的商和余的项的个数*/ 
	for(int i = limit; i >= 0; i--) {
		if (fabs(d[i])>=0.05) cnt1++;
		if (fabs(a[i])>=0.05) cnt2++;
	}
	cout << cnt1;
	for(int i = limit; i >= 0; i--) {
		if (fabs(d[i])>=0.05) printf(" %d %.1f",i,d[i]);
	}
	if (cnt1==0) printf(" 0 0.0");
	
	cout << endl;
	
	cout << cnt2;
	for(int i = limit; i >= 0; i--) {
		if (fabs(a[i])>=0.05) printf(" %d %.1f",i,a[i]);
	}
	if (cnt2==0) printf(" 0 0.0");
	
	return 0;
}