团体程序设计天梯赛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

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
#include <stack>
#include <queue>
#include <climits>
#include <functional>
#include <vector>
#include <map>
#include <set>

using namespace std;

#define LL long long

double c[100050], d[100050], ans[100050], e[100050];
int n, m, ma1, ma2;

int main()
{
    int x;
    double v;
    scanf("%d",&n);
    for(int i=0;i<=100000;i++) c[i]=d[i]=0;
    for(int i=1;i<=n;i++)
    {
        scanf("%d%lf",&x,&v), c[x]=v;
        if(i==1) ma1=x;
    }
    scanf("%d",&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%lf",&x,&v), d[x]=v;
        if(i==1) ma2=x;
    }
    for(int i=1;i<=ma1-ma2+1;i++)
    {
        int k=-1;
        for(int j=100000;j>=ma2;j--)
            if(c[j])
            {
                k=j;
                break;
            }
        if(k==-1) break;
        ans[k-ma2]=c[k]/d[ma2];
        double temp=ans[k-ma2];
        for(int j=ma1;j>=0;j--) e[j]=0;
        for(int j=ma2;j>=0;j--) e[j+k-ma2]=temp*d[j];
        for(int j=100000;j>=0;j--) c[j]-=e[j];
    }
    int sum=0;
    for(int i=100000;i>=0;i--)
        if(fabs(ans[i])>=0.05) sum++;
    if(sum==0) printf("0 0 0.0\n");
    else
    {
        printf("%d",sum);
        for(int i=100000;i>=0;i--)
            if(fabs(ans[i])>=0.05) printf(" %d %.1lf",i,ans[i]);
        printf("\n");
    }
    sum=0;
    for(int i=100000; i>=0;i--)
        if(fabs(c[i])>=0.05) sum++;
    if(sum==0) printf("0 0 0.0\n");
    else
    {
        printf("%d",sum);
        for(int i=100000;i>=0;i--)
            if(fabs(c[i])>=0.05) printf(" %d %.1lf",i,c[i]);
        printf("\n");
    }
    return 0;
}