sol1.c File Reference

Problem 26 solution More...

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>

Include dependency graph for sol1.c:

Macros
#define	MAX_DENO   2000
	limit of unit fractions
#define	MAX_LEN    (MAX_DENO + 10)
	length of resulting recurring fraction number

Functions
int	compare (const void *a, const void *b)
	comparison function for use with internal qsort algorithm More...
int	main (int argc, char *argv[])
	Main function. More...

Detailed Description

Problem 26 solution

Author

Krishna Vedala

Function Documentation

compare()

int compare	(	const void *	a,
		const void *	b
	)

comparison function for use with internal qsort algorithm

{
    return (*(unsigned short *)a - *(unsigned short *)b);
}

main()

int main	(	int	argc,
		char *	argv[]
	)

Main function.

{
    unsigned short max_digits = 0, max_idx_number = 0;
 
    clock_t start_time = clock();
    short deno;
#ifdef _OPENMP
#pragma omp for
#endif
    for (deno = 2; deno < MAX_DENO; deno++)
    {
        unsigned short remainders[MAX_LEN];
        unsigned short rem = 1, *rem_ptr = remainders;
        memset(remainders, (unsigned short)-1,
               MAX_LEN * sizeof(unsigned short));
        // remainders[0] = 1;
        // printf("1/%-4u\t ", deno);
        unsigned short index = 0, num_digits;
 
        while (rem != 0)
        {
            rem = (rem * 10) % deno;
            if (rem == 0)
            {
                index = 0;
                break;
            }
            rem_ptr = (unsigned short *)bsearch(
                &rem, remainders, MAX_LEN, sizeof(unsigned short), compare);
            // printf("%2d, ", rem);
            // printf("(%14p), ", rem_ptr);
            if (rem_ptr != NULL)
                break;
            remainders[index] = rem;
            rem_ptr = remainders;
            index++;
        }
 
        num_digits = index - (rem_ptr - remainders);
        // printf("\n\t(%14p, %14p, %4u, %4u)\n", rem_ptr, remainders, index,
        // num_digits);
#ifdef _OPENMP
#pragma omp critical
        {
#endif
            if (num_digits > max_digits)
            {
                max_digits = num_digits;
                max_idx_number = deno;
                // printf("\t (%u, %u)\n ", max_digits, max_idx_number);
            }
#ifdef _OPENMP
        }
#endif
    }
    clock_t end_time = clock();
 
    printf("Time taken: %.4g ms\n",
           1e3 * (double)(end_time - start_time) / CLOCKS_PER_SEC);
    printf("Maximum digits: %hu\t Denominator: %hu\n", max_digits,
           max_idx_number);
 
    return 0;
}

Here is the call graph for this function: