Can't Copy and Paste this?

Click here for a copy-and-paste friendly version of this code!

//**************************************

//     

//INCLUDE files for :Compute Euler's Tot

//     ient function (phi)

//**************************************

//     

/* +++Date last modified: 05-Jul-1997 */

/*

** SNIPMATH.H - Header file for SNIPPETS math functions and macros

*/

#ifndef SNIPMATH__H

#define SNIPMATH__H

#include 

#include "sniptype.h"

#include "round.h"

/*

** Callable library functions begin here

*/

voidSetBCDLen(int n); /* Bcdl.C */

longBCDtoLong(char *BCDNum); /* Bcdl.C */

voidLongtoBCD(long num, char BCDNum[]);/* Bcdl.C */

double bcd_to_double(void *buf, size_t len, /* Bcdd.C */

int digits);

int double_to_bcd(double arg, char *buf, /* Bcdd.C */

size_t length, size_t digits );

DWORDncomb1 (int n, int m);/* Combin.C*/

DWORDncomb2 (int n, int m);/* Combin.C*/

voidSolveCubic(double a, double b, double c, /* Cubic.C*/

double d, int *solutions,

double *x);

DWORDdbl2ulong(double t); /* Dbl2Long.C */

longdbl2long(double t);/* Dbl2Long.C */

double dround(double x); /* Dblround.C */

/* Use #defines for Permutations and Combinations -- Factoryl.C */

#define log10P(n,r) (log10factorial(n)-log10factorial((n)-(r)))

#define log10C(n,r) (log10P((n),(r))-log10factorial(r))

double log10factorial(double N); /* Factoryl.C */

double fibo(unsigned short term);/* Fibo.C */

double frandom(int n);/* Frand.C*/

double ipow(double x, int n);/* Ipow.C */

int ispow2(int x);/* Ispow2.C*/

longdouble ldfloor(long double a);/* Ldfloor.C */

int initlogscale(long dmax, long rmax);/* Logscale.C */

longlogscale(long d); /* Logscale.C */

floatMSBINToIEEE(float f); /* Msb2Ieee.C */

floatIEEEToMSBIN(float f); /* Msb2Ieee.C */

int perm_index (char pit[], int size);/* Perm_Idx.C */

int round_div(int n, int d); /* Rnd_Div.C */

longround_ldiv(long n, long d);/* Rnd_Div.C */

double rad2deg(double rad); /* Rad2Deg.C */

double deg2rad(double deg); /* Rad2Deg.C */

#include "pi.h"

#ifndef PHI

#define PHI ((1.0+sqrt(5.0))/2.0) /* the golden number*/

#define INV_PHI (1.0/PHI) /* the golden ratio */

#endif

/*

** File: ISQRT.C

*/

    struct int_sqrt {

    unsigned sqrt,

    frac;

};

void usqrt(unsigned long x, struct int_sqrt *q);

#endif /* SNIPMATH__H */

code:

Can't Copy and Paste this?

Click here for a copy-and-paste friendly version of this code!

Terms of Agreement:   

By using this code, you agree to the following terms...   

1) You may use this code in your own programs (and may compile it into a program and distribute it in compiled format for langauges that allow it) freely and with no charge.   

2) You MAY NOT redistribute this code (for example to a web site) without written permission from the original author. Failure to do so is a violation of copyright laws.   

3) You may link to this code from another website, but ONLY if it is not wrapped in a frame. 

4) You will abide by any additional copyright restrictions which the author may have placed in the code or code's description.  

//     

/* +++Date last modified: 05-Jul-1997 */

/*

** ETPHI.c

**

** Compute Euler's Totient function, phi().

** phi(n) = the number of positive integers less than n

** which have no common factors with n

**

** Written by M. Stapleton of Graphic Bits

** public Domain

**

** Sep 7 1996

**

** Uses ISQRT.C, also in SNIPPETS

*/

#include "snipmath.h"

#undef ULONG

typedef unsigned long ULONG;

ULONG intsqrt(ULONG n);

static void dofact(ULONG i);

ULONG etphi(ULONG n);

static ULONG num, ph;

void dofact(ULONG i)

    {

    long p;

    for (p = 0; (num % i) == 0; p++)

    num /= i;

    if (p)

        {

        ph *= i - 1;

        while(--p)

        ph *= i;

    }

}

ULONG etphi(ULONG n)

    {

    ULONG i;

    struct int_sqrt mi;

    if (n<2)

    return 0;

    num = n;

    ph = 1;

    usqrt(n, &mi);

    dofact(2);

    dofact(3);

    for (i = 5; (i <= mi.sqrt) && (num >= i); i += 6)

        {

        dofact(i);

        if (num < i + 2)

        break;

        dofact(i + 2);

    }

    if (num > 1)

    ph *= num - 1;

    return ph;

}

#ifdef TEST

#include 

#include 

int main(int argc, char *argv[])

    {

    long n = 0;

    if (argc > 1)

    n = atol(argv[1]);

    printf("etphi(%ld) = %ld\n", n, etphi(n));

    return 0;

}

#endif