PTA-L2-018 多项式A除以B (25分)

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;
}