|
|
A057142
|
|
Occurrences of most frequently occurring number in 1-to-n multiplication table.
|
|
8
|
|
|
1, 2, 2, 3, 3, 4, 4, 4, 4, 4, 4, 6, 6, 6, 6, 6, 6, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 16, 16, 16, 16, 16, 16, 16, 16
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
EXAMPLE
|
M(n) is the array in which m(x,y)= x*y for x = 1 to n and y = 1 to n. In m(5), the most frequently occurring number is 4. It occurs 3 times, so a(5) = 3.
|
|
PROG
|
(Haskell)
import Data.List (group, sort)
a057142 n = head $ reverse $ sort $ map length $ group $
sort [u * v | u <- [1..n], v <- [1..n]]
(PARI) T(n, f=factor(n))=my(k=#f~); f[, 1]=primes(k+1)[2..k+1]~; f[1, 1]=6; factorback(f)
listA025487(Nmax)=vecsort(concat(vector(logint(Nmax, 2), n, select(t->t<=Nmax, if(n>1, [factorback(primes(#p), Vecrev(p))|p<-partitions(n)], [1, 2])))))
ct(n, k)=sumdiv(n, d, max(d, n/d)<=k)
a(n)=if(n==1, return(1)); my(v=listA025487(n^2), r, t); for(i=1, #v, t=ct(v[i], n); if(t>r, r=t)); r \\ Charles R Greathouse IV, Feb 05 2022
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
More terms from Larry Reeves (larryr(AT)acm.org), Apr 18 2001
|
|
STATUS
|
approved
|
|
|
|