local mce, mfl, msq, mmi, mma, mab = math.ceil, math.floor, math.sqrt, math.min, math.max, math.abs local function getprimes(x) local primes = {} local allnums = {} for i = 1, x do allnums[i] = true end for i = 2, x do if allnums[i] then table.insert(primes, i) local lim = mfl(x / i) for j = 2, lim do allnums[j * i] = false end end end return primes end local function getdivisorparts(x, primes) local prime_num = #primes local tmp = {} local lim = mce(msq(x)) local primepos = 1 local dv = primes[primepos] while primepos <= prime_num and dv <= lim do if x % dv == 0 then x = mfl(x / dv) local cnt = 1 while x % dv == 0 do x = mfl(x / dv) cnt = cnt + 1 end tmp[dv] = cnt lim = mce(msq(x)) end if primepos == prime_num then break end primepos = primepos + 1 dv = primes[primepos] end if x ~= 1 then tmp[x] = 1 end return tmp end local primes = getprimes(mce(msq(31 * 100000000000))) local predvp = {0} for i = 2, 31 do predvp[i] = getdivisorparts(i, primes) end local q = io.read("*n") for iq = 1, q do local x = io.read("*n") local v = false for i = 1, #primes do local p = primes[i] if x % p ~= 0 then v = p break end end local ans = 1 * v local dvp = getdivisorparts(x, primes) local tot = 1 for k, v in pairs(dvp) do tot = tot * (1 + v) end for j = 2, ans - 1 do local dvp2 = j <= 31 and predvp[j] or getdivisorparts(j, primes) local tot2 = tot for k, v in pairs(dvp2) do if dvp[k] then tot2 = mfl(tot2 * (dvp[k] + v + 1) / (dvp[k] + 1)) else tot2 = tot2 * (v + 1) end end if tot * 2 == tot2 then ans = j break end end print(ans * x) end -- print(os.clock())