[seqfan] Re: A120511 a(n) = min{j : A006949(j) = n}
Frank Ruskey
ruskey at cs.uvic.ca
Sat Sep 12 01:42:07 CEST 2009
Please note the following equation from an updated version of
the Deugau-Ruskey paper listed in A120511 that appeared in
JIS: http://www.cs.uwaterloo.ca/journals/JIS/VOL12/Ruskey/ruskey6.pdf
See page 25.
p(n) = s*ceil(log_k n) + (kn-d-1)/(k-1)
Here, d is the sum of the digits of the k-ary expression of n-1.
In your case s = 1 and k = 2. Nice to see that somebody is
actually looking at these sequences!
Best Regards,
Frank
Max Alekseyev wrote:
> On Fri, Sep 11, 2009 at 8:11 AM, Georgi Guninski <guninski at guninski.com> wrote:
>
>>> I guess you meant a(n) = a(ceil(n/2)) + n.
>>>
>> thank you.
>>
>> assuming A120511 is efficiently computable (your proof or the potential
>> doubling formulas i found) is A006949 efficiently computable by
>> something like a binary search in A120511?
>
> It is possible.
>
> Actually, from the recurrence for A120511(n) we can derive even an
> explicit formula:
>
> Let z(n) be the number of zero bits in the binary representation of n. Then
>
> A120511(n) = 2n + z(n) - k - [n==2^k],
>
> where k = valuation(n,2), i.e., the maximum power of 2 dividing n.
>
> Here is a corresponding PARI/GP code:
>
> { A120511(n) = local(t,k); t=binary(n); k=valuation(n,2); 2*n + #t -
> sum(i=1,#t,t[i]) - k - (n==2^k) }
>
> Note that k <= z(n) <= log2(n)-1, implying that
> 2n-1 <= A120511(n) <= 2n + log2(n) - 1
>
> Since A006949(m) equals the largest n such that A120511(n) <= m (and
> thus A120511(n+1) > m), from 2n-1 <= A120511(n) <= m it follows that
> A006949(m) <= (m+1)/2.
> Similarly, from
> m < A120511(n+1) < 2(n+1) + log2(n+1) - 1 <= 2(n+1) + log2((m+1)/2+1) - 1,
> it follows that
> A006949(m) >= (m - log2(m+3)) / 2.
>
> Therefore,
> | A006949(m) - m/2 | <= log2(m+3)/2
> that gives an interval of just logarithmic length to search the value
> of A006949(m).
>
> Regards,
> Max
>
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/
>
--
----------------------
Frank Ruskey e-mail: (last_name)(AT)cs(DOT)uvic(DOT)ca
Dept. of Computer Science office: 250-472-5794
University of Victoria fax: 250-472-5708
Victoria, B.C. V8W 3P6 CANADA WWW: http://www.cs.uvic.ca/~(last_name)
--
This message has been scanned for viruses and
dangerous content by MailScanner, and is
believed to be clean.
More information about the SeqFan
mailing list