/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C COMPUTE MERTENS FUNCTION (using Deleglise and Rivat's algorithm) C C 12/07/15 (DKC) C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/ #include <math.h> long long newlocl(unsigned long long x, unsigned int u, char *mob, int *M); long long newhicl(unsigned long long x, unsigned int u, char *mob, int *M); int newrivl(unsigned long long x, unsigned int u, char *mobb, int *M) { int f1,f2,f3; f1=(int)newlocl(x,u,mobb,M); f2=(int)newhicl(x,u,mobb,M); f3=M[u-1]-f1-f2; return(f3); }