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.

A276533 Least prime p with A271518(p) = n.

Original entry on oeis.org

5, 2, 19, 127, 17, 67, 163, 41, 89, 101, 131, 313, 257, 211, 227, 461, 241, 401, 613, 337, 433, 353, 577, 467, 863, 887, 617, 787, 601, 569, 761, 641, 823, 673, 857, 1217, 881, 1091, 1289, 977, 1427, 1097, 1801, 929, 1153, 953, 1321, 1049, 1747, 1409
Offset: 1

Views

Author

Zhi-Wei Sun, Dec 12 2016

Keywords

Comments

Conjecture: a(n) exists for any positive integer n.
In contrast, it is known that for each prime p the number of ordered integral solutions to the equation x^2 + y^2 + z^2 + w^2 = p is 8*(p+1).
In 1998 J. Friedlander and H. Iwaniec proved that there are infinitely many primes p of the form w^2 + x^4 = w^2 + (x^2)^2 + 0^2 + 0^2 with w and x nonnegative integers. Since x^2 + 3*0 + 5*0 is a square, we see that A271518(p) > 0 for infinitely many primes p.

Examples

			a(1) = 5 since 5 is the first prime which can be written in a unique way as x^2 + y^2 + z^2 + w^2 with x,y,z,w nonnegative integer and x + 3*y + 5*z a square; in fact, 5 = 1^2 + 0^2 + 0^2 + 2^2 with 1 + 3*0 + 5*0 = 1^2.
a(2) = 2 since 2 = 1^2 + 0^2 + 0^2 + 1^2 with 1 + 3*0 + 5*0 = 1^2, and 2 = 1^2 + 1^2 + 0^2 + 0^2 with 1 + 3*1 + 5*0 = 2^2.
		

Crossrefs

Programs

  • Mathematica
    SQ[n_]:=SQ[n]=IntegerQ[Sqrt[n]];
    Do[m=0;Label[aa];m=m+1;r=0;Do[If[SQ[Prime[m]-x^2-y^2-z^2]&&SQ[x+3y+5z],r=r+1;If[r>n,Goto[aa]]],{x,0,Sqrt[Prime[m]]},{y,0,Sqrt[Prime[m]-x^2]},{z,0,Sqrt[Prime[m]-x^2-y^2]}];If[r