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 Input2 1 2.4 0 3.2 2 2 1.5 1 0.5Sample Output
3 2 1.5 1 2.9 0 3.2
算法:
- 算法开始。
- 构造结构体term类型,包含in类型的ni,double类型的ai。
- 读入所有数据,第一个多项式项数为n1,第二个多项式项数为n2,分别存储于term类型的ptr1和ptr2。
- 如果i不小于n1而且j不小于n2,那么继续,否则跳到
- 如果ptr[i].ni > ptr2[j].ni,ptr3[k++]=ptr1[i++],跳回第四步。
- 如果ptr[i].ni == ptr2[j].ni,ptr3[k].ni=ptr1[i].ni,ptr3[k].ai=ptr1[i].ai+ptr2[j].ai。否则跳到第九步。
- 如果fabs(ptr3[k].ai)>FLT_EPSILON,那么k++。
- i++, j++。跳回第四步。
- ptr3[k++]=ptr2[j++],跳回第四步。
- 如果i不小于n1,跳到第十四步。
- ptr3[k++]=ptr1[i++]。跳回第十步。
- 如果j不小于n2,跳到第十四步。
- ptr3[k++]=ptr2[j++]。跳回第十二步。
- 如果k等于零,输出零。跳到第十八步。
- 输出k。i=0。
- 如果i不小于k,则跳到第十八步。
- 如果fabs(ptr3[i].ai)>FLT_EPSILON,输出ptr3[i].ni和ptr3[i].ai。回到第十六步。
- 算法结束。
#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;
}