fork download

copy

#include <gmp.h>
#include <stdio.h>
 
//
// 一般項
//
mp_bitcnt_t gmp_fibonacci_f(unsigned int n_, mpz_t y_, mp_bitcnt_t prec_) {
 
	mpf_t r5;
	mpf_t ltmp;
	mpf_t rtmp;
	mp_bitcnt_t prec;
 
	mpf_init2(r5, prec_);
	mpf_init2(ltmp, prec_);
	mpf_init2(rtmp, prec_);
 
	mpf_sqrt_ui(r5, 5);
 
	mpf_add_ui(rtmp, r5, 1);
	mpf_div_ui(rtmp, rtmp, 2);
	mpf_pow_ui(ltmp, rtmp, n_);
	mpf_ui_sub(rtmp, 1, rtmp);
	mpf_pow_ui(rtmp, rtmp, n_);
	mpf_sub(ltmp, ltmp, rtmp);
	mpf_div(ltmp, ltmp, r5);
	mpf_set_d(rtmp, 0.5);
	mpf_add(ltmp, ltmp, rtmp);
	mpz_set_f(y_, ltmp);
 
	prec = mpf_get_prec(r5);
 
	mpf_clear(ltmp);
	mpf_clear(rtmp);
	mpf_clear(r5);
 
	return prec;
}
 
mp_bitcnt_t gmp_zgetbits(mpz_t y_) {
        mp_bitcnt_t prev_bits = 0;
        mp_bitcnt_t bits;
 
        bits = mpz_scan1(y_, prev_bits);
        if (bits == -1) {
                return 1;
        }
        prev_bits = 1;
 
        for (;;) {
                bits = mpz_scan1(y_, prev_bits);
                if (bits == -1) {
                        break;
                }
                prev_bits = bits + 1;
        }
 
        return prev_bits;
}
 
//
// 必要な精度を試しに取得
//
mp_bitcnt_t gmp_try_fibonacci_f(unsigned int n_, mpz_t y_) {
 
	mp_bitcnt_t prec = 1;
	mpz_t try_y;
	mpz_init(try_y);
 
	for ( ; ; ++prec) {
		prec = gmp_fibonacci_f(n_, try_y, prec);
 
		if (!mpz_cmp(try_y, y_)) {
			break;
		}
	}
 
	mpz_clear(try_y);
	return prec;
}
 
int main(int argc, char* argv[]) {
 
	mpz_t gy;
	mpz_t fy;
	unsigned int res_bits;
	unsigned int req_bits;
	unsigned int est_res_bits;
 
	mpz_init(gy);
	mpz_init(fy);
 
	fprintf(stdout, "n,res,req,est_res\n");
	for (unsigned int i = 0; i < 10000; ++i) {
 
		mpz_fib_ui(gy, i);
		res_bits = gmp_zgetbits(gy);
		gmp_fibonacci_f(i, fy, 32);
		est_res_bits = gmp_zgetbits(fy);
		req_bits     = gmp_try_fibonacci_f(i, gy);
		fprintf(stdout, "%d,%d,%d,%d\n", i, res_bits, req_bits, est_res_bits);
	}
 
	mpz_clear(gy);
	mpz_clear(fy);
	return 0;
}
 
 

/home/TYxDZS/ccjoN1iO.o: In function `gmp_fibonacci_f':
prog.c:(.text+0x28): undefined reference to `__gmpf_init2'
prog.c:(.text+0x33): undefined reference to `__gmpf_init2'
prog.c:(.text+0x3e): undefined reference to `__gmpf_init2'
prog.c:(.text+0x4b): undefined reference to `__gmpf_sqrt_ui'
prog.c:(.text+0x5b): undefined reference to `__gmpf_add_ui'
prog.c:(.text+0x6b): undefined reference to `__gmpf_div_ui'
prog.c:(.text+0x79): undefined reference to `__gmpf_pow_ui'
prog.c:(.text+0x89): undefined reference to `__gmpf_ui_sub'
prog.c:(.text+0x97): undefined reference to `__gmpf_pow_ui'
prog.c:(.text+0xa5): undefined reference to `__gmpf_sub'
prog.c:(.text+0xb3): undefined reference to `__gmpf_div'
prog.c:(.text+0xc3): undefined reference to `__gmpf_set_d'
prog.c:(.text+0xd1): undefined reference to `__gmpf_add'
prog.c:(.text+0xdc): undefined reference to `__gmpz_set_f'
prog.c:(.text+0xe4): undefined reference to `__gmpf_get_prec'
prog.c:(.text+0xef): undefined reference to `__gmpf_clear'
prog.c:(.text+0xf7): undefined reference to `__gmpf_clear'
prog.c:(.text+0xff): undefined reference to `__gmpf_clear'
/home/TYxDZS/ccjoN1iO.o: In function `gmp_zgetbits':
prog.c:(.text+0x131): undefined reference to `__gmpz_scan1'
prog.c:(.text+0x14b): undefined reference to `__gmpz_scan1'
/home/TYxDZS/ccjoN1iO.o: In function `gmp_try_fibonacci_f':
prog.c:(.text+0x177): undefined reference to `__gmpz_init'
prog.c:(.text+0x1a1): undefined reference to `__gmpz_cmp'
prog.c:(.text+0x1ad): undefined reference to `__gmpz_clear'
/home/TYxDZS/ccjoN1iO.o: In function `main':
prog.c:(.text.startup+0x23): undefined reference to `__gmpz_init'
prog.c:(.text.startup+0x2b): undefined reference to `__gmpz_init'
prog.c:(.text.startup+0x57): undefined reference to `__gmpz_fib_ui'
prog.c:(.text.startup+0xb7): undefined reference to `__gmpz_clear'
prog.c:(.text.startup+0xbf): undefined reference to `__gmpz_clear'
collect2: error: ld returned 1 exit status

Standard output is empty