# calculate φ(x) for 1 <= x <= M
M = 10**7+10
*phi, = range(M+1)
for x in range(2, M+1):
    if phi[x] == x:
        for y in range(x, M+1, x):
            phi[y] = phi[y] // x * (x-1)

for i in range(1,M+1):
    phi[i] += phi[i-1]
def sol():
    n = int(input())
    ans = (-phi[n] + 1) + n*(n-1)
    print(ans)

t = int(input())
for i in range(t):
    sol()