PTA L2-018 多项式A除以B

这仍然是一道关于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 <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <cmath>

using namespace std;

const int N = 1000010;

double a[N],b[N];
double s[N],y[N];

void print(double t[],int mt)
{
    int cnt = 0;
    for(int i = mt;i >= 0;i--)
        if(abs(t[i]) + 0.05 > 0.1)//系数不为0
            cnt++;
    if(!cnt)
    {
        puts("0 0 0.0");
        return;
    }
    printf("%d",cnt);
    for(int i = mt;i >= 0;i--)//0.05
        if(abs(t[i]) + 0.05 > 0.1)//系数不为0
            printf(" %d %.1lf",i,t[i]);
    puts("");
}

int main()
{
    int n = 0;
    scanf("%d",&n);
    int ma = 0,mb = 0;
    for(int i = 0; i < n;i++)
    {
        int x = 0,z = 0;
        scanf("%d%d",&z,&x);
        a[z] = x;
        ma = max(ma,z);
    }
    scanf("%d",&n);
    for(int i = 0; i < n;i++)
    {
        int x = 0,z = 0;
        scanf("%d%d",&z,&x);
        b[z] = x;
        mb = max(mb,z);
    }

    int ms = 0;
    while(ma >= 0 && ma >= mb)//指数判断
    {
        int z = ma - mb;
        double x = a[ma] / b[mb];
        s[z] = x;
        ms = max(ms,z);

        for(int i = mb; i >= 0;i--)
        {
            if(abs(b[i]) + 0.05 > 0.1)//系数不为0
            {
                a[i + z] -= b[i] * x;
            }
        }
        
        while(ma >= 0 && abs(a[--ma]) + 0.05 < 0.1);//系数不为0
        
    }
    print(s,ms);
    print(a,ma);
    return 0;
}

结果：