[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A270269
Prime numbers with locations of right angle turns in the Ulam square spiral that are vertices of isosceles right triangles.
1
3, 5, 7, 31, 37, 43, 8011, 8101, 8191, 920641, 921601, 922561, 3894703, 3896677, 3898651, 5902471, 5904901, 5907331, 7450171, 7452901, 7455631, 7482961, 7485697, 7488433, 36066031, 36072037, 36078043, 37155121, 37161217, 37167313, 39759331, 39765637, 39771943
OFFSET
1,1
COMMENTS
See the illustration for more information.
Subsequence of A172979. This sequence is probably infinite.
An interesting property: the sequence of the differences between prime numbers that are vertices for each triangle is the sequence {2, 6, 90, 960, 1974, 2430, 2730, 2736, 6006, 6096, 6306, ...} = A087277: numbers n such that the three second-degree cyclotomic polynomials x^2 + 1, x^2 - x + 1 and x^2 + x + 1 are simultaneously prime.
For example:
2 = 5 - 3 = 7 - 5;
6 = 37 - 31 = 43 - 37;
90 = 8101 - 8011 = 8191 - 8101.
Consequence: a(3n) + A087277(n) is a square. The corresponding sequence of the squares is {3^2, 7^2, 91^2, 961^2, 1975^2, 2431^2, 2731^2, 2737^2, 6007^2, ...}.
Examples:
a(3) + A087277(1) = 7 + 2 = 3^2;
a(6) + A087277(2) = 43 + 6 = 7^2;
a(9) + A087277(3) = 8191 + 90 = 91^2.
MAPLE
nn:=20000:T:=array(1..nn):a0:=1:kk:=0:
for p from 1 to nn do :
a1:=a0+floor(p/2):a0:=a1:
if isprime(a1)
then
kk:=kk+1:T[kk]:=a1:
else
fi:
od:
for n from 1 to kk-2 do:
d1:=T[n+2]-T[n+1]:d2:=T[n+1]-T[n]:
if d1=d2
then
printf("%d %d %d \n", T[n], T[n+1], T[n+2]):
else
fi:
od:
CROSSREFS
Sequence in context: A288891 A358699 A322749 * A247104 A006378 A162714
KEYWORD
nonn
AUTHOR
Michel Lagneau, Mar 14 2016
STATUS
approved