|
|
A067631
|
|
If n is composite then a(n) is the standard deviation of the prime factors of n, rounded off to the nearest integer (rounding up if there's a choice), with each factor counted according to its frequency of occurrence in the prime factorization. If n is 1 or prime then a(n)=0.
|
|
1
|
|
|
0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 4, 1, 0, 0, 1, 0, 2, 3, 6, 0, 1, 0, 8, 0, 3, 0, 2, 0, 0, 6, 11, 1, 1, 0, 12, 7, 2, 0, 3, 0, 5, 1, 15, 0, 0, 0, 2, 10, 6, 0, 1, 4, 3, 11, 19, 0, 1, 0, 21, 2, 0, 6, 5, 0, 9, 14, 3, 0, 1, 0, 25, 1, 10, 3, 6, 0, 1, 0, 28, 0, 2, 8, 29, 18, 5, 0, 1, 4, 12, 20, 32, 10, 0, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
2,9
|
|
COMMENTS
|
The (sample) standard deviation sigma of {x_1,...,x_n} is calculated from sigma^2 = 1/(n-1) * sum_{1,...,n}(x_i - mu)^2, where mu denotes the average of {x_1,...,x_n}.
|
|
LINKS
|
|
|
EXAMPLE
|
24 = 2^3 * 3^1, so the corresponding average = (2 + 2 + 2 + 3)/ 4 = 2.25 and the standard deviation is [(1/3){3 * (2-2.25)^2 + (3-2.25)^2}]^0.5 = 0.5, which rounds to 1. So a(24) = 1.
|
|
MATHEMATICA
|
<<Statistics`NormalDistribution` f[n_] := Flatten[Table[ #[[1]], {#[[2]]}]&/@FactorInteger[n]]; a[n_] := If[PrimeQ[n]||n==1, 0, Floor[StandardDeviation[f[n]]+1/2]]
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|