n,p=map(int,input().split())
a=list(map(int,input().split()))
a.sort()
c=[{} for i in range(40)]
L=10**10
g=0
for v in a:
  for i in range(1,40):
    if p**i<L:
      y=v%(p**i)
      if y not in c[i]:
        c[i][y]=0
      g+=c[i][y]
      c[i][y]+=1
    else:
      break
print(g)