Пусть
![$$\ell(n)=\left\lfloor\log_2 n\right\rfloor$$ $$\ell(n)=\left\lfloor\log_2 n\right\rfloor$$](https://dxdy-02.korotkov.co.uk/f/5/3/7/537bd544bf52dcef583f0e05e91668c582.png)
Пусть
![$$T(n,k)=\left\lfloor\frac{n}{2^k}\right\rfloor\operatorname{mod}2$$ $$T(n,k)=\left\lfloor\frac{n}{2^k}\right\rfloor\operatorname{mod}2$$](https://dxdy-03.korotkov.co.uk/f/a/5/0/a5016bd8aaa841d829469eb247d4abd182.png)
Здесь
![$T(n,k)$ $T(n,k)$](https://dxdy-02.korotkov.co.uk/f/d/7/7/d77ae0e85c7707d2449584121bb2c55d82.png)
это
![$(k+1)$ $(k+1)$](https://dxdy-01.korotkov.co.uk/f/8/e/f/8efe9ff4209e9ab5e98c62cd39393f0e82.png)
-й справа бит в двоичном представлении
![$n$ $n$](https://dxdy-02.korotkov.co.uk/f/5/5/a/55a049b8f161ae7cfeb0197d75aff96782.png)
.
Пусть
![$a(n)$ $a(n)$](https://dxdy-03.korotkov.co.uk/f/e/b/1/eb1e4d26da404d5c7d56055c19d365a082.png)
это последовательность состоящая из натуральных чисел, такая, что мы начинаем с
![$A:=0$ $A:=0$](https://dxdy-02.korotkov.co.uk/f/9/3/2/932b051bfee310339aae5df53463d05d82.png)
и затем для
![$k=0..\ell(n)$ $k=0..\ell(n)$](https://dxdy-01.korotkov.co.uk/f/4/e/d/4ed5e1cb1f45979b77f72f5db0cb67f282.png)
мы пошагово повторяем следующее:
- Если
, то
; в противном случае
;
.
Тогда
![$a(n)$ $a(n)$](https://dxdy-03.korotkov.co.uk/f/e/b/1/eb1e4d26da404d5c7d56055c19d365a082.png)
это результирующее значение
![$A$ $A$](https://dxdy-02.korotkov.co.uk/f/5/3/d/53d147e7f3fe6e47ee05b88b166bd3f682.png)
.
Например для
![$n=18=10010_2$ $n=18=10010_2$](https://dxdy-01.korotkov.co.uk/f/c/0/e/c0e2be10c3041aa98efad9544455ccfe82.png)
будем иметь:
1.
![$A:=0$ $A:=0$](https://dxdy-02.korotkov.co.uk/f/9/3/2/932b051bfee310339aae5df53463d05d82.png)
;
2.
![$T(n,0)=0$ $T(n,0)=0$](https://dxdy-02.korotkov.co.uk/f/1/6/f/16f51e78e6680903a56f53980c0716c082.png)
,
![$A:=A+1=1$ $A:=A+1=1$](https://dxdy-03.korotkov.co.uk/f/e/e/a/eea604a5671e2d1046099c8b8a00b89982.png)
,
![$A:=A+1=2$ $A:=A+1=2$](https://dxdy-03.korotkov.co.uk/f/2/1/2/212ccd0d485e8ff79e1ec28dca39394182.png)
;
3.
![$T(n,1)=1$ $T(n,1)=1$](https://dxdy-02.korotkov.co.uk/f/1/4/d/14d621e8bf5015d609b737785820ae0e82.png)
,
![$A:=\left\lfloor\frac{A}{2}\right\rfloor=1$ $A:=\left\lfloor\frac{A}{2}\right\rfloor=1$](https://dxdy-01.korotkov.co.uk/f/4/6/7/4670ff09ea6169e406f512deb1b3111f82.png)
,
![$A:=A+1=2$ $A:=A+1=2$](https://dxdy-03.korotkov.co.uk/f/2/1/2/212ccd0d485e8ff79e1ec28dca39394182.png)
;
4.
![$T(n,2)=0$ $T(n,2)=0$](https://dxdy-03.korotkov.co.uk/f/e/c/d/ecd9ac974ae15ba66618d26e1f570f1c82.png)
,
![$A:=A+1=3$ $A:=A+1=3$](https://dxdy-04.korotkov.co.uk/f/3/d/f/3df954934b3304ce221f9810bcac920f82.png)
,
![$A:=A+1=4$ $A:=A+1=4$](https://dxdy-04.korotkov.co.uk/f/3/e/d/3edf9eec136960aa25c09957b299c5f182.png)
;
5.
![$T(n,3)=0$ $T(n,3)=0$](https://dxdy-04.korotkov.co.uk/f/3/7/a/37aab15219eaf1551830fb567fb66fc782.png)
,
![$A:=A+1=5$ $A:=A+1=5$](https://dxdy-04.korotkov.co.uk/f/7/5/9/7598f656082971dc6efc4b3fc6191cb982.png)
,
![$A:=A+1=6$ $A:=A+1=6$](https://dxdy-02.korotkov.co.uk/f/1/d/7/1d7e4c5073c228aa477ac88d85e95fce82.png)
;
6.
![$T(n,4)=1$ $T(n,4)=1$](https://dxdy-03.korotkov.co.uk/f/2/e/6/2e68ff1798223fb570d19b2a2bffafd582.png)
,
![$A:=\left\lfloor\frac{A}{2}\right\rfloor=3$ $A:=\left\lfloor\frac{A}{2}\right\rfloor=3$](https://dxdy-01.korotkov.co.uk/f/8/2/f/82f9c6a2993c7756f76c976124b7534a82.png)
,
![$A:=A+1=4$ $A:=A+1=4$](https://dxdy-04.korotkov.co.uk/f/3/e/d/3edf9eec136960aa25c09957b299c5f182.png)
.
Отсюда
![$a(18)=4$ $a(18)=4$](https://dxdy-04.korotkov.co.uk/f/7/5/d/75d8a8b62dc7d11e80cd20d34bdb4a1f82.png)
.
Пусть
Гипотеза:
если
либо
;
если
либо
;
в противном случае.
Чтобы проверить данную гипотезу вы можете использовать следующую программу на PARI:
Код:
a(n) = my(A=0); for(i=0, logint(n, 2), if(bittest(n, i), A\=2, A++); A++); A
R1(n) = my(v); v=vector(n, i, sum(k=2^(n-1), 2^n-1, a(k)==i))
R(n, k) = if(k==1, 1, if(k<=n, R(n-1, k-1) + R(n-1, 2*(k-1)) + R(n-1, 2*k-1)))
R2(n) = my(v); v=vector(n, i, R(n, i))
test(n) = R1(n)==R2(n)
Существует ли способ это как-то доказать? Возможна ли замкнутая форма для
![$R(n,k)$ $R(n,k)$](https://dxdy-02.korotkov.co.uk/f/5/9/b/59b05fb654c67441dcf92df6bd43572682.png)
?