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

Author

Michel Lagneau, Apr 06 2010

Keywords

• nonn
• easy

Comments

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

Examples

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

Links

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

Crossrefs

Cf. A007947.

Programs

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}]

Formula

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