1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
|
#include <stdio.h>
#include <stdlib.h>
#define MAX_DEGREE 100
typedef struct Polynomial* polynomial;
struct Polynomial
{
int coeff_array[MAX_DEGREE + 1];
int high_power;
};
int max(int x, int y)
{
return x > y ? x : y;
}
void zero_polynomial(polynomial poly)
{
int i;
for(i = 0; i <= MAX_DEGREE; i++)
poly->coeff_array[i] = 0;
poly->high_power == 0;
}
void add_polynomial(const polynomial poly1, const polynomial poly2,
polynomial poly_sum)
{
int i;
poly_sum->high_power = max(poly1->high_power, poly2->high_power);
for(i = 0; i <= poly_sum->high_power; i++) {
poly_sum->coeff_array[i] = poly1->coeff_array[i]
+ poly2->coeff_array[i];
}
}
void mult_polynomial(const polynomial poly1, const polynomial poly2,
polynomial poly_prod)
{
int i, j;
poly_prod->high_power = poly1->high_power + poly2->high_power;
for(i = 0; i <= poly1->high_power; i++) {
for(j = 0; j<= poly2->high_power; j++)
/* += used here */
poly_prod->coeff_array[i + j] += poly1->coeff_array[i]
* poly2->coeff_array[j];
}
}
void print_poly(const polynomial poly)
{
int i;
for(i = 0; i <= poly->high_power; i++)
if(poly->coeff_array[i] != 0)
printf("%dx^%d + ", poly->coeff_array[i], i);
printf("\n");
}
void test()
{
int i;
polynomial poly1, poly2, poly_sum, poly_prod;
poly1 = malloc(sizeof(struct Polynomial));
poly2 = malloc(sizeof(struct Polynomial));
poly_sum = malloc(sizeof(struct Polynomial));
poly_prod = malloc(sizeof(struct Polynomial));
for(i = 0; i <= 20; i++) {
if(i % 2 == 0) {
poly1->coeff_array[i] = i;
poly1->high_power = i;
} else {
poly2->coeff_array[i] = i;
poly2->high_power = i;
}
}
print_poly(poly1);
print_poly(poly2);
/* test sum */
add_polynomial(poly1, poly2, poly_sum);
print_poly(poly_sum);
/* test mult foo */
mult_polynomial(poly1, poly2, poly_prod);
print_poly(poly_prod);
}
int main()
{
test();
return 0;
}
|
