#include<stdio.h>
#include<math.h>

long long n, b, f, nr, ex, min, fminim,contor;

int prim(long long x)
{
  if (x==1 || x==0) return 0;
  if (x==2) return 1;
  for (long long i=3; i<=x/2; i+=2)
    if (x%i==0) return 0;
  return 1;
}

void descomp(long long x)
{
  long long y=x, k; min=1000;
  fminim=0;
  if (y%2==0)
    fminim=2;
  for (long long i=3; i<=y; i+=2)
    {
      if (y%i==0)
	if (prim(i))
	  fminim=i;
	  
    }
}


int main()
{
  long long i,z,k,p;
  freopen("zero2.in","r",stdin);
  freopen("zero2.out","w",stdout);
  long long nr;
  for (z=1; z<=10; z++)
    {
      scanf("%lld %lld",&n, &b);
      descomp(b);
      fminim=b;
      k=floor(n/fminim)-1;p=1;	
      contor+=(k*(k+1)/2*pow(fminim,p) + (k+1)*(n-(k+1)*pow(fminim,p)+1));
      printf("%lld\n",contor);
    }
  return 0;
}