Given a rational number n
n >= 0, with denominator strictly positive
as a string (example: “2/3” in Ruby, Python, Clojure, JS, CS, Go)
or as two strings (example: “2” “3” in Haskell, Java, CSharp, C++, Swift, Dart)
or as a rational or decimal number (example: 3/4, 0.67 in R)
or two integers (Fortran)
decompose this number as a sum of rationals with numerators equal to one and without repetitions (2/3 = 1/2 + 1/6 is correct but not 2/3 = 1/3 + 1/3, 1/3 is repeated).
The algorithm must be “greedy”, so at each stage the new rational obtained in the decomposition must have a denominator as small as possible. In this manner the sum of a few fractions in the decomposition gives a rather good approximation of the rational to decompose.
2/3 = 1/3 + 1/3 doesn’t fit because of the repetition but also because the first 1/3 has a denominator bigger than the one in 1/2 in the decomposition 2/3 = 1/2 + 1/6.
Example:
(You can see other examples in “Sample Tests:”)
decompose(“21/23”) or “21” “23” or 21/23 should return
[“1/2”, “1/3”, “1/13”, “1/359”, “1/644046”] in Ruby, Python, Clojure, JS, CS, Haskell, Go
“[1/2, 1/3, 1/13, 1/359, 1/644046]” in Java, CSharp, C++, Dart
“1/2,1/3,1/13,1/359,1/644046” in C, Swift, R
Notes
The decomposition of 21/23 as
21/23 = 1/2 + 1/3 + 1/13 + 1/598 + 1/897
is exact but don’t fulfill our requirement because 598 is bigger than 359. Same for
21/23 = 1/2 + 1/3 + 1/23 + 1/46 + 1/69 (23 is bigger than 13)
or
21/23 = 1/2 + 1/3 + 1/15 + 1/110 + 1/253 (15 is bigger than 13).
The rational given to decompose could be greater than one or equal to one, in which case the first “fraction” will be an integer (with an implicit denominator of 1).
If the numerator parses to zero the function “decompose” returns [] (or “”, or “[]”).
The number could also be a decimal which can be expressed as a rational.
Example:
0.6 in Ruby, Python, Clojure,JS, CS, Julia, Go
“66” “100” in Haskell, Java, CSharp, C++, C, Swift, Scala, Kotlin, Dart, …
0.67 in R.
Ref: http://en.wikipedia.org/wiki/Egyptian_fraction
#include <stddef.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
char *series_sum(const size_t n)
{
double a=1, b=1, s = 0;
char* st = (char *)malloc( sizeof(char) * 6 );
if( n == 0 )
{
strcpy( st, "0.00" );
return st;
}
for( size_t i = 0; i < n; i++, b += 3 )
s += a / b;
int sum = (int)(s * 1000);
if( sum % 10 >= 5 )
{
sum /= 10;
sum++;
}
else
sum /= 10;
st[0] = sum / 100 + '0';
st[1] = '.';
st[2] = sum / 10 % 10 + '0';
st[3] = sum % 10 + '0';
st[4] = '\0';
return st;
}