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