OFFSET
1,2
COMMENTS
Equivalent descriptions of the centrality of n: 1) Probability that a randomly chosen product in the multiplication table for positive integers (A003991; see also A061017) is a multiple of n.
2) Probability taken over all exponential numerical bases that if the last digit of a number represents n, the number is a multiple of n. (For example, in base 10, the probability of a number that ends in 5 being a multiple of 5 is 1. Over all possible bases, the fraction of numbers ending in 5 that are multiples of 5 is the centrality of 5, 9/25 or .36.)
An infinite number of integers have the same centrality as at least one other integer. The only such examples in the first 114 terms of the sequence are 64 and 120, which share a centrality of .0625; they are listed in numerical order.
LINKS
T. D. Noe, Table of n, a(n) for n=1..1000
FORMULA
EXAMPLE
The number 6 has a gcd of 1 with all numbers congruent to 1 or 5 modulo 6, 2 with all numbers congruent to 2 or 4 mod 6, 3 with all 3 mod 6 numbers and 6 with all numbers congruent to 0 mod 6. Its average gcd with other integers is 2.5 (A018804(6)/6), which represents 5/12 or .41666... of 6. This places 6 fifth in centrality among the integers, behind 1 (whose centrality is 1), 2 (.75), 3 (5/9 or .555...) and 4 (.5); it is therefore listed fifth in the sequence.
MATHEMATICA
maxTerms = 100; Clear[c, s]; c[n_] := c[n] = Sum[d*EulerPhi[n/d], {d, Divisors[n] }]/n^2; s[terms_] := s[terms] = Sort[Range[terms], c[#1] >= c[#2] & ][[1 ;; maxTerms]]; s[terms = maxTerms]; s[terms += maxTerms]; While[s[terms] != s[terms - maxTerms], terms += maxTerms]; A080997 = s[terms] (* Jean-François Alcover, Feb 19 2015 *)
CROSSREFS
KEYWORD
nice,nonn
AUTHOR
Matthew Vandermast, Feb 28 2003
STATUS
approved