/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C                                                                             C
C  COMPUTE NUMBER OF FRACTIONS						      C
C  12/19/14 (DKC)							      C
C                                                                             C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/
#include <math.h>
#include <stdlib.h>
unsigned int haros5(unsigned int N, unsigned int M, unsigned int *R,
		    unsigned int H, unsigned int K, unsigned int HP,
		    unsigned int KP, unsigned int D);
unsigned int lagrange1(unsigned int N, unsigned int D, unsigned int O);
void mertensi(unsigned int N, unsigned int *count, unsigned int *S) {
unsigned int L,count1,count2,ctemp1,ctemp2;
if (N<7) {
   count[0]=0;
   count[1]=0;
   return;
   }
if (N==7) {
   count[0]=0;
   count[1]=1;
   return;
   }
if (N==8) {
   count[0]=0;
   count[1]=2;
   return;
   }
if (N==9) {
   count[0]=0;
   count[1]=3;
   return;
   }
if (N==10) {
   count[0]=0;
   count[1]=4;
   return;
   }
if (N==11) {
   count[0]=0;
   count[1]=5;
   return;
   }
if (N==12) {
   count[0]=0;
   count[1]=5;
   return;
   }
if ((N>12)&&(N<129)) {
   count1=haros5(N,0,S,0,1,1,N,12);
   count1=count1-2;
   count[0]=count1;
   L=lagrange1(1,12,N);
   count2=haros5(N,0,S,1,12,L>>16,L&0xffff,6);
   count2=count2-2;
   count[1]=count2;
   return;
   }
if (N>50000) {
   count[0]=0;
   count[1]=0;
   return;
   }
ctemp1=haros5(N,0,S,0,1,1,N,128);	      // towards 1/128
L=lagrange1(1,128,N);
count1=haros5(N,0,S,1,128,L>>16,L&0xffff,64); // towards 1/64
count1=count1+ctemp1-2;

L=lagrange1(1,64,N);
ctemp1=haros5(N,0,S,1,64,L>>16,L&0xffff,128);  // towards 3/128
count1=count1+ctemp1-1;
L=lagrange1(3,128,N);
ctemp1=haros5(N,0,S,3,128,L>>16,L&0xffff,32);  // towards 1/32
count1=count1+ctemp1-1;

L=lagrange1(1,32,N);
ctemp1=haros5(N,0,S,1,32,L>>16,L&0xffff,128);  // towards 5/128
count1=count1+ctemp1-1;
L=lagrange1(5,128,N);
ctemp1=haros5(N,0,S,5,128,L>>16,L&0xffff,64);  // towards 3/64
count1=count1+ctemp1-1;

L=lagrange1(3,64,N);
ctemp1=haros5(N,0,S,3,64,L>>16,L&0xffff,128);  // towards 7/128
count1=count1+ctemp1-1;
L=lagrange1(7,128,N);
ctemp1=haros5(N,0,S,7,128,L>>16,L&0xffff,16);  // towards 1/16
count1=count1+ctemp1-1;

L=lagrange1(1,16,N);
ctemp1=haros5(N,0,S,1,16,L>>16,L&0xffff,128);  // towards 9/128
count1=count1+ctemp1-1;
L=lagrange1(9,128,N);
ctemp1=haros5(N,0,S,9,128,L>>16,L&0xffff,12);  // towards 1/12
count1=count1+ctemp1-2;
count[0]=count1;
//
//
//
L=lagrange1(1,12,N);
ctemp2=haros5(N,0,S,1,12,L>>16,L&0xffff,128);  // towards 11/128
L=lagrange1(11,128,N);
count2=haros5(N,0,S,11,128,L>>16,L&0xffff,32);	// towards 3/32
count2=count2+ctemp2-1;

L=lagrange1(3,32,N);
ctemp2=haros5(N,0,S,3,32,L>>16,L&0xffff,128);  // towards 13/128
count2=count2+ctemp2-1;
L=lagrange1(13,128,N);
ctemp2=haros5(N,0,S,13,128,L>>16,L&0xffff,64);	 // towards 7/64
count2=count2+ctemp2-1;

L=lagrange1(7,64,N);
ctemp2=haros5(N,0,S,7,64,L>>16,L&0xffff,128);	// towards 15/128
count2=count2+ctemp2-1;
L=lagrange1(15,128,N);
ctemp2=haros5(N,0,S,15,128,L>>16,L&0xffff,8);	// towards 1/8
count2=count2+ctemp2-1;

L=lagrange1(1,8,N);
ctemp2=haros5(N,0,S,1,8,L>>16,L&0xffff,128);   // towards 17/128
count2=count2+ctemp2-1;
L=lagrange1(17,128,N);
ctemp2=haros5(N,0,S,17,128,L>>16,L&0xffff,64);	// towards 9/64
count2=count2+ctemp2-1;

L=lagrange1(9,64,N);
ctemp2=haros5(N,0,S,9,64,L>>16,L&0xffff,128);  // towards 19/128
count2=count2+ctemp2-1;
L=lagrange1(19,128,N);
ctemp2=haros5(N,0,S,19,128,L>>16,L&0xffff,32);	// towards 5/32
count2=count2+ctemp2-1;

L=lagrange1(5,32,N);
ctemp2=haros5(N,0,S,5,32,L>>16,L&0xffff,128);  // towards 21/128
count2=count2+ctemp2-1;
L=lagrange1(21,128,N);
ctemp2=haros5(N,0,S,21,128,L>>16,L&0xffff,6);  // towards 1/6
count2=count2+ctemp2-3;
count[1]=count2;
return;
}