1002. A+B for Polynomials

本文介绍了一种解决多项式加法问题的算法实现方法,包括输入格式、输出格式及具体的算法流程。通过构造结构体存储多项式的指数和系数,实现了两个多项式的相加,并输出结果。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

This time, you are supposed to find A+B where A and B are two polynomials.

Input

Each input file contains one test case. Each case occupies 2 lines, and each line contains the information of a polynomial: K N1 aN1 N2 aN2 ... NK aNK, where K is the number of nonzero terms in the polynomial, Ni and aNi (i=1, 2, ..., K) are the exponents and coefficients, respectively. It is given that 1 <= K <= 10,0 <= NK < ... < N2 < N1 <=1000.

Output

For each test case you should output the sum of A and B in one line, with the same format as the input. Notice that there must be NO extra space at the end of each line. Please be accurate to 1 decimal place.

Sample Input
2 1 2.4 0 3.2
2 2 1.5 1 0.5
Sample Output
3 2 1.5 1 2.9 0 3.2

 算法:

  1. 算法开始。
  2. 构造结构体term类型,包含in类型的ni,double类型的ai。
  3. 读入所有数据,第一个多项式项数为n1,第二个多项式项数为n2,分别存储于term类型的ptr1和ptr2。
  4. 如果i不小于n1而且j不小于n2,那么继续,否则跳到
  5. 如果ptr[i].ni > ptr2[j].ni,ptr3[k++]=ptr1[i++],跳回第四步。
  6. 如果ptr[i].ni == ptr2[j].ni,ptr3[k].ni=ptr1[i].ni,ptr3[k].ai=ptr1[i].ai+ptr2[j].ai。否则跳到第九步。
  7. 如果fabs(ptr3[k].ai)>FLT_EPSILON,那么k++。
  8. i++, j++。跳回第四步。
  9. ptr3[k++]=ptr2[j++],跳回第四步。
  10. 如果i不小于n1,跳到第十四步。
  11. ptr3[k++]=ptr1[i++]。跳回第十步。
  12. 如果j不小于n2,跳到第十四步。
  13. ptr3[k++]=ptr2[j++]。跳回第十二步。
  14. 如果k等于零,输出零。跳到第十八步。
  15. 输出k。i=0。
  16. 如果i不小于k,则跳到第十八步。
  17. 如果fabs(ptr3[i].ai)>FLT_EPSILON,输出ptr3[i].ni和ptr3[i].ai。回到第十六步。
  18. 算法结束。
算法时间复杂度为O(n),空间复杂度O(1)。注意:如果非零项数为零,那么只要输出项数零就可以了,不需要输出等于零的那一项。
下列代码编写于VS2015,修改scanf_s()函数就可以在PAT上AC。
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <float.h>
#define MAX 100
typedef struct term {
	int ni;
	double ai;
}term;
int main(void) {
	int i, j, k, n1, n2;
	term* ptr1 = (term*)calloc(MAX, sizeof(term));
	term* ptr2 = (term*)calloc(MAX, sizeof(term));
	term* ptr3 = (term*)calloc(MAX, sizeof(term));
	scanf_s("%d", &n1);
	for (i = 0; i < n1; i++) {
		scanf_s("%d %lf", &(ptr1[i].ni), &(ptr1[i].ai));
	}
	scanf_s("%d", &n2);
	for (i = 0; i < n2; i++) {
		scanf_s("%d %lf", &(ptr2[i].ni), &(ptr2[i].ai));
	}
	for (i = 0, j = 0, k = 0; i < n1&&j < n2;) {
		if (ptr1[i].ni > ptr2[j].ni) {
			ptr3[k++] = ptr1[i++];
		}
		else if (ptr1[i].ni == ptr2[j].ni) {
			ptr3[k].ni = ptr1[i].ni;
			ptr3[k].ai = ptr1[i].ai + ptr2[j].ai;
			if (fabs(ptr3[k].ai) > FLT_EPSILON) {
				k++;
			}
			i++;
			j++;
		}
		else {
			ptr3[k++] = ptr2[j++];
		}
	}
	while (i < n1) {
		ptr3[k++] = ptr1[i++];
	}
	while (j < n2) {
		ptr3[k++] = ptr2[j++];
	}
	if (!k) {
		printf("%d", 0);
	}
	else {
		printf("%d", k);
		for (i = 0; i < k; i++) {
			if (fabs(ptr3[i].ai) > FLT_EPSILON) {
				printf(" %d %.1lf", ptr3[i].ni, ptr3[i].ai);
			}
		}
	}
	free(ptr1);
	free(ptr2);
	free(ptr3);
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值