x,p = map(int,input().split())
a = list(map(int,input().split()))
q = p; m = max(a); ans = 0
while q<m:
    r = [0]*q
    for v in a: r[v%q] += 1
    for v in r: ans += v*(v-1)//2
    q *= p
print(ans)