**
7-9 多项式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
- 做这种题有种返璞归真的感觉,有点理解为陈景润的伟大了-.-。
- 大概思路就是用一个数组来存储,数组的下标代表次幂。
- 然后最高次幂挨个相除后再两个数相减如此循环到最高次幂小于除数的最高次幂。
#include<bits/stdc++.h>
using namespace std;
const int maxn=3e3+10;
double c1[maxn],c2[maxn],c3[maxn];
int nonNegativeNum(double c[],int st)//统计非负项个数
{
int cnt=0;
for(int i=st;i>=0;--i)
if(abs(c[i])>=0.05)cnt++;
return cnt;
}
void printPoly(double c[],int st)
{
printf("%d",nonNegativeNum(c,st));
if(nonNegativeNum(c,st)==0)printf(" 0 0.0");
for(int i=st;i>=0;--i)
if(abs(c[i])+0.05>=0.1)printf(" %d %.1lf",i,c[i]);
}
int main()
{
int max1=-1,max2=-1;
int m;
scanf("%d",&m);
for(int i=0;i<m;++i)
{
int t;
scanf("%d",&t);
max1=max(max1,t);
scanf("%lf",&c1[t]);
}
int n;
scanf("%d",&n);
for(int i=0;i<n;++i)
{
int t;scanf("%d",&t);
max2=max(max2,t);//max函数的使用
scanf("%lf",&c2[t]);
}
int t1=max1,t2=max2;
while(t1>=t2)
{
double x=c1[t1]/c2[t2];//最高次幂的商的系数
c3[t1-t2]=x;
for(int i=t1,j=t2;j>=0;--j,--i)
c1[i]-=c2[j]*x;
while(abs(c1[t1])<1e-6) t1--;//如果该项是0,那么最高次幂降1;
}
printPoly(c3,max1-max2);
printf("\n");
printPoly(c1,t1);
return 0;
}