|
|
A280873
|
|
Numbers whose binary expansion does not begin 10 and does not contain 2 adjacent 0's; Ahnentafel numbers of X-chromosome inheritance of a male.
|
|
6
|
|
|
0, 1, 3, 6, 7, 13, 14, 15, 26, 27, 29, 30, 31, 53, 54, 55, 58, 59, 61, 62, 63, 106, 107, 109, 110, 111, 117, 118, 119, 122, 123, 125, 126, 127, 213, 214, 215, 218, 219, 221, 222, 223, 234, 235, 237, 238, 239, 245, 246, 247, 250, 251, 253, 254, 255
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
The number of ancestors at generation m from whom a living individual may have received an X chromosome allele is F_m, the m-th term of the Fibonacci Sequence.
The starting offset is zero (with a(0) = 0) for the same reason that we have A003714(0) = 0. Indeed, b(n) = A054429(A003714(n)) for n >= 0 yields the terms of this sequence, but in different order.
(End)
|
|
LINKS
|
|
|
FORMULA
|
|
|
MAPLE
|
gen[0]:= {0, 1, 3}:
gen[1]:= {6, 7}:
for n from 2 to 10 do
gen[n]:= map(t -> 2*t+1, gen[n-1]) union
map(t -> 2*t, select(type, gen[n-1], odd))
od:
sort(convert(`union`(seq(gen[i], i=0..10)), list)); # Robert Israel, Oct 11 2017
|
|
MATHEMATICA
|
male = {1, 3}; generations = 8;
Do[x = male[[i - 1]]; If[EvenQ[x],
male = Append[ male, 2*x + 1] ,
male = Flatten[Append[male, {2*x, 2*x + 1}]]]
, {i, 3, Fibonacci[generations + 1]}]; male
|
|
PROG
|
(PARI)
isA003754(n) = { n=bitor(n, n>>1)+1; n>>=valuation(n, 2); (n==1); }; \\ After Charles R Greathouse IV's Feb 06 2017 code.
isA280873(n) = (isA003754(n) && isA004760(n));
n=0; k=0; while(k <= 10946, if(isA280873(n), write("b280873.txt", k, " ", n); k=k+1); n=n+1; ); \\ Antti Karttunen, Oct 11 2017
(Python)
yield 1
if ((x & 1) and (x > 1)):
yield 2*x
yield 2*x+1
def take(n, g):
'''Returns a list composed of the next n elements returned by generator g.'''
z = []
if 0 == n: return(z)
for x in g:
z.append(x)
if n > 1: n = n-1
else: return(z)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
a(0) = 0 prepended and more descriptive alternative name added by Antti Karttunen, Oct 11 2017
|
|
STATUS
|
approved
|
|
|
|