L2-018. 多项式A除以B

L2-018. 多项式A除以B

时间限制

400 ms

内存限制

65536 kB

代码长度限制

8000 B

判题程序

Standard

作者

陈越

这仍然是一道关于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

---

提交代码

大体多项式除法的思路是，先用a多项式的最高次数项除b多项式的最高次数项，得到一项多项式，然后用a多项式减去这一项乘b

更新a继续除，如果a的最高次数项的系数小于0.05就被认作系数为零这项没有了，就要更新a的最高次数项，继续除

具体操作就是先存好a，b两个多项式，并记录下。如果a的最高次数大于等于b的最高次数就执行循环（如果b的最高次数都大于a的最高次数了，那么a除b肯定是0了），每次循环中，算出a，b最高次项系数比，指数比，如果求出来的这个系数比大于等于0.05那么就可以作为商的一项存入答案数组，并且更新a（更新过程就是a-这一项乘b）如果求出来的系数比小于0.05，就需要减少a的最高次数；然后在保证a的最高次数大于等于b的最高次数的前提下找到系数大于等于0.05的那一项，并把这个作为a的新的最高次数，完成新的a的最高次数的更新；直到循环结束

如果答案数组个数为零，就是零多项式，否则正常输出

code：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <map>
#include <cmath>
using namespace std;
struct poly{
    int e;
    double c;
}p[10000],ans[10000];//p记录除数多项式，ans记录商或余数

int main(){
    int n,m,ant = 0,e;//m记录被除数多项式指数的最高次
    double c;
    map<int,double>q;
    cin >> n;
    for(int i = 0; i < n; i++){
        cin >> e >> c;
        q[e] = c;
        if(i == 0) m = e;
    }
    cin >> n;
    for(int i = 0; i < n; i++){
        cin >> p[i].e >> p[i].c;
    }
    while(m >= p[0].e){
        double change = q[m] / p[0].c;//系数bi
        int diff = m - p[0].e;//指数比
        if(fabs(change) >= 0.05){
            ans[ant].e = diff;
            ans[ant++].c = change;
            for(int i = 0; i < n; i++){//更新a多项式
                q[p[i].e + diff] -= change * p[i].c;//通过b多项式加上指数比去找a多项式的项进行减法
            }
        }
        else m--;
        while(m >= p[0].e && fabs(q[m]) < 0.05){//这一步的作用是防止上一步m--后的项虽然次数减少了，可能系数小于0.05
            m--;
        }
    }
    cout << ant;
    if(!ant) cout << " 0 0.0";
    for(int i = 0; i < ant; i++){
        printf(" %d %.1f",ans[i].e,ans[i].c);
    }
    cout << endl;
    ant = 0;
    while(m >= 0){
        if(fabs(q[m]) >= 0.05){
            ans[ant].e = m;
            ans[ant++].c = q[m];
        }
        m--;
    }
    cout << ant;
    if(!ant) cout << " 0 0.0";
    for(int i = 0; i < ant; i++){
        printf(" %d %.1f",ans[i].e,ans[i].c);
    }
    return 0;
}