多项式A除以B[天梯赛]

输入样例：
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 <iostream>
#include <cstdio> 
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;
const int N = 1000010;
double A[N], B[N]; //A[i]表示：多项式A 指数为i的一项系数为A[i];B[i]同理 
double C[N]; //存放商的结果 C[i]表示：多项式C 指数为i的一项系数为C[i]
int A_max_e, B_max_e; //多项式A和B的最高次幂 
int main()
{
	int n;
	int e, c; //e为当前项的指数，c为当前项的系数
	scanf("%d", &n);
	for(int i = 0; i < n; i ++)
	{
		scanf("%d%d", &e, &c);
		if(!i) A_max_e = e; //题目规定指数按递减顺序给出，因此第一个为最高次幂 
		A[e] = c; 
	}
	int m;
	scanf("%d", &m);
	for(int i = 0; i < m; i ++)
	{
		scanf("%d%d", &e, &c);
		if(!i) B_max_e = e;
		B[e] = c; 
	}
	
	//模拟多项式除法
	/*
	 每次将被除数(A)最高次幂的一项消掉，得到商的一部分
	 该部分的
	 	指数：被除数(A)最高次幂 - 除数(B)的最高次幂
		底数：被除数(A)最高次幂对应的系数 / 除数(B)的最高次幂对应的系数
	 之后将该部分 * 除数， 再让被除数(A)减去， 将最后的结果保存在被除数(A)中---即被除数(A)始终变化，除数(B)不变化 
	*/
	for(int i = A_max_e; i >= B_max_e; i --) //当被除数的最高次幂 < 除数最高次幂 时，不在有商，此时的被除数即为余数 
	{
		C[i - B_max_e] = A[i] / B[B_max_e];
		
		/*
			更新被除数
			C[i - B_max_e] * B[j]为商 * 除数的结果 
			遍历范围为 B_max_e ~ 0, 对应除数的指数范围，即使除数没有对应的指数的一项，那么数组C对应的值为0，不影响结果
			A数组的变化为 i - B_max_e + j, i为当前被除数(A)的最高次幂，则被除数的指数范围最大为为i ~ 0,因此有 - B_max_e + j 的偏移量 
		*/
		for(int j = B_max_e; j >= 0; j --) 
			A[i - B_max_e + j] -= C[i - B_max_e] * B[j];  
	}
	int cnt1 = 0, cnt2 = 0; //商和余数的项数
	
	//商
	for(int i = A_max_e - B_max_e; i >= 0; i --) //商的最高次幂为第一次操作得到的那部分
	{
		if(abs(C[i]) >= 0.1) cnt1 ++; //题目要求保留一位小数，因此小于0.1的视为0，题目要求不输出 
	} 
	
	//余数
	for(int i = B_max_e - 1; i >= 0; i --) //余数中包含的最高次幂最大为 B_max_e(除数最高次幂) - 1 
	{
		if(abs(A[i]) >= 0.1) cnt2 ++;
	} 
	
	//输出结果
	if(!cnt1) printf("0 0 0.0\n");
	else
	{
		printf("%d", cnt1);
		for(int i = A_max_e - B_max_e; i >= 0; i --)
		{
			if(abs(C[i]) >= 0.1) printf(" %d %.1f", i, C[i]); 
		}	
		printf("\n");
	} 
	
	if(!cnt2) printf("0 0 0.0\n");
	else
	{
		printf("%d", cnt2);
		for(int i = B_max_e - 1; i >= 0; i --)
		{
			if(abs(A[i]) >= 0.1) printf(" %d %.1f", i, A[i]); 
		}	
		printf("\n");
	} 
	return 0; 
}

欢迎批评指正！！！