【GPLT 二阶题目集】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 <iostream>
#define MAXSIZE 10000
using namespace std;

double sswr(double x) //四舍五入
{
    if (x > -0.05 && x < 0.05)
        return 0;
    return x;
}

int main()
{
    double c1[MAXSIZE] = { 0 }, c2[MAXSIZE] = { 0 };
    double cshang[MAXSIZE] = { 0 }; //下标为指数，值为系数
    int e1_max, e2_max; //最大指数
    int n; cin >> n;
    if (n == 0) //某多项式为零项式，直接输出结果
    {
        printf("0 0 0.0\n0 0 0.0");
        return 0;
    }
    else
    {
        for (int i = 0; i < n; i++)
        {
            int e, c;
            cin >> e >> c;
            if (i == 0)
                e1_max = e;
            c1[e] = c;
        }
    }
    cin >> n;
    if (n == 0)
    {
        printf("0 0 0.0\n0 0 0.0");
        return 0;
    }
    else
    {
        for (int i = 0; i < n; i++)
        {
            int e, c;
            cin >> e >> c;
            if (i == 0)
                e2_max = e;
            c2[e] = c;
        }
    }
    bool flag = true;
    int eshang_max = 0; //商的最大指数
    while (e1_max >= e2_max) //被除式（余数）的最大指数需要比除式的大
    {
        //商新增一项
        int zhishu;
        double xishu;
        double ctemp[MAXSIZE] = { 0 };
        zhishu = e1_max - e2_max;
        if (flag)
        {
            eshang_max = zhishu;
            flag = false;
        }
        xishu = c1[e1_max] / c2[e2_max];
        cshang[zhishu] = xishu;
        for (int j = e2_max; j >= 0; j--)
        { //该项与除式相乘得出临时多项式
            ctemp[zhishu + j] = xishu * c2[j];
        }
        for (int j = e1_max; j >= 0; j--)
        { //当前被除式-临时多项式得到新的被除式
            c1[j] -= ctemp[j];
        }
        for (int j = e1_max; j >= 0; j--)
        { //更新被除式的最大指数
            if (c1[j] != 0)
            {
                e1_max = j;
                break;
            }
            if (j == 0)
                e1_max = 0;
        }
    }
    int n_shang = 0, n_yu = 0; //统计商和余数的项数
    for (int i = eshang_max; i >= 0; i--)
    {
        if (sswr(cshang[i]) != 0)
            n_shang++;
        else
            cshang[i] = 0;
    }
    for (int i = e1_max; i >= 0; i--)
    {
        if (sswr(c1[i]) != 0)
            n_yu++;
        else
            c1[i] = 0;
    }
    if (n_shang == 0) //商为0的情况特殊处理
    {
        printf("0 0 0.0");
    }
    else //正常输出商
    {
        printf("%d", n_shang);
        for (int i = eshang_max; i >= 0; i--)
        {
            if (cshang[i] != 0)
                printf(" %d %.1f", i, cshang[i]);
        }
    }
    cout << endl;
    if (n_yu == 0) //余数为0的情况特殊处理
    {
        printf("0 0 0.0");
    }
    else //正常输出余数
    {
        printf("%d", n_yu);
        for (int i = e1_max; i >= 0; i--)
        {
            if (c1[i] != 0)
                printf(" %d %.1f", i, c1[i]);
        }
    }
    cout << endl;
    return 0;
}

 注意事项：

注意四舍五入后再确定商或余数的项数。

刚开始没仔细看题目，没注意到给定多项式是按先A后B的顺序，一直有测试点答案错误。

如有问题，欢迎提出。