L2-018 多项式A除以B (25分)（C语言）

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

解题思路：

我们设置三个数组，分别为abc，用来模拟多项式ab以及最后的商，用数组的下标来表示阶数，而最后的余数就是多项式a最后无法再除时所剩余的式子，我们就用这些数组来模拟这个除法的过程。
这题其实不难，主要是不知道什么是多项式除法，如果我写的过程看不懂，可以去自己再搜一下。

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>
//长头发啦长头发啦
int main(void)
{
    double a[10010],b[10010],c[10010];
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(c,0,sizeof(c));
    int n,m,x,i=0,tmp1=0,tmp2=0;
    double y;//记录商的每一项的系数
    scanf("%d",&n);
    for(int i=0; i<n; i++)
    {
        scanf("%d%lf",&x,&y);
        a[x]+=y;
        tmp1=tmp1>x?tmp1:x;//找到最大的阶数
    }
    scanf("%d",&m);
    for(int i=0; i<m; i++)
    {
        scanf("%d%lf",&x,&y);
        b[x]+=y;
        tmp2=tmp2>x?tmp2:x;
    }
    if(tmp1<tmp2)//不能省略，不然测试点二和四过不了
    {
        printf("0 0 0.0\n");
        printf("%d",n);
        for(int i=tmp1; i>=0; i--)
        {
            if(fabs(a[i])>0.05)
                printf(" %d %.1lf\n",i,a[i]);
        }
        printf("\n");
    }
    else
    {
        int tp=tmp1;
        while(tmp1>=tmp2)
        {
            y=a[tmp1]/b[tmp2];//商的系数，阶数用i+tmp1-tmp2来表示
            c[tmp1-tmp2]+=y;
            for(i=tmp2; i>=0; i--)
            {
                if(fabs(b[i])>0.05)
                    a[i+tmp1-tmp2]-=y*b[i];
            }
            int ep=0;
            for(i=tmp1; i>=0; i--)
            {
                if(fabs(a[i])>0.05)
                {
                    ep=i;
                    break;
                }
            }
            tmp1=ep;//除了一次，所以阶数要变了
        }
        int count1=0,count2=0;
        for(int i=tp; i>=0; i--)
        {
            if(fabs(a[i])>0.05)
                count1++;//a多项式计数器，计数不为零的项
        }
        for(int i=tp; i>=0; i--)
        {
            if(fabs(c[i])>0.05)
                count2++;//商计数器，计算不为零的项
        }
        int i=0;
        if(count2==0)//如果商为零
            printf("0 0 0.0\n");
        else
        {
            printf("%d",count2);//商不为零，按照递减顺序输出不为零的所有项
            for(count2=0,i=tp; i>=0; i--)
            {
                if(count2!=0)
                {
                    count2++;
                    printf(" ");
                }
                if(fabs(c[i])>0.05)
                    printf(" %d %.1lf",i,c[i]);
            }
            printf("\n");
        }
        if(count1==0)//如果余数为零
            printf("0 0 0.0\n");
        else
        {
            printf("%d",count1);//余数不为零，同上
            for(count1=0,i=tp; i>=0; i--)
            {
                if(count1!=0)
                {
                    printf(" ");
                    count1++;
                }
                if(fabs(a[i])>0.05)
                    printf(" %d %.1lf",i,a[i]);//因为精度问题，就像本题之中的1/27，所以题中的所有>0全部替换成大于0.05，这个值不唯一，在一个不是太大也不是太小的区间之内就行。
            }

        }
    }
}