#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <stdint.h>
#include <limits.h>
#include <float.h>

#define _CRT_SECURE_NO_WARNINGS
#define LN10    2.302585092994045684018
#define RLN10  0.43429448190325182765

#define  MAX_LEN  128
......................

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <assert.h>

#define MAX_Q  1000000000 
// 返回最大下标
int calc_cfrac(double x, int a[], int p[], int q[], int max_idx)
{
   double ipart;
   int i;
   
......................

#include <stdio.h>
#include <stdlib.h>
#include <float.h>
#include <stdint.h>
#include <math.h>
#include <mpfr.h>
#include <mpreal.h>
#include <time.h>
#include <inttypes.h>

#define LG2  0.30102999566398119521
#define RLN2 1.44269504088896340735
......................

#include <stdlib.h>
#include <stdio.h>
#include <stdint.h>

typedef union
{
   double f;
   uint32_t i[2];
}Packet;

#define M  ((1<<20)-1)
void test( double f)
......................

#include <stdio.h>
#include <stdlib.h>
#include <assert.h>
#include <math.h>
#include <limits.h>

int primes[] =
{
    2, 3, 5, 7, 11, 13, 17, 19,
    23, 29, 31, 37, 41, 43, 47, 53,
    59, 61, 67, 71, 73, 79, 83, 89,
    97, 101, 103, 107, 109, 113, 127, 131,
......................