cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A079428 Least coprime number m > n with same number of divisors as n, a(1) = 1.

This page as a plain text file.
%I A079428 #9 Mar 26 2025 08:30:44
%S A079428 1,3,5,9,7,35,11,15,25,21,13,175,17,15,22,81,19,175,23,63,22,27,29,
%T A079428 385,49,27,34,45,31,1001,37,45,34,35,38,1225,41,39,46,189,43,715,47,
%U A079428 45,52,51,53,4375,121,63,55,63,59,385,57,135,58,65,61,7007,67,65,68,729,69
%N A079428 Least coprime number m > n with same number of divisors as n, a(1) = 1.
%H A079428 Amiram Eldar, <a href="/A079428/b079428.txt">Table of n, a(n) for n = 1..10000</a>
%F A079428 tau(a(n)) = tau(n), where tau = A000005.
%F A079428 a(A000040(k)) = A079427(A000040(k)) = A000040(k+1).
%t A079428 a[n_] := Module[{d = DivisorSigma[0, n], m = n + 1}, While[! CoprimeQ[m, n] || DivisorSigma[0, m] != d, m++]; m]; a[1] = 1; Array[a, 100] (* _Amiram Eldar_, Mar 26 2025 *)
%o A079428 (PARI) a(n) = if(n == 1, 1, my(d = numdiv(n), m = n + 1); while(gcd(m, n) > 1 || numdiv(m) != d, m++); m); \\ _Amiram Eldar_, Mar 26 2025
%Y A079428 Cf. A000005, A000040, A079427.
%K A079428 nonn
%O A079428 1,2
%A A079428 _Reinhard Zumkeller_, Jan 08 2003