团体程序设计天梯赛 多项式A除以B

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

解题思路
比赛的时候一点思路都没有，然后下午突然想到自己高中一直在用的多项式除法，已无力吐槽自己迟钝的脑细胞[是不是讲的太委婉了]。然后就开始敲直到现在改代码改到想吐血………..
首先当你知道多项式除法的规则之后，就可以确定这是模拟除法了，我是直接数组暴力过的，有好方法请戳我(⊙o⊙)哦。
将 A 和 B 多项式的系数分别放置在数组 A[] 和 B[] 中，数组的下标代表每一项的指数（初始化为0），因为除数 B[] 整个过程不可变，所以开辟数组 C[] 来存储每次商与除数的乘积，数组 ans[] 存储商，数组下标同样代表指数。便可通过 A[] 与 C[] 同指数相加减来更新数组 A[]，即更新为每一步的余数，以进行下一步相除运算操作。直至 A[] 中的最大指数小于 B[] 的最大指数相等，结束相除，此时数组 A[] 为最后的余数。

代码实现

#include <iostream>
#include<cstdio>
#include<cmath>
#include<stdlib.h>
#include<cstring>
using namespace std;
double A[1000],C[1000],ans[1000];
int B[1000];
int judge(double x)
{
    if(abs((int)x)>0||abs((int)(x*100)%100)>=5||abs((int)(x*10)%10)>0)
        return 1;
    else
        return 0;
}
int main()
{
    int p,q,i,j,mark;
    int numA,numB;
    int x,y;
    int maxA=0,maxB=0;
    memset(A,0,sizeof(A));
    memset(B,0,sizeof(B));
    memset(C,0,sizeof(C));
    memset(ans,0,sizeof(ans));
    scanf("%d",&numA);
    while(numA--)
    {
        scanf("%d %d",&x,&y);
        A[x]=(double)y;
        maxA=max(maxA,x);
    }
    int t=maxA;
    scanf("%d",&numB);
    while(numB--)
    {
        scanf("%d %d",&x,&y);
        B[x]=y;
        maxB=max(maxB,x);
    }
    while(maxA>=maxB)
    {
        if(A[maxA]==0)break;
        for(i=maxA; i>=0; i--)
        {
            if(A[maxA]!=0)
            {
                p=maxA,q=maxB;
                double xi=(double)A[p]/(double)B[q];
                ans[p-q]=xi;
                for(j=maxB; j>=0; j--)
                {
                    C[j+p-q]=B[j]*xi;
                }
                for(j=0; j<=maxA; j++)
                {
                    A[j]=A[j]-C[j];
                    if(A[j]!=0)
                        mark=j;
                }
                maxA=mark;
            }
            break;
        }
    }
    int countt=0;
    for(i=t; i>=0; i--)
        if(judge(ans[i]))
            countt++;
    if(countt==0)
        printf("0 0 0.0\n");
    else
    {
        printf("%d",countt);
        for(i=t; i>=0; i--)
            if(judge(ans[i]))
                printf(" %d %.1lf",i,ans[i]);
        printf("\n");
    }
    countt=0;
    for(i=maxA; i>=0; i--)
        if(judge(A[i]))
            countt++;
    if(countt==0)
        printf("0 0 0.0\n");
    else
    {
        printf("%d",countt);
        for(i=maxA; i>=0; i--)
            if(judge(A[i]))
                printf(" %d %.1lf",i,A[i]);
        printf("\n");
    }
    return 0;
}