A192577 Numbers n such that the arithmetic mean of the unitary divisors of n is a prime number.
3, 5, 6, 9, 12, 13, 25, 37, 48, 61, 73, 81, 121, 157, 193, 277, 313, 361, 397, 421, 457, 541, 613, 625, 661, 673, 733, 757, 768, 841, 877, 997, 1093, 1153, 1201, 1213, 1237, 1321, 1381, 1453, 1621, 1657, 1753, 1873, 1933, 1993, 2017, 2137, 2341, 2401, 2473
Offset: 1
Keywords
Examples
48 has unitary divisors 1, 3, 16, 48 and (1+3+16+48)/4 = 17 is prime, therefore 48 is in the sequence.
Links
- Klaus Brockhaus, Table of n, a(n) for n = 1..10000
- A. Roldan Martinez, Numeros y hoja de calculo
Programs
-
Magma
UnitaryDivisors:=func< n | [ d: d in Divisors(n) | Gcd(d, n div d) eq 1 ] >; [ n: n in [1..2500] | IsPrime(k) and s mod #U eq 0 where k is s div #U where s is &+U where U is UnitaryDivisors(n) ]; // Klaus Brockhaus, Jul 09 2011
-
PARI
usigma(n)= {local(f, u=1); f=factor(n); for(i=1, matsize(f)[1], u*=(1+ f[i, 1]^f[i, 2])); return(u)} ud(n)= {local (f, u); f=factor(n); u=2^(matsize(f)[1]); return(u) } { for (n=2, 10^4, c=usigma(n)/ud(n); if (c==truncate(c),if(isprime(c), print1(n, ", ")))) } \\ Antonio Roldán, Oct 08 2012
Comments