这仍然是一道关于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 <iostream>
#define MAXSIZE 10000
using namespace std;
double sswr(double x) //四舍五入
{
if (x > -0.05 && x < 0.05)
return 0;
return x;
}
int main()
{
double c1[MAXSIZE] = { 0 }, c2[MAXSIZE] = { 0 };
double cshang[MAXSIZE] = { 0 }; //下标为指数,值为系数
int e1_max, e2_max; //最大指数
int n; cin >> n;
if (n == 0) //某多项式为零项式,直接输出结果
{
printf("0 0 0.0\n0 0 0.0");
return 0;
}
else
{
for (int i = 0; i < n; i++)
{
int e, c;
cin >> e >> c;
if (i == 0)
e1_max = e;
c1[e] = c;
}
}
cin >> n;
if (n == 0)
{
printf("0 0 0.0\n0 0 0.0");
return 0;
}
else
{
for (int i = 0; i < n; i++)
{
int e, c;
cin >> e >> c;
if (i == 0)
e2_max = e;
c2[e] = c;
}
}
bool flag = true;
int eshang_max = 0; //商的最大指数
while (e1_max >= e2_max) //被除式(余数)的最大指数需要比除式的大
{
//商新增一项
int zhishu;
double xishu;
double ctemp[MAXSIZE] = { 0 };
zhishu = e1_max - e2_max;
if (flag)
{
eshang_max = zhishu;
flag = false;
}
xishu = c1[e1_max] / c2[e2_max];
cshang[zhishu] = xishu;
for (int j = e2_max; j >= 0; j--)
{ //该项与除式相乘得出临时多项式
ctemp[zhishu + j] = xishu * c2[j];
}
for (int j = e1_max; j >= 0; j--)
{ //当前被除式-临时多项式得到新的被除式
c1[j] -= ctemp[j];
}
for (int j = e1_max; j >= 0; j--)
{ //更新被除式的最大指数
if (c1[j] != 0)
{
e1_max = j;
break;
}
if (j == 0)
e1_max = 0;
}
}
int n_shang = 0, n_yu = 0; //统计商和余数的项数
for (int i = eshang_max; i >= 0; i--)
{
if (sswr(cshang[i]) != 0)
n_shang++;
else
cshang[i] = 0;
}
for (int i = e1_max; i >= 0; i--)
{
if (sswr(c1[i]) != 0)
n_yu++;
else
c1[i] = 0;
}
if (n_shang == 0) //商为0的情况特殊处理
{
printf("0 0 0.0");
}
else //正常输出商
{
printf("%d", n_shang);
for (int i = eshang_max; i >= 0; i--)
{
if (cshang[i] != 0)
printf(" %d %.1f", i, cshang[i]);
}
}
cout << endl;
if (n_yu == 0) //余数为0的情况特殊处理
{
printf("0 0 0.0");
}
else //正常输出余数
{
printf("%d", n_yu);
for (int i = e1_max; i >= 0; i--)
{
if (c1[i] != 0)
printf(" %d %.1f", i, c1[i]);
}
}
cout << endl;
return 0;
}
注意事项:
注意四舍五入后再确定商或余数的项数。
刚开始没仔细看题目,没注意到给定多项式是按先A后B的顺序,一直有测试点答案错误。
如有问题,欢迎提出。