login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A182580 Position of first occurrence of n in A182576. 1
0, 1, 2, 5, 11, 27, 29, 46, 76, 167, 133, 348, 521, 1160, 1364, 3005, 3571, 9348, 9349, 20820, 24476, 64062, 64079, 123651, 167761, 439203, 439204, 1149850, 1149851, 3010348, 3010349, 7881195, 7881196, 20633238, 20633239, 54018520, 54018521, 141422323, 141422324 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
Conjecture: for k > 12, a(2k) = a(2k-1) + 1.
LINKS
EXAMPLE
First occurrence of 4 in A182576(n) has index 11, so a(4)=11.
PROG
(C)
#include <stdio.h> // GCC // ~880 secs
typedef unsigned long long U64;
U64 fibs[992], first[992], prpr=0, prev=1, current, n, sq;
long long ftop, fpos, terms;
int main(int argc, char **argv)
{
for (ftop=0; ftop<900; ++ftop) {
fibs[ftop] = current = prpr+prev;
if (current<prev) break;
prpr = prev, prev = current;
}
printf("%llu\n%llu\n%llu\n%llu\n", ftop, prev, prpr, (U64)-1-prev);
for (n=0; n < ((U64)1<<32); ++n) {
fpos=ftop-1, terms=0;
for (sq=n*n; sq>0 && fpos>=0; ++terms, sq-=fibs[fpos], --fpos) {
while (fibs[fpos]>sq && fpos>0) --fpos;
}
if (sq>0) { printf("Error!"); exit(1); }
if (first[terms]==0)
first[terms]=n, printf("%2llu %llu\n", terms, n);
}
for (n=0; n<ftop; ++n) printf("%llu, ", first[n]);
return 0;
}
CROSSREFS
Cf. A182576.
Sequence in context: A000664 A364597 A242766 * A067922 A212567 A339755
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, May 05 2012
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 08:27 EDT 2024. Contains 371964 sequences. (Running on oeis4.)