[go: up one dir, main page]

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”).

A093506
Variation on Golomb's sequence starting with (1,2): a(n)=length of n-th run of consecutive integers with same parity.
5
1, 2, 4, 3, 5, 7, 9, 6, 8, 10, 11, 13, 15, 17, 19, 12, 14, 16, 18, 20, 22, 24, 21, 23, 25, 27, 29, 31, 33, 35, 37, 26, 28, 30, 32, 34, 36, 39, 41, 43, 45, 47, 49, 51, 53, 38, 40, 42, 44, 46, 48, 50, 52, 54, 56, 55, 57, 59, 61, 63, 65, 67, 69, 71, 73, 75, 58, 60, 62, 64, 66, 68, 70
OFFSET
1,2
COMMENTS
A permutation of positive integers.
FORMULA
Does a(n)=n+o(n)? Does a(n)=n for infinitely many values of n?
EXAMPLE
Runs of odd or even terms are : (1),(2,4),(3,5,7,9),(6,8,10),(11,13,15,17,19),...and computing the number of integers in each run we get : 1,2,4,3,5,...the sequence itself.
PROG
(GNU bc)
maxarray=2^16;
print oe[1]=a[1]=1, ", ", a[2]=2, ", ", oe[k=0]=a[p=n=3]=4, ", ";
for(max=100; p<maxarray&&n<max; p++) {
for(i=0*k=!k; i<a[p]&&n<max; i++) {
an=oe[k]+=2;
if (++n<maxarray) a[n]=an;
print an, ", "
};
}; /* Carl R. White, Jan 05 2013 */
(Haskell)
a093506 n = a093506_list !! (n-1)
a093506_list = 1 : 2 : f 1 [1] [3, 5..] [4, 6..]
where f 0 (z:zs) odds evens = orun ++ f 1 (zs ++ orun) odds' evens
where (orun, odds') = splitAt z odds
f 1 (z:zs) odds evens = erun ++ f 0 (zs ++ erun) odds evens'
where (erun, evens') = splitAt z evens
-- Reinhard Zumkeller, Jan 06 2013
CROSSREFS
Cf. A001462.
Cf. A187790 (inverse), A187792 (fixed points).
Sequence in context: A348357 A199779 A374800 * A307485 A238980 A327143
KEYWORD
nonn
AUTHOR
Benoit Cloitre, May 14 2004
EXTENSIONS
Sequence corrected by Carl R. White, Jan 06 2013
STATUS
approved