Last active
October 18, 2016 04:23
-
-
Save CuyiGuaton/eb98e70676c676c4a51f5401027da30a to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#include <cstdlib> | |
#include <gmp.h> | |
#include <iostream> | |
using namespace std; | |
void Div2(mpz_t u, mpz_t v, mpz_t A, mpz_t B) { | |
mpz_t aux; | |
mpz_init(aux); | |
if (mpz_even_p(u) != 0 && mpz_even_p(v) != 0) { | |
mpz_tdiv_q_ui(u, u, 2); | |
mpz_tdiv_q_ui(v, v, 2); | |
} else { | |
//(s+B)/2 | |
mpz_add(aux, u, B); | |
mpz_tdiv_q_ui(u, aux, 2); | |
//(t-A)/2 | |
mpz_sub(aux, v, A); | |
mpz_tdiv_q_ui(v, aux, 2); | |
} | |
mpz_clear(aux); | |
} | |
int main(int argc, char const *argv[]) { | |
mpz_t a; | |
mpz_t b; | |
mpz_t d; | |
mpz_t s; | |
mpz_t t; | |
mpz_t s1; | |
mpz_t s2; | |
mpz_t t1; | |
mpz_t t2; | |
mpz_t u; | |
mpz_t v; | |
mpz_t r; | |
mpz_t A; | |
mpz_t B; | |
mpz_t aux; | |
mpz_init(a); | |
mpz_init(b); | |
mpz_init(d); | |
mpz_init(s); | |
mpz_init(t); | |
mpz_init(s1); | |
mpz_init(s2); | |
mpz_init(t1); | |
mpz_init(t2); | |
mpz_init(u); | |
mpz_init(v); | |
mpz_init(r); | |
mpz_init(A); | |
mpz_init(B); | |
mpz_init(aux); | |
mpz_set_str(a, argv[1], 10); | |
mpz_set_str(b, argv[2], 10); | |
if (argc != 3) | |
cout << "Error, debe ingresar los parametros a y b de los que sacar el " | |
"mcd(a,b)"; | |
else if (mpz_cmp_ui(a, 0) == 0 || mpz_cmp_ui(b, 0) == 0 ) | |
cout << "Error, a o b no pueden ser 0"; | |
else { | |
mpz_set_ui(d, 1); | |
// si a y b son par. | |
while (mpz_even_p(a) != 0 && mpz_even_p(b) != 0) { | |
mpz_tdiv_q_ui(a, a, 2); // a=a/2 trunca el resultado | |
mpz_tdiv_q_ui(b, b, 2); // b=b/2 | |
mpz_mul_ui(d, d, 2); // d=d*2 | |
} | |
mpz_set(A, a); | |
mpz_set(B, b); | |
mpz_set_ui(s1, 1); | |
mpz_set_ui(t1, 0); | |
mpz_set_ui(s2, 0); | |
mpz_set_ui(t2, 1); | |
// mientras a sea par | |
while (mpz_even_p(a) != 0) { | |
mpz_tdiv_q_ui(a, a, 2); // a=a/2 | |
Div2(s1, t1, A, B); | |
} | |
// mientras b sea par | |
while (mpz_even_p(b) != 0) { | |
mpz_tdiv_q_ui(b, b, 2); // b=b/2 | |
Div2(s2, t2, A, B); | |
} | |
// si b>a | |
if (mpz_cmp(b, a) > 0) { | |
mpz_set(r, b); | |
mpz_set(b, a); | |
mpz_set(a, r); | |
mpz_set(u, s2); | |
mpz_set(s2, s1); | |
mpz_set(s1, u); | |
mpz_set(v, t2); | |
mpz_set(t2, t1); | |
mpz_set(t1, v); | |
} | |
// mientras b!=0 | |
while (mpz_cmp_ui(b, 0) != 0) { | |
mpz_sub(r, a, b); // r= a-b | |
mpz_sub(u, s1, s2); // u = s1-s2 | |
mpz_sub(v, t1, t2); // v=t1-t2 | |
// mpz_out_str(stdout, 10, r); esto es para depurar los valores de r | |
// cout<<endl; | |
while (mpz_odd_p(r) == 0 && mpz_cmp_ui(r, 0) != 0) { // mientras r sea par | |
mpz_tdiv_q_ui(r, r, 2); // r=r/2 | |
Div2(u, v, A, B); | |
} | |
if (mpz_cmp(r, b) >= 0) { | |
mpz_set(a, r); | |
mpz_set(s1, u); | |
mpz_set(t1, v); | |
} else { | |
mpz_set(a, b); | |
mpz_set(b, r); | |
mpz_set(s1, s2); | |
mpz_set(s2, u); | |
mpz_set(t1, t2); | |
mpz_set(t2, v); | |
} | |
} | |
mpz_mul(d, d, a); | |
mpz_set(s, s1); | |
mpz_set(t, t1); | |
cout << endl; | |
mpz_out_str(stdout, 10, d); | |
cout << " = "; | |
mpz_out_str(stdout, 10, s); | |
cout << "*"; | |
mpz_out_str(stdout, 10, A); | |
cout << " + "; | |
mpz_out_str(stdout, 10, t); | |
cout << "*"; | |
mpz_out_str(stdout, 10, B); | |
cout << endl; | |
mpz_clear(a); | |
mpz_clear(b); | |
mpz_clear(d); | |
mpz_clear(s); | |
mpz_clear(t); | |
mpz_clear(s1); | |
mpz_clear(s2); | |
mpz_clear(t1); | |
mpz_clear(t2); | |
mpz_clear(u); | |
mpz_clear(v); | |
mpz_clear(r); | |
mpz_clear(A); | |
mpz_clear(B); | |
} | |
return 0; | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment