这仍然是一道关于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
具体过程参考代码注释
完整代码:
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
double a[maxn],b2[maxn],c[maxn];
int b1[maxn];
int main()
{
// ios::sync_with_stdio(false);
// cin.tie(0);
int n,first;
cin>>n;
for(int i=0;i<n;i++){
int x;
cin>>x;
cin>>a[x];//a[x]表示第一个多项式指数为x的项的系数
if(!i) first=x;//first记录第一项的指数,即:多项式的最高次幂
}
int m;
cin>>m;
for(int i=0;i<m;i++) cin>>b1[i]>>b2[i];//b1[i]和b2[i]分别表示第二个多项式第i项的指数和系数
for(int i=first;i>=b1[0];i--){,
int t=i-b1[0];//t为当前被除数和除数的最高次幂之差,即:当前商的指数
c[t]=a[i]/b2[0];//c[t]为商中指数为t的这一项的系数(因为每次都是让除数*商消掉被除数的最高次幂那一项,所以商的系数=除数的系数/被除数的系数)
for(int j=0;j<m;j++){
a[t+b1[j]]-=c[t]*b2[j];//得出商之后,被除数的每一项都要减去商*除数作为下一次的被除数,不能再除时,a[i]就保存了余数中指数为i的系数
}
}
int cnt1=0,cnt2=0;//cnt1和cnt2统计商和余数中非零项的个数
for(int i=first;i>=0;i--){
if(fabs(c[i])>=0.05) cnt1++;
if(fabs(a[i])>=0.05) cnt2++;
}
cout<<cnt1;
for(int i=first;i>=0;i--){
if(fabs(c[i])>=0.05){
printf(" %d %.1f",i,c[i]);
}
}
if(cnt1==0) cout<<" 0 0.0";
cout<<endl;
cout<<cnt2;
for(int i=first;i>=0;i--){
if(fabs(a[i])>=0.05){
printf(" %d %.1f",i,a[i]);
}
}
if(cnt2==0) cout<<" 0 0.0"<<endl;
return 0;
}