这仍然是一道关于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
题目分析: 模拟多项式除法过程更新除数和余数的指数和系数,详解看柳神博客L2-018. 多项式A除以B -PAT团体程序设计天梯赛GPLT
题目解答:
#include<bits/stdc++.h>
using namespace std;
double c1[3000],c2[3000],c3[3000];
int count(double c[],int s){
int cnt=0;
for(int i=s;i>=0;i--){
if(abs(c[i])+0.05>=0.1) cnt++;
}
return cnt;
}
void polymult(double c[],int s){
if(count(c,s))
cout<<count(c,s);
else cout<<"0 0 0.0";
for(int i=s;i>=0;i--){
if(abs(c[i])+0.05>=0.1)
printf(" %d %.1lf",i,c[i]);
}
}
int main(){
int n;
cin>>n;
int max1=0,max2=0;
for(int i=0;i<n;i++){
int a;
double b;
cin>>a>>b;
c1[a]=b;
if(a>max1){
max1=a;
}
}
int k;
cin>>k;
for(int i=0;i<k;i++){
int a;
double b;
cin>>a>>b;
c2[a]=b;
if(a>max2){
max2=a;
}
}
int t1=max1;
int t2=max2;
while(t1>=t2){
double c=c1[t1]/c2[t2];
c3[t1-t2]=c;
for(int i=t1,j=t2;j>=0;j--,i--) c1[i]-=c2[j]*c;
while(abs(c1[t1])<0.000001) t1--;
}
polymult(c3,max1-max2);
cout<<endl;
polymult(c1,t1);
return 0;
}