思路:
这道题真的把我的心态搞崩了。中间一些小细节写错了,没注意,调bug调到眼都花了!
(1)把被除式、除式按某个字母作降幂排列,并把所缺的项用零补齐.
(2)用被除式的第一项除以除式第一项,得到商式的第一项.
(3)用商式的第一项去乘除式,把积写在被除式下面(同类项对齐),消去相等项,把不相等的项结合起来.
(4)把减得的差当作新的被除式,再按照上面的方法继续演算,直到余式为零或余式的次数低于除式的次数时为止,被除式=除式×商式+余式。若余式为零,说明这个多项式能被另一个多项式整除
q(x)商 f(x)被除数 g(x)除数 r(x)余
具体看代码~~
#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<math.h>
using namespace std;
double c1[100005], c2[100005], c3[100005];
int f(double a[], int e)//判断多项式非0的项数
{
int ret = 0;
for (int i = e; i >= 0; i--)
{
if (fabs(a[i]) >= 0.05)
ret++;
}
return ret;
}
void solve(double a[], int e)//打印多项式
{
cout << f(a, e);
if (f(a, e) == 0)
cout << " 0 0.0";//零多项式
else
{
for (int i = e; i >= 0; i--)
{
if (fabs(a[i]) >= 0.05)
printf(" %d %.1f", i, a[i]);
}
}
cout << endl;
}
int main()
{
int n;
while (cin >> n)
{
int maxe1, maxe2;//A B的最大指数
for (int i = 0; i < n; i++)
{
int e;
double c;
cin >> e >> c;
if (i == 0)
{
maxe1 = e;
}
c1[e] = c;//指数为e,系数为c
}
int m;
cin >> m;
for (int i = 0; i < m; i++)
{
int e;
double c;
cin >> e >> c;
if (i == 0)
{
maxe2 = e;
}
c2[e] = c;
}
int d = maxe1 - maxe2;
while (maxe2 <= maxe1)
{
double v = c1[maxe1] / c2[maxe2];
c3[maxe1 - maxe2] = v;//商
for (int i = maxe1,j = maxe2; i >= 0 && j >= 0; j--, i--)
{
c1[i] -= c2[j] * v;//余数 指数相同的对应位置系数相减
}
while (fabs(c1[maxe1]) < 0.05)//四舍五入为0的项去掉
maxe1--;
}
solve(c3, d);//打印
solve(c1, maxe1);
}
system("pause");
return 0;
}