/*
* GCD -- compute the greatest common divisor of two integers recursively
*
* Usage: gcd
*
* Inputs: Integers m, n
* Outputs: gcd(m, n) = g -- m, n as above; g is gcd(m, n)
* Errors: - m, n not an integer
*
* Assumptions: - the arguments are not digits followed by non-digits
* (unchecked)
* - at least 1 of m, n is non-zero (checked)
*
* Matt Bishop, ECS 36A
* October 17, 2019 -- wrote original version
*/
#include <stdio.h>
/*
* GCD -- compute the GCD of m and n recursively
* Parameters: m, n 2 integers of which at least 1 is non-zero
* Returns: gcd of m and n
* Assumptions: - at least 1 of m, n is non-zero (not checked)
* - neither m, n are negative (not checked)
* Algorithm: Euclid's Algorithm, using the relation:
* gcd(m, n) = gcd(n, n - m)
*
*/
int gcd(int m, int n)
{
/* base case */
if (n == 0) return(m);
/* recursive case */
return(gcd(n, m % n));
}
/*
* main ro
*/
int main(void)
{
int rv; /* success of failure of read */
int ch; /* character to discard */
int m, n; /* numbers to compute the gcd of */
int badct = 0; /* number of bad arguments */
/*
* loop through the argument list
*/
printf("> ");
while((rv = scanf("%d,%d", &m, &n)) != EOF){
if (rv != 2){
fprintf(stderr, "Invalid input\n");
/* discard the rest of the line */
while((ch = getchar()) != '\n' && ch != EOF)
;
/* quit on EOF; loop on error */
if (ch == EOF)
return(badct);
badct++;
printf("> ");
continue;
}
/*
* now check error conditions and other things
*/
/* if both m, n are 0, the gcd is undefined */
if (m == 0 && n == 0){
fprintf(stderr, "need at least 1 non-zero argument\n");
printf("> ");
continue;
}
/* handle negative numbers; gcd of negatives is same */
/* as gcd of positives with same absolute values */
if (n < 0) n = -n;
if (m < 0) m = -m;
/*
* all's right with the world
* compute and print the GCD
*/
printf("gcd(%d, %d) = %d\n", m, n, gcd(m, n));
/* reprompt */
printf("> ");
}
/*
* return exit code of the number of failures
*/
return(badct);
}
|
ECS 36A, Programming & Problem Solving Version of April 2, 2024 at 12:13PM
|
You can get the raw source code here. |