PTA-L2-018 多项式A除以B (25分)
这道题，多项式除法。
先开始我不知道多项除法的操作过程。
其实就是每次要约掉被除数的最高项，依次做下去，知道除数的最高项大于操作过后的被除数的最高项为止。
我们用a[]来保存被除数
b[]来保存除数。
考虑两种情况：
1.如果多项式a的最高项指数小于b的最高项指数，那么这个商就为0，余数就是a
2.如果多项式a的最高项指数大于等于b的最高项指数，那么就可以进行操作，每次商的最高指数就是当前操作的多项式a的最高指数-多项式b的最高指数，相应的可以得到其系数，然后系数累加到我们商的多项式中即可c[]。
然后就是用这个得到的当前的商（指数+系数部分)与多项式b[]相乘，用当前多项式a[]减去这个积，更新多项式a[]，其实就是每次我们的余数多项式部分，一直循环到a的最高项指数小于b的最高项指数即可。
最后就是输出部分，
商存在c[]中
余数存在a[]中。
我们统计一下有多少个非零的需要输出的项，这里注意一下，我们保留一位小数输出，那么临界值就是0.05，大于等于0.05的可以输出，小于0.05的四舍五入后都看为0，不输出。
注意空格部分即可。
代码部分：
#include <bits/stdc++.h>
#define mst(a, n) memset(a, n, sizeof(a))
using namespace std;
const int N = 1e5 + 10;
const int M = 55;
const int INF = 1e6 + 10;
const double eps = 0.05;
typedef long long ll;
double a[N], b[N];
double c[N];
int n, m;
int max1, max2;
double k1, k2;
int main()
{
int i;
int x, y;
cin >> n;
for (i = 1; i <= n; i++)
{
scanf ("%d%d", &x, &y);
a[x] = y;
if (x > max1)
{
max1 = x;
}
}
cin >> m;
for (i = 1; i <= m; i++)
{
scanf ("%d%d", &x, &y);
b[x] = y;
if (x > max2)
{
max2 = x;
}
}
if (max2 > max1)
{
cout << "0 0 0.0\n";
cout << n;
for (i = max1; i >= 0; i--)
{
if (fabs(a[i]) >= eps)
{
printf (" %d %.1f", i, a[i]);
}
}
cout << endl;
}
else
{
int t = max1;
while (max1 >= max2)
{
double tt = a[max1] / b[max2];
c[max1 - max2] += tt;
for (i = max2; i >= 0; i--)
{
if (fabs(b[i]) >= eps)
{
a[i + max1 - max2] -= tt * b[i];
}
}
int maxx = 0;
for (i = t; i >= 0; i--)
{
if (fabs(a[i]) >= eps)
{
if (i > maxx)
{
maxx = i;
break;
}
}
}
if (i < 0)
{
break;
}
max1 = maxx;
}
int cnt1 = 0;
int cnt2 = 0;
for (i = t; i >= 0; i--)
{
if (fabs(c[i]) >= eps)
{
cnt1++;
}
if (fabs(a[i]) >= eps)
{
cnt2++;
}
}
if (!cnt1)
{
cout << "0 0 0.0\n";
}
else
{
int flag = 0;
cout << cnt1 << " ";
for (i = t; i >= 0; i--)
{
if (fabs(c[i]) >= eps)
{
if (flag)
{
cout << " ";
}
flag = 1;
printf ("%d %.1f", i, c[i]);
}
}
cout << endl;
}
if (!cnt2)
{
cout << "0 0 0.0\n";
}
else
{
int flag = 0;
cout << cnt2 << " ";
for (i = t; i >= 0; i--)
{
if (fabs(a[i]) >= eps)
{
if (flag)
{
cout << " ";
}
flag = 1;
printf ("%d %.1f", i, a[i]);
}
}
cout << endl;
}
}
return 0;
}