PTA 多项式A除以B （25 分）

7-10 多项式A除以B （25 分）
这仍然是一道关于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>
using namespace std;
const int N=1e5+5;
int b[N];
double a[N],c[N],d[N];
int main()
{
    int n,m,f1,n1=0,n2=0;
    scanf("%d",&n);
    for(int i=0,x; i<n; i++)
    {
        scanf("%d",&x);
        scanf("%lf",&a[x]);
        if(!i)f1=x;
    }
    scanf("%d",&m);
    for(int i=0; i<m; i++)scanf("%d%lf",&b[i],&d[i]);
    for(int i=f1,t; i>=b[0]; i--)
    {
        t=i-b[0];
        c[t]=a[i]/d[0];
        for(int j=0; m>j; j++)a[t+b[j]]=a[t+b[j]]-c[t]*d[j];
    }
    for(int i=f1; i>=0; i--)
    {
        if(fabs(a[i])>=0.05)n2++;
        if(fabs(c[i])>=0.05)n1++;
    }
    printf("%d",n1);
    for(int i=f1; i>=0; i--)if(fabs(c[i])>=0.05)printf(" %d %.1f",i,c[i]);
    if(n1==0)printf(" 0 0.0");
    printf("\n");
    printf("%d",n2);
    for(int i=f1; i>=0; i--)if(fabs(a[i])>=0.05)printf(" %d %.1f",i,a[i]);
    if(n2==0)printf(" 0 0.0");
    return 0;
}