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,这个值不唯一,在一个不是太大也不是太小的区间之内就行。
}
}
}
}