import sys
input = lambda: sys.stdin.readline().rstrip()
ii = lambda: int(input())
mi = lambda: map(int, input().split())
li = lambda: list(mi())
INF = 2**63-1
mod = 998244353


# calculate φ(x) for 1 <= x <= M
M = 10**7 + 1
phi = list(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 x in range(2, M+1):
    phi[x] += phi[x - 1] 

def solve():
    #i, jについて，互いに素なら答えは1
    #互いに素でないなら，1では行けない．i → 1 → jと行くと，2が達成できる
    #よって，答えは n * (n - 1) - Σφ(x)
    n = ii()

    print(n * (n - 1) - phi[n] + 1)
    
for _ in range(ii()):
    solve()