|
|
A362434
|
|
Numbers that can be written as A000045(i) + j^2 for i,j>=0 in 4 ways.
|
|
3
|
|
|
17, 5185, 1669265, 537497857, 173072640401, 55728852710977, 17944517500293905, 5778078906241926145, 1860523463292399924497, 599082777101246533761601, 192902793703138091471310737, 62114100489633364207228295425, 20000547454868240136636039815825, 6440114166367083690632597592399937
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
A000045(1) and A000045(2) are counted separately, even though they both are 1.
|
|
LINKS
|
|
|
FORMULA
|
With a = A000045(6*k-1) and b = A000045(6*k) and x = 1 + b^2/4, we have
Conjecture: a(n) = 1 + A000045(6*n)^2/4.
|
|
EXAMPLE
|
|
|
MAPLE
|
N:= 10^6: # to get terms <= N
V:= Array(0..N, datatype=integer[1]):
for i from 0 do
f:= combinat:-fibonacci(i);
if f > N then break fi;
s:= floor(sqrt(N-f));
J:=[seq(f+i^2, i=0..s)];
V[J]:= V[J] +~ 1;
od:
select(i -> V[i] >= 4, [$1..N]);
|
|
PROG
|
(PARI) A362503(n) = my(f, s=0); for(i=0, oo, if(n<f=fibonacci(i), return(s), if(issquare(n-f), s++)));
lista(nn) = my(v=List([]), x, y, t); for(i=3, log(sqrt(5)*nn+1.5)\log((1+sqrt(5))/2), x=fibonacci(i); for(j=1, i-1, y=x-fibonacci(j); fordiv(y, d, if(d>sqrtint(y), break); t=y/d-d; if(t%2==0, for(k=0, j-1, if(issquare(t^2+4*(x-fibonacci(k))), listput(v, x+t^2/4))))))); v=Set(v); for(i=1, #v, if(v[i]>nn, break); if(A362503(v[i])==4, print1(v[i], ", "))); // Jinyuan Wang, Apr 24 2023
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|