L2-018 多项式A除以B (25 分)

L2-018 多项式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

思路

开始以为是没学过的知识，在OI Wiki上研究了很久多项式求逆和多项式除法，开始怀疑人生觉得自己是不是完全没学会数论。。。

后来看了L2-018 多项式A除以B (25分)的解法才明白这其实就是道模拟，，，就直接模拟竖式除法的运算过程就可以。注意特判零多项式的情况。

代码

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define PI acos(-1)
#define pb push_back
using namespace std;
typedef pair<int, int> P;
typedef long long ll;
const int N = 5e3 + 19;
const ll mod = 1e9 + 7;
const double eps = 1e-12;

double a[N], b[N], c[N];

int main()
{
    int n = 0, m = 0;
    int maxn = -INF, maxm = -INF;
    int pos;
    double num;
    
    scanf("%d", &n);
    for(int i = 0; i < n; i++)
    {
        scanf("%d %lf", &pos, &num);
        a[pos] = num;
        maxn = max(maxn, pos);
    }
    scanf("%d", &m);
    for(int i = 0; i < m; i++)
    {
        scanf("%d %lf", &pos, &num);
        b[pos] = num;
        maxm = max(maxm, pos);
    }
    
    int len = maxn - maxm;
    
    while(maxn >= maxm)
    {
        double num = a[maxn] / b[maxm];//
        c[maxn - maxm] = num;
        for(int i = maxn, j = maxm; i >= 0 && j >= 0; i--, j--)
        {
            a[i] -= num * b[j];
        }
        while(abs(a[maxn]) < eps)
        {
            maxn--;
        }
    }
    
    int cnt = 0;
    for(int i = 0; i <= len; i++)
        if(abs(c[i]) + 0.05 >= 0.1)
            cnt++;
    printf("%d", cnt);
    if(cnt == 0)
    {
        printf(" 0 0.0");
    }
    for(int i = len; i >= 0; i--)
    {
        if(abs(c[i]) + 0.05 >= 0.1)
        {
            printf(" %d %.1f", i, c[i]);
        }
    }
    puts("");
    
    cnt = 0;
    for(int i = 0; i < maxm; i++)
        if(abs(a[i]) + 0.05 >= 0.1)
            cnt++;
    printf("%d", cnt);
    if(cnt == 0)
    {
        printf(" 0 0.0");
    }
    for(int i = maxm - 1; i >= 0; i--)
    {
        if(abs(a[i]) + 0.05 >= 0.1)
        {
            printf(" %d %.1f", i, a[i]);
        }
    }
    puts("");
    
    return 0;
}