A093770 Non-perfect powers k for which q = A051903(k)/A051904(k) is an integer, A051904(k) > 1.

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

Amiram Eldar, Table of n, a(n) for n = 1..10000

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

Cf. A001597, A051903, A051904, A093769.

Sequence in context: A010089 A216142 A340303 * A190115 A203394 A144716

Adjacent sequences: A093767 A093768 A093769 * A093771 A093772 A093773

Keyword

nonn

Author

Labos Elemer, Apr 16 2004

Status

approved