/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C                                                                             C
C  COMPUTE SUM OF EULER'S PHI FUNCTION                                        C
C  05/22/14 (DKC)							      C
C                                                                             C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/
#include <math.h>
unsigned int euclid(unsigned int d, unsigned int e);
void uplolim(unsigned int g, unsigned int *sums) {
unsigned int h,j,k,l,sum,losum,hisum;
losum=sums[0];
hisum=sums[1];
sum=1;
for (j=2; j<g; j++) {
   k=euclid(g,j);
   if (k==1)
      sum=sum+1;
   }
l=sum/4;
h=l;
if (sum!=l*4)
   h=h+1;
losum=losum+l;
hisum=hisum+h;
sums[0]=losum;
sums[1]=hisum;
return;
}