login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A307640 Least number k such that n divides gcd(sigma(k), phi(k), tau(k)). 0
1, 3, 18, 15, 3344, 45, 24128, 30, 882, 3344, 1012736, 126, 1953792, 24128, 16200, 168, 452263936, 2016, 1852571648, 3344, 40768, 1012736, 27007123456, 420, 1490000, 1953792, 103968, 24128, 2739920699392, 30096, 8348342681600, 840, 9114624, 452263936, 6163776, 2016 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
For each n >= 1 there are infinitely many numbers s such that n divides sigma(s), phi(s) and tau(s).
From Dirichlet's theorem there are infinitely many numbers m for which the numbers p = n*m + 1 are prime. Then sigma(p^(n-1)), phi(p^(n-1)) and tau(p^(n-1)) numbers are divisible by n.
REFERENCES
Laurențiu Panaitopol, Alexandru Gica, Arithmetic problems and number theory. Ideas and methods of solving, Ed. Gil, Zalău, 2006, ch. 13, p. 79, pr. 18. (in Romanian).
LINKS
EXAMPLE
For n = 2, sigma(3) = 4, phi(3) = 2, tau(3) = 4 are divisible by 2.
For n = 5, sigma(3344) = 7440, phi (3344) = 1440, tau (3344) = 20 are divisible by 5 and by 10.
For n = 11, sigma(1012736) = 2161632 = 11 * 196512, phi(1012736) = 11 * 43008, tau(1012736) = 11 * 4 are divisible by 11.
MATHEMATICA
Array[Block[{i = 1}, While[Mod[GCD[DivisorSigma[1, i], EulerPhi@ i, DivisorSigma[0, i]], #] != 0, i++]; i] &, 16] (*Adaptation after A222713*)
PROG
(Magma) for m in [1..16] do
for n in [1..2000000] do
if IsIntegral(SumOfDivisors(n)/m) and IsIntegral(EulerPhi(n)/m) and IsIntegral(NumberOfDivisors(n)/m) then
m, n;
break;
end if;
end for;
end for;
(PARI) isok(n, k) = ! frac(gcd(sigma(k), gcd(eulerphi(k), numdiv(k)))/n);
a(n) = my(k=1); while(!isok(n, k), k++); k; \\ Michel Marcus, Apr 20 2019
(PARI) a(n) = {if(n==1, return(1)); my(res = oo, f = factor(n), hpf = f[#f~, 1]); forprime(p = 2, oo, if(p ^ (hpf - 1) > res, return(res)); forstep(i = p ^ (hpf - 1), res, p ^ (hpf - 1), if(isok(n, i), res = min(res, i); next(2) ) ) ) } \\ uses isok from above \\ David A. Corneth, Apr 22 2019
CROSSREFS
Sequence in context: A077104 A073815 A212994 * A339631 A195998 A291167
KEYWORD
nonn
AUTHOR
Marius A. Burtea, Apr 19 2019
EXTENSIONS
More terms from David A. Corneth, Apr 21 2019
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 23 13:51 EDT 2024. Contains 371914 sequences. (Running on oeis4.)