|
Search: id:A112622
|
|
|
| A112622 |
|
If p^b(p,n) is the highest power of the prime p dividing n, then a(n) = product_{p|n} b(p,n)^b(p,n). |
|
+0 1
|
|
| 1, 1, 1, 4, 1, 1, 1, 27, 4, 1, 1, 4, 1, 1, 1, 256, 1, 4, 1, 4, 1, 1, 1, 27, 4, 1, 27, 4, 1, 1, 1, 3125, 1, 1, 1, 16, 1, 1, 1, 27, 1, 1, 1, 4, 4, 1, 1, 256, 4, 4, 1, 4, 1, 27, 1, 27, 1, 1, 1, 4, 1, 1, 4, 46656, 1, 1, 1, 4, 1, 1, 1, 108, 1, 1, 4, 4, 1, 1, 1, 256, 256, 1, 1, 4, 1, 1, 1, 27, 1, 4, 1, 4
(list; graph; listen)
|
|
|
OFFSET
|
2,4
|
|
|
EXAMPLE
|
45 = 3^2 * 5^1. So a(45) = 2^2 * 1^1 = 4.
|
|
MATHEMATICA
|
f[n_] := Block[{fi = Last@Transpose@FactorInteger@n}, Times @@ (fi^fi)]; Rest@Array[f, 93] (* Robert G. Wilson v *)
|
|
PROGRAM
|
(PARI) a(n)=local(v, r, i); v=factorint(n); r=1; for(i=1, matsize(v)[1], r*=v[i, 2]^v[i, 2]); r (Herrgesell)
|
|
CROSSREFS
|
Sequence in context: A037291 A063851 A124777 this_sequence A119591 A010125 A010324
Adjacent sequences: A112619 A112620 A112621 this_sequence A112623 A112624 A112625
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Leroy Quet (qq-quet(AT)mindspring.com), Dec 25 2005
|
|
EXTENSIONS
|
More terms from Robert G. Wilson v (rgwv(AT)rgwv.com) and Lambert Herrgesell (zero815(AT)googlemail.com), Dec 27 2005
|
|
|
Search completed in 0.002 seconds
|