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.

A370682 Numbers k such that A370681(k) = A370681(k+1).

This page as a plain text file.
%I A370682 #7 Feb 27 2024 01:32:39
%S A370682 1,5,11,19,71,203,271,979,991,2755,5976,6407,7004,8560,9143,33552,
%T A370682 46583,75856,87955,90287,101955,243271,265587,269648,457135,510192,
%U A370682 512783,551172,559656,1804555,1917639,2902175,4339952,6783820,9842748,10463980,17428447,21677503
%N A370682 Numbers k such that A370681(k) = A370681(k+1).
%H A370682 Amiram Eldar, <a href="/A370682/b370682.txt">Table of n, a(n) for n = 1..120</a>
%e A370682 5 is a term since A370681(5) = A370681(6) = 4.
%t A370682 s[n_] := Module[{d = Reverse[Select[Divisors[n], CoprimeQ[#, n/#] &]]}, Total[(-1)^(Range[Length[d]] + 1)*d]]; Select[Range[10^5], s[#] == s[# + 1] &]
%o A370682 (PARI)  s(n) = {my(d = Vecrev(select(x->(gcd(x, n/x) == 1), divisors(n)))); sum(i=1, #d, (-1)^(i+1)*d[i]);}
%o A370682 lista(nmax) = {my(s1 = s(1), s2); for(n = 2, nmax, s2 = s(n); if(s1 == s2, print1(n-1,", ")); s1=s2);}
%Y A370682 Cf. A064125, A333261, A370681.
%K A370682 nonn
%O A370682 1,2
%A A370682 _Amiram Eldar_, Feb 26 2024