login
A093770
Non-perfect powers k for which q = A051903(k)/A051904(k) is an integer, A051904(k) > 1.
2
10800, 16200, 18000, 21168, 31752, 40500, 43200, 45000, 49392, 52272, 67500, 72000, 73008, 78408, 84672, 98000, 109512, 111132, 124848, 137200, 145800, 155952, 172800, 172872, 187272, 191664, 197568, 209088, 228528, 233928, 242000
OFFSET
1,1
LINKS
FORMULA
Solutions to integer values of q = A051903(x)/A051904(x), when A051904(x) > 1.
MATHEMATICA
ffi[x_] :=Flatten[FactorInteger[x]] ep[x_] :=Table[Part[ffi[x], 2*w], {w, 1, lf[x]}] lf[x_] :=Length[FactorInteger[x]] Do[s1=Min[ep[n]]; s2=Max[ep[n]]; If[ !Equal[s1, 1]&&IntegerQ[q=(s2/s1)]&& Equal[Union[Table[IntegerQ[Part[ep[n], j]/s1], {j, 1, lf[n]}]], {False, True}], Print[n]], {n, 2, 10000000}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Apr 16 2004
STATUS
approved