OFFSET
1,1
COMMENTS
Conjecture 1: If A062402(n) = A000203(A000010(n)) = sigma(phi(n)) is a prime p for some n, then p is Mersenne prime (A000668); a(n) > 0 for all n.
Conjecture 2: a(n) = the smallest number k such that phi(k) has exactly A000043(n)-1 divisors; see A276044.
Conjecture 3: a(n) = the smallest number k such that phi(k) has exactly A000043(n)-1 prime factors (counted with multiplicity); see A275969.
Equals A002181 (index in A002202 of (intersection of A023194 and A002202)). - Michel Marcus, Feb 12 2017
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..13
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
PROG
(Magma) A281627:=func<n|exists(r){k:k in[1..1000000] | SumOfDivisors(EulerPhi(k)) eq n}select r else 0>; Set(Sort([A281627(n): n in [SumOfDivisors(EulerPhi(n)): n in[1..1000000] | IsPrime(SumOfDivisors(EulerPhi(n)))]]))
(PARI) terms() = {v = readvec("b023194.txt"); for(i=1, #v, if (istotient(v[i], &n), print1(n/2, ", ")); ); } \\ Michel Marcus, Feb 12 2017
(PARI) f(p) = {my(s = invsigma(p), t, minv = 0); for(i = 1 , #s, t = invphi(s[i]); for(j = 1, #t, if(minv == 0, minv = t[j]); if(t[j] < minv, minv = t[j]))); minv; } \\ using Max Alekseyev's invphi.gp
list(pmax) = forprime(p = 1, pmax, if(isprime(2^p-1), print1(f(2^p-1), ", "))); \\ Amiram Eldar, Dec 23 2024
CROSSREFS
Cf. A053576 (includes the first 13 known terms of this sequence).
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Feb 11 2017
EXTENSIONS
a(8) from Michel Marcus, Feb 12 2017
a(9)-a(12) from Amiram Eldar, Dec 23 2024
STATUS
approved