A176031 a(n) = n^rad(n)

1, 4, 27, 16, 3125, 46656, 823543, 64, 729, 10000000000, 285311670611, 2985984, 302875106592253, 11112006825558016, 437893890380859375, 256, 827240261886336764177, 34012224, 1978419655660313589123979

Offset

1,2

Comments

rad(n) is the product of the primes dividing n (A007947).

Links

T. D. Noe, Table of n, a(n) for n = 1..100

Formula

a(n)= n^A007947(n).

Example

a(1) = 1^1 = 1 ; a(2) = 2^2 = 4 ; a(3) = 3^3 = 27 ; a(4) = 4^2 = 16.

Maple

with(numtheory):for n from 1 to 200 do:t1:= ifactors(n)[2] : p :=mul(t1[i][1], i=1..nops(t1)):x:=n^p:print(x):od:

Mathematica

Table[n^(Times @@ Transpose[FactorInteger[n]][[1]]), {n, 20}]

Crossrefs

Cf. A007947.

Sequence in context: A351365 A200768 A176147 * A093845 A351747 A250095

Adjacent sequences: A176028 A176029 A176030 * A176032 A176033 A176034

Keyword

nonn,easy

Author

Michel Lagneau, Apr 06 2010

Status

approved