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.

A131559 Indices of records in A078571.

This page as a plain text file.
%I A131559 #21 Sep 22 2024 06:15:25
%S A131559 1,3,7,8,14,41,270,277,1595,5899,6320,6668,117221,180037,295677,
%T A131559 587152,703625,4397877,12151726,29573235
%N A131559 Indices of records in A078571.
%C A131559 Numbers k such that total number of prime factors (with multiplicity) of the average of k-th twin prime pair is a record.
%C A131559 The records are 1, 2, 3, 4, 5, 7, 9, 10, 12, 14, 15, 17, 19, ... and the associated twin primes recorded in A086827. - _R. J. Mathar_, Oct 24 2007
%F A131559 a(n) = {k such that A001222(A014574(k)) is a record}.
%e A131559 a(1) = 1 because A078571(1) = 2, since the average of twin prime pair (3,5) is 4 = 2^2 is a semiprime with two prime factors (with multiplicity).
%e A131559 a(2) = 3 because A078571(3) = smallomega(12) = 3.
%e A131559 a(3) = 7 because A078571(7) = smallomega(60) = 4.
%e A131559 a(4) = 8 because A078571(8) = smallomega(72) = 5.
%e A131559 a(5) = 14 because A078571(14) = smallomega(192) = 7.
%e A131559 a(6) = 41 because A078571(41) = smallomega(1152) = smallomega(2^7 * 3^2) = 9.
%o A131559 (PARI) lista(nn) = {my(r=0, c=0, p=2); forprime(q=3, nn, if(q-p==2, c++; if(bigomega(p+1)>r, r=bigomega(p+1); print1(c, ", "))); p=q); } \\ _Jinyuan Wang_, Apr 01 2020
%Y A131559 Cf. A001221, A001359, A006512, A014574, A078571.
%K A131559 nonn,more
%O A131559 1,2
%A A131559 _Jonathan Vos Post_, Oct 02 2007
%E A131559 More terms from _R. J. Mathar_, Oct 24 2007
%E A131559 a(13)-a(20) from _Jinyuan Wang_, Apr 01 2020