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.

Showing 1-6 of 6 results.

A127490 Primes in A127489.

Original entry on oeis.org

2927, 342889, 2760049, 5070049, 321039529, 1146272329, 9936508189, 14927726089, 37769930377, 79361849257, 222647590729, 246738487669, 259512888949, 636714780049, 818804190829, 1204024293601, 1265682003121
Offset: 1

Views

Author

Artur Jasinski, Jan 16 2007

Keywords

Crossrefs

Cf. A127489.

Programs

  • Mathematica
    Select[Table[CoefficientList[Expand[(x-Prime[n])*(x-Prime[n+1])*(x-Prime[n+2])* (x-Prime[n+3])*(x-Prime[n+4])],x][[2]],{n,1,130}], PrimeQ[ # ]&]

Extensions

Edited by Stefan Steinerberger, Jul 18 2007

A127491 Primes which are half of the absolute coefficients [x^2] of the 5th-order polynomials with prime roots as defined in A127489.

Original entry on oeis.org

310733, 426871, 15722159, 166492163, 177861107, 270396557, 342955763, 406947461, 1606837039, 1908243773, 2902193117, 3386269021, 5441167877, 6953015807, 7671152921, 10005413687, 10979785673, 14774655421, 16546239937
Offset: 1

Views

Author

Artur Jasinski, Jan 16 2007

Keywords

Comments

The polynomials are of the form (x-prime(i))*(x-prime(i+1))*..*(x-prime(i+4)). The quadratic terms have coefficients which are of the form -sum_{j

Examples

			The first contribution is from the 11th polynomial, (x-prime(11)) *(x-prime(12)) *(x-prime(13)) *(x-prime(14)) *(x-prime(15)) = x^5 -199x^4 +15766x^3 -621466x^2 +12185065x -95041567,
where the coefficient of [x^2] is -621466. Its sign-reversed half is 310733, a prime.
		

Programs

  • Maple
    isA127491 := proc(k)
        local x,j,p ;
        mul( x-ithprime(k+j),j=0..4) ;
        expand(%) ;
        abs(coeff(%,x,2)/2) ;
        isprime(%)
    end proc:
    A127491k := proc(n)
        option remember ;
        if n = 0 then
            0;
        else
            for k from procname(n-1)+1 do
                if isA127491(k) then
                    return k ;
                end if;
            end do:
        end if;
    end proc:
    A127491 := proc(n)
        option remember ;
        local k ;
        k := A127491k(n) ;
        mul( x-ithprime(k+j),j=0..4) ;
        expand(%) ;
        abs(coeff(%,x,2)/2) ;
    end proc:
    seq(A127491(n),n=1..60) ; # R. J. Mathar, Apr 23 2023

Extensions

Entries replaced to comply with the definition. - R. J. Mathar, Sep 26 2011

A127492 Indices m of primes such that Sum_{k=0..2, k

Original entry on oeis.org

2, 10, 17, 49, 71, 72, 75, 145, 161, 167, 170, 184, 244, 250, 257, 266, 267, 282, 286, 301, 307, 325, 343, 391, 405, 429, 450, 537, 556, 561, 584, 685, 710, 743, 790, 835, 861, 904, 928, 953
Offset: 1

Author

Artur Jasinski, Jan 16 2007

Keywords

Comments

Let p_0 .. p_4 be five consecutive primes, starting with the m-th prime. The index m is in the sequence if the absolute value [x^0] of the polynomial (x-p_0)*[(x-p_1)*(x-p_2) + (x-p_2)*(x-p_3) + (x-p_3)*(x-p_4)] + (x-p_1)*[(x-p_2)*(x-p_3) + (x-p_3)*(x-p_4)] + (x-p_2)*(x-p_3)*(x-p_4) is two times a prime. The correspondence with A127491: the coefficient [x^2] of the polynomial (x-p_0)*(x-p_1)*..*(x-p_4) is the sum of 10 products of a set of 3 out of the 5 primes. Here the sum is restricted to the 6 products where the two largest of the 3 primes are consecutive. - R. J. Mathar, Apr 23 2023

Programs

  • Maple
    isA127492 := proc(k)
        local x,j ;
        (x-ithprime(k))* mul( x-ithprime(k+j),j=1..2)
        +(x-ithprime(k))* mul( x-ithprime(k+j),j=2..3)
        +(x-ithprime(k))* mul( x-ithprime(k+j),j=3..4)
        +(x-ithprime(k+1))* mul( x-ithprime(k+j),j=2..3)
        +(x-ithprime(k+1))* mul( x-ithprime(k+j),j=3..4)
        +(x-ithprime(k+2))* mul( x-ithprime(k+j),j=3..4) ;
        p := abs(coeff(expand(%/2),x,0)) ;
        if type(p,'integer') then
            isprime(p) ;
        else
            false ;
        end if ;
    end proc:
    for k from 1 to 900 do
        if isA127492(k) then
            printf("%a,",k) ;
        end if ;
    end do: # R. J. Mathar, Apr 23 2023
  • Mathematica
    a = {}; Do[If[PrimeQ[(Prime[x] Prime[x + 1]Prime[x + 2] + Prime[x] Prime[x + 2]Prime[x + 3] + Prime[x] Prime[x + 3] Prime[x + 4] + Prime[x + 1] Prime[x + 2]Prime[x + 3] + Prime[x + 1] Prime[x + 3]Prime[x + 4] + Prime[x + 2] Prime[x + 3] Prime[x + 4])/2], AppendTo[a, x]], {x, 1, 1000}]; a
    prQ[{a_,b_,c_,d_,e_}]:=PrimeQ[(a b c+a c d+a d e+b c d+b d e+c d e)/2]; PrimePi/@Select[ Partition[ Prime[Range[1000]],5,1],prQ][[;;,1]] (* Harvey P. Dale, Apr 21 2023 *)

Extensions

Definition simplified by R. J. Mathar, Apr 23 2023
Edited by Jon E. Schoenfield, Jul 23 2023

A127493 Indices k such that the coefficient [x^1] of the polynomial Product_{j=0..4} (x-prime(k+j)) is prime.

Original entry on oeis.org

1, 5, 8, 9, 22, 29, 45, 49, 60, 69, 87, 89, 90, 107, 114, 124, 125, 131, 134, 138, 145, 156, 171, 183, 188, 191, 203, 204, 207, 212, 219, 255, 261, 290, 298, 303, 329, 330, 343, 344, 349, 354, 378, 397, 398, 400, 403, 454, 456, 466, 474, 515, 549, 560, 570, 578
Offset: 1

Author

Artur Jasinski, Jan 16 2007

Keywords

Comments

A fifth-order polynomial with 5 roots which are the five consecutive primes from prime(k) onward is defined by Product_{j=0..4} (x-prime(k+j)). The sequence is a catalog of the cases where the coefficient of its linear term is prime.
Indices k such that e4(prime(k), prime(k+1), ..., prime(k+4)) is prime, where e4 is the elementary symmetric polynomial summing all products of four variables. - Charles R Greathouse IV, Jun 15 2015

Examples

			For k=2, the polynomial is (x-3)*(x-5)*(x-7)*(x-11)*(x-13) = x^5-39*x^4+574*x^3-3954*x^2+12673*x-15015, where 12673 is not prime, so k=2 is not in the sequence.
For k=5, the polynomial is x^5-83*x^4+2710*x^3-43490*x^2+342889*x-1062347, where 342889 is prime, so k=5 is in the sequence.
		

Programs

  • Maple
    isA127493 := proc(k)
        local x,j ;
        mul( x-ithprime(k+j),j=0..4) ;
        expand(%) ;
        isprime(coeff(%,x,1)) ;
    end proc:
    A127493 := proc(n)
        option remember ;
        if n = 1 then
            1;
        else
            for a from procname(n-1)+1 do
                if isA127493(a) then
                    return a;
                end if;
            end do:
        end if;
    end proc:
    seq(A127493(n),n=1..60) ; # R. J. Mathar, Apr 23 2023
  • Mathematica
    a = {}; Do[If[PrimeQ[(Prime[x] Prime[x + 1]Prime[x + 2]Prime[x + 3] + Prime[x] Prime[x + 2]Prime[x + 3]Prime[x + 4] + Prime[x] Prime[x + 1]Prime[x + 3]Prime[x + 4] + Prime[x] Prime[x + 1]Prime[x + 2]Prime[x + 4] + Prime[x + 1] Prime[x + 2]Prime[x + 3]Prime[x + 4])], AppendTo[a, x]], {x, 1, 1000}]; a
  • PARI
    e4(v)=sum(i=1,#v-3, v[i]*sum(j=i+1,#v-2, v[j]*sum(k=j+1,#v-1, v[k]*vecsum(v[k+1..#v]))))
    pr(p, n)=my(v=vector(n)); v[1]=p; for(i=2,#v, v[i]=nextprime(v[i-1]+1)); v
    is(n,p=prime(n))=isprime(e4(pr(p,5)))
    v=List(); n=0; forprime(p=2,1e4, if(is(n++,p), listput(v,n))); Vec(v) \\ Charles R Greathouse IV, Jun 15 2015

Extensions

Definition and comment rephrased and examples added by R. J. Mathar, Oct 01 2009

A125270 Coefficient of x^2 in polynomial whose zeros are 5 consecutive primes starting with the n-th prime.

Original entry on oeis.org

1358, 3954, 10478, 22210, 43490, 78014, 129530, 206650, 324350, 466270, 621466, 853742, 1132130, 1436690, 1870850, 2388050, 2886370, 3440410, 4133410, 4904906, 5926654, 7195670, 8425430, 9792950, 11040910, 12098990, 13917898, 16097810
Offset: 1

Author

Artur Jasinski, Jan 16 2007

Keywords

Comments

Sums of all distinct products of 3 out of 5 consecutive primes, starting with the n-th prime; value of 3rd elementary symmetric function on the 5 consecutive primes.

Programs

  • Mathematica
    a = {}; Do[AppendTo[a, (Prime[x] Prime[x + 1] Prime[x + 2] + Prime[x] Prime[x + 1] Prime[x + 3] + Prime[x] Prime[x + 1] Prime[x + 4] + Prime[x] Prime[x + 2] Prime[x + 3] + Prime[x] Prime[x + 2] Prime[x + 4] + Prime[x] Prime[x + 3] Prime[x + 4] + Prime[x + 1] Prime[x + 2] Prime[x + 3] + Prime[x + 1] Prime[x + 2] Prime[x + 4] + Prime[x + 1] Prime[x + 3] Prime[x + 4] + Prime[x + 2] Prime[x + 3] Prime[x + 4])], {x, 1, 100}]; a
    fcp[{p_,q_,r_,s_,t_}]:=p*q(r+s+t)+(p+q)r(s+t)+(p+q+r)s*t; fcp/@Partition[ Prime[ Range[40]],5,1] (* Harvey P. Dale, Sep 05 2014 *)

Formula

Let p = Prime(n), q = Prime(n+1), r = Prime(n+2), s = Prime(n+3) and t = Prime(n+4). Then a(n) = p q (r+s+t) + (p + q) r (s + t) + (p + q + r) s t.

Extensions

Edited and corrected by Franklin T. Adams-Watters, Jan 23 2007

A191216 Arithmetic derivative of prime(n) * prime(n+1) * prime(n+2) * prime(n+3) * prime(n+4) * prime(n+5).

Original entry on oeis.org

361, 230456, 1005768, 3462570, 11006128, 25925028, 61456764, 127697940, 249379116, 448408452, 740850012, 1263239320, 1914568816, 2884222410, 4371191782, 6287341056, 8758591370, 11640682466, 15938770638, 21721208748, 29153150298, 38784336168, 49888704100, 62506263054, 76188213990, 95511276660, 118760260290, 150724895476, 187405610004, 243040520764
Offset: 1

Author

Giorgio Balzarotti, May 26 2011

Keywords

Crossrefs

Programs

  • Maple
    der:=n->n*add(op(2,p)/op(1,p),p=ifactors(n)[2]);
    seq(dif(ithprime(i)*ithprime(i+1)*ithprime(i+2)*ithprime(i+3)*ithprime(i+4)*ithprime(i+5)),i=1..30);

Formula

a(n) = (prime(n) * prime(n+1) * prime(n+2) * prime(n+3) * prime(n+4) * prime(n+5))' where f' is the arithmetic derivative (see A003415) of f.
Showing 1-6 of 6 results.