7-9 多项式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

题目分析： 模拟多项式除法过程更新除数和余数的指数和系数，详解看柳神博客L2-018. 多项式A除以B -PAT团体程序设计天梯赛GPLT

题目解答：

#include<bits/stdc++.h>
using namespace std;
double c1[3000],c2[3000],c3[3000];
int count(double c[],int s){
    int cnt=0;
    for(int i=s;i>=0;i--){
        if(abs(c[i])+0.05>=0.1) cnt++;
    }
    return cnt;
}
void polymult(double c[],int s){
    if(count(c,s))
    cout<<count(c,s);
    else cout<<"0 0 0.0";
    for(int i=s;i>=0;i--){
        if(abs(c[i])+0.05>=0.1)
        printf(" %d %.1lf",i,c[i]);
    }
}
int main(){
    int n;
    cin>>n;
    int max1=0,max2=0;
    for(int i=0;i<n;i++){
        int a;
        double b;
        cin>>a>>b;
        c1[a]=b;
        if(a>max1){
            max1=a;
        }
    }
    int k;
    cin>>k;
    for(int i=0;i<k;i++){
        int a;
        double b;
        cin>>a>>b;
        c2[a]=b;
        if(a>max2){
            max2=a;
        }
    }
    int t1=max1;
    int t2=max2;
    while(t1>=t2){
        double c=c1[t1]/c2[t2];
        c3[t1-t2]=c;
        for(int i=t1,j=t2;j>=0;j--,i--) c1[i]-=c2[j]*c;
        while(abs(c1[t1])<0.000001) t1--;
    }
    polymult(c3,max1-max2);
    cout<<endl;
    polymult(c1,t1);
    return 0;
}