A074964 Numbers k such that Max ( sigma(x*y) : 1 <= x <= k, 1 <= y <= k ) = sigma(k^2).
1, 2, 3, 4, 6, 8, 12, 18, 24, 60
Offset: 1
Programs
-
Haskell
a074964 n = a074964_list !! (n-1) a074964_list = filter (\x -> a074963 x == a065764 x) [1..] -- Reinhard Zumkeller, Nov 14 2011
-
Maple
with(numtheory): s := proc(n) option remember: return sigma(n): end: a:= proc(n) option remember: if(n=0)then return 0: fi: return max(a(n-1),seq(s(x*n),x=1..n)): end: for n from 1 to 100 do if(a(n)=s(n^2))then printf("%d, ", n): end: od: # Nathaniel Johnston, Sep 26 2011
-
PARI
isok(k) = vecmax(setbinop((x,y)->sigma(x*y), [1..k])) == sigma(k^2); \\ Michel Marcus, Feb 03 2022
Comments