diff options
author | Torbjorn Granlund <tg@gmplib.org> | 2019-09-15 14:47:43 +0200 |
---|---|---|
committer | Torbjorn Granlund <tg@gmplib.org> | 2019-09-15 14:47:43 +0200 |
commit | bc2774c8baa5b4aad105e7c49f49076035cc9a22 (patch) | |
tree | 302073c7f9c23cc7a92b6f60932f9e5ddf1cd78f /mpz/bin_ui.c | |
parent | 0ce395bd6fba963e927960865526b617a7be5270 (diff) | |
download | gmp-bc2774c8baa5b4aad105e7c49f49076035cc9a22.tar.gz |
Whitespace cleanup.
Diffstat (limited to 'mpz/bin_ui.c')
-rw-r--r-- | mpz/bin_ui.c | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/mpz/bin_ui.c b/mpz/bin_ui.c index 753c2b102..bab3e0769 100644 --- a/mpz/bin_ui.c +++ b/mpz/bin_ui.c @@ -175,25 +175,25 @@ rek_raising_fac4 (mpz_ptr r, mpz_ptr p, mpz_ptr P, unsigned long int k, unsigned Assume k = 4i then compute: p = (n+1)(n+4i)/2 - i - (n+1+1)(n+4i)/2 = p + i + (n+4i)/2 - (n+1+1)(n+4i-1)/2 = p + i + ((n+4i)-(n+1+1))/2 = p + i + (n-n+4i-2)/2 = p + 3i-1 + (n+1+1)(n+4i)/2 = p + i + (n+4i)/2 + (n+1+1)(n+4i-1)/2 = p + i + ((n+4i)-(n+1+1))/2 = p + i + (n-n+4i-2)/2 = p + 3i-1 P = (p + i)*(p+3i-1)/2 = (n+1)(n+2)(n+4i-1)(n+4i)/8 n' = n + 2 i' = i - 1 - (n'-1)(n')(n'+4i'+1)(n'+4i'+2)/8 = P + (n'-1)(n')(n'+4i'+1)(n'+4i'+2)/8 = P (n'-1)(n'+4i'+2)/2 - i' - 1 = p (n'-1+2)(n'+4i'+2)/2 - i' - 1 = p + (n'+4i'+2) (n'-1+2)(n'+4i'+2-2)/2 - i' - 1 = p + (n'+4i'+2) - (n'-1+2) = p + 4i' + 1 (n'-1+2)(n'+4i'+2-2)/2 - i' = p + 4i' + 2 p' = p + 4i' + 2 = (n'+1)(n'+4i')/2 - i' - p' - 4i' - 2 = p - (p' - 4i' - 2 + i)*(p' - 4i' - 2+3i-1)/2 = P - (p' - 4i' - 2 + i' + 1)*(p' - 4i' - 2 + 3i' + 3 - 1)/2 = P - (p' - 3i' - 1)*(p' - i')/2 = P - (p' - 3i' - 1 + 4i' + 1)*(p' - i' + 4i' - 1)/2 = P + (4i' + 1)*(p' - i')/2 + (p' - 3i' - 1 + 4i' + 1)*(4i' - 1)/2 - (p' + i')*(p' + 3i' - 1)/2 = P + (4i')*(p' + p')/2 + (p' - i' - (p' + i'))/2 - (p' + i')*(p' + 3i' - 1)/2 = P + 4i'p' + (p' - i' - p' - i')/2 - (p' + i')*(p' + 3i' - 1)/2 = P + 4i'p' - i' + p' - 4i' - 2 = p + (p' - 4i' - 2 + i)*(p' - 4i' - 2+3i-1)/2 = P + (p' - 4i' - 2 + i' + 1)*(p' - 4i' - 2 + 3i' + 3 - 1)/2 = P + (p' - 3i' - 1)*(p' - i')/2 = P + (p' - 3i' - 1 + 4i' + 1)*(p' - i' + 4i' - 1)/2 = P + (4i' + 1)*(p' - i')/2 + (p' - 3i' - 1 + 4i' + 1)*(4i' - 1)/2 + (p' + i')*(p' + 3i' - 1)/2 = P + (4i')*(p' + p')/2 + (p' - i' - (p' + i'))/2 + (p' + i')*(p' + 3i' - 1)/2 = P + 4i'p' + (p' - i' - p' - i')/2 + (p' + i')*(p' + 3i' - 1)/2 = P + 4i'p' - i' P' = P + 4i'p' - i' And compute the product P * P' * P" ... @@ -417,7 +417,7 @@ mpz_bin_ui (mpz_ptr r, mpz_srcptr n, unsigned long int k) } #endif else - { /* k = 1 */ + { /* k = 1 */ mpz_add_ui (r, ni, 1); } } |