GPLT多项式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

#include<bits/stdc++.h>
#define el '\n'
const int N = 1e5 + 10;
using namespace std;
int main()
{   double a[N],b[N],c[N];
	int n,m,i,j,ma,mb,na=0,nc=0,x;
	cin>>n;
	for(i=0;i<n;i++){
	cin>>x;getchar();cin>>a[x];
	if(!i)ma=x;	
	}
	cin>>m;
	for(i=0;i<m;i++){
	cin>>x;getchar();cin>>b[x];
	if(!i)mb=x;		
	}
	for(i=ma;i>=mb;i--){
		c[i-mb]=a[i]/b[mb];
		for(j=mb;j>=0;j--)
		a[i-mb+j]-=c[i-mb]*b[j];
	}
	for(i=ma-mb;i>=0;i--)
	   if(fabs(c[i])>=0.05)nc++;
	   else c[i]=0;
	for(i=mb-1;i>=0;i--)
	   if(fabs(a[i])>=0.05) na++;
	   else  a[i]=0;
	   
	if(nc==0)cout<<"0 0 0.0"<<el;
	else{
		cout<<nc;
		for(i=ma-mb;i>=0;i--)
        if(fabs(c[i])>1e-6)printf(" %d %.1lf",i,c[i]);
		cout<<el;
	}
	if(na==0)cout<<"0 0 0.0"<<el;
	else{
		cout<<na;
		for(i=mb-1;i>=0;i--)
            if(fabs(a[i])>1e-6)printf(" %d %.1lf",i,a[i]);
		cout<<el;
	}
	return 0;
}