#include <iostream.h>
#include <time.h>
#include <stdlib.h>

/* for some reason this wasn't defined */

#define min(a,b) ((a) < (b) ? (a) : (b))

unsigned S2(unsigned x, unsigned yin);

unsigned JSHpi(unsigned xin, unsigned yin)
{
//cache most common answers
#define cache_size 32
  static unsigned f[cache_size]=
    {0,1,2,2,3,3,4,4,4,4,5,5,6,6,6,6,7,7,8,8,8,8,9,9,9,9,9,9,10,10,11,11};
  if(xin==yin && (xin<cache_size))return f[xin-1];
  return (xin-S2(xin,yin)-1);
}

unsigned S2(unsigned x, unsigned yin)
{
  unsigned sum=0, i, sum1, sum2;
  const unsigned dx=1;

  for(unsigned i=2; i*i<=yin; i++)
  {
    sum2 = ( JSHpi(i,i) - JSHpi(i-dx,i-dx));
    if (sum2!=0){
      int tmp = (i-dx)*(i-dx),
      tmp2=x/i;
      sum1 = ( JSHpi(x/i,min(tmp,tmp2)) -
          JSHpi(i-dx,i-dx));
      sum+=(sum1 * sum2);
    }
  } return sum;
}

unsigned long long pi(unsigned long long x) { return JSHpi(x, x); }