这仍然是一道关于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
代码长度限制
16 KB
时间限制
400 ms
内存限制
64 MB
分析:
我们可以根据高等数学中学的多项式除法进行模拟计算,再在里面加一点细节,就能过掉了;
代码:
#include <bits/stdc++.h>
using namespace std;
const int N=3000;
double c1[N],c2[N],c3[N];
int num(double c[],int start)
{
int cnt=0;
for(int i=start;i>=0;i--)
{
if(abs(c[i])+0.05>=0.1) cnt++;
}
return cnt;
}
void PrintPloy(double c[],int start)
{
int cnt=num(c,start);
cout<<cnt;
if(cnt==0) cout<<" 0 0.0";
for(int i=start;i>=0;i--)
if(abs(c[i])+0.05>=0.1)
printf(" %d %.1lf",i,c[i]);
cout<<endl;
}
int main()
{
int n,m;
cin>>n;
int t=0,mx1=-1,mx2=-1;
for(int i=0;i<n;i++)
{
cin>>t>>c1[t];
mx1=max(t,mx1);
}
cin>>m;
for(int i=0;i<m;i++)
{
cin>>t>>c2[t];
mx2=max(mx2,t);
}
int t1=mx1,t2=mx2;
while(t1>=t2)
{
double c=c1[t1]/c2[t2];
c3[t1-t2]=c;
for(int i=t1,j=t2;j>=0;i--,j--)
c1[i]-=c2[j]*c;
while(abs(c1[t1])<0.000001) t1--;
}
PrintPloy(c3,mx1-mx2);
PrintPloy(c1,t1);
return 0;
}