PTA多项式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
代码长度限制
16 KB
时间限制
400 ms
内存限制
64 MB

答案:

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

#define MAX 3000//注意此处如果太小了，就会出现数组越界的问题，从而导致段错误。
#define ZERO 1E-6

int Digit(double p[], int start);
void Print(double p[], int start);
int Read(double p[]);
double* Calc(double p1[], double p2[], int max1, int max2);

int main()
{
	double p1[MAX], p2[MAX], * p;
	int max1, max2;
	max1 = Read(p1);
	max2 = Read(p2);
	p = Calc(p1, p2, max1, max2);
	Print(p, max1 - max2);
	printf("\n");
	Print(p1, max1);
	free(p);
	return 0;
}

int Digit(double p[], int start)
{
	int cnt = 0;
	for (int i = start; i >= 0; i--)
	{
		if (fabs(p[i]) + 0.05 >= 0.1)cnt++;//注意千万不要写成if (fabs(p[i] + 0.05 >= 0.1)),否则错误很难被发现。
	}
	return cnt;
}

void Print(double p[], int start)
{
	int digit;
	digit = Digit(p, start);
	printf("%d", digit);
	if (digit == 0)printf(" 0 0.0");
	for (int i = start; i >= 0; i--)
	{
		if (fabs(p[i]) + 0.05 >= 0.1)printf(" %d %.1f", i, p[i]);
	}
}

int Read(double p[])
{
	int n, e;//expon
	double coef;
	int max = -1;
	scanf("%d", &n);
	for (int i = 0; i < n; i++)
	{
		scanf("%d %lf", &e, &coef);
		p[e] = coef;
		max = (e > max) ? e : max;
	}
	return max;
}

double* Calc(double p1[], double p2[], int max1, int max2)
{
	int e1 = max1, e2 = max2;
	double coef;
	int i, j;
	double* p = (double*)malloc((max1 + 1) * sizeof(double));
    //注意此处，如果用double p[MAX],则在函数结束时，编译器会自动释放内存。
    //而动态分配内存的话，编译器不会自动释放，需要手动释放。

	while (e1 >= e2)
	{
		coef = p1[e1] / p2[e2];
		p[e1 - e2] = coef;
		for (i = e1, j = e2; j >= 0; i--, j--)
		{
			p1[i] -= p2[j] * coef;
		}
		while (fabs(p1[e1]) < ZERO)e1--;
	}
	return p;
}

不过让我疑问的一点就是，为什么此程序在VS2022里得不到正确的输出，这个问题留着以后再解答吧。