/* -*- coding: utf-8 -*- * * 1611.cc: No.1611 Minimum Multiple with Double Divisors - yukicoder */ #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; /* constant */ const int MAX_P = 32; /* typedef */ typedef long long ll; typedef vector vi; typedef pair pii; typedef vector vpii; /* global variables */ bool primes[MAX_P + 1]; vpii pdss[MAX_P]; int pfs[MAX_P]; /* subroutines */ int gen_primes(int maxp, vi &pnums) { memset(primes, true, sizeof(primes)); primes[0] = primes[1] = false; int p; for (p = 2; p * p <= maxp; p++) if (primes[p]) { pnums.push_back(p); for (int q = p * p; q <= maxp; q += p) primes[q] = false; } for (; p <= maxp; p++) if (primes[p]) pnums.push_back(p); return (int)pnums.size(); } bool prime_decomp(int n, vi &pnums, vpii& pds) { pds.clear(); int pn = pnums.size(); for (int i = 0; i < pn; i++) { int pi = pnums[i]; if (pi * pi > n) { if (n > 1) pds.push_back(pii(n, 1)); return true; } if (n % pi == 0) { int fi = 0; while (n % pi == 0) n /= pi, fi++; pds.push_back(pii(pi, fi)); } } return false; } int pfact(ll x, int p) { int f = 0; while (x % p == 0) x /= p, f++; return f; } /* main */ int main() { vi pnums; int pn = gen_primes(MAX_P, pnums); for (int q = 2; q < MAX_P; q++) prime_decomp(q, pnums, pdss[q]); int tn; scanf("%d", &tn); while (tn--) { ll x; scanf("%lld", &x); for (int i = 0; i < pn; i++) pfs[pnums[i]] = pfact(x, pnums[i]); int q; for (q = 2; q < MAX_P; q++) { int k0 = 1, k1 = 1; for (auto &pd: pdss[q]) { int f = pfs[pd.first]; k0 *= f + 1; k1 *= f + 1 + pd.second; } if (k1 == k0 * 2) break; } printf("%lld\n", x * q); } return 0; }