Я прекрасно понимаю, что Вы можете меня послать курить матчасть, но если существует такая возможность, объясните "на пальцах" (на примере), как практически применить Ваши формулы.
Ну, например, на
PARI/GP указанные формулы можно реализовать так:
Код:
? aux(t,m) = sum(i=0,t\(m+1), binomial(t-m*i,i) * (-1)^i * 2^(t-(m+1)*i) )
%1 = (t,m)->sum(i=0,t\(m+1),binomial(t-m*i,i)*(-1)^i*2^(t-(m+1)*i))
? ncomp(n,m) = aux(n,m) - aux(n-1,m)
%2 = (n,m)->aux(n,m)-aux(n-1,m)
? matrix(10,10,n,m,ncomp(n,m))
%3 =
[1 1 1 1 1 1 1 1 1 1]
[1 2 2 2 2 2 2 2 2 2]
[1 3 4 4 4 4 4 4 4 4]
[1 5 7 8 8 8 8 8 8 8]
[1 8 13 15 16 16 16 16 16 16]
[1 13 24 29 31 32 32 32 32 32]
[1 21 44 56 61 63 64 64 64 64]
[1 34 81 108 120 125 127 128 128 128]
[1 55 149 208 236 248 253 255 256 256]
[1 89 274 401 464 492 504 509 511 512]
Как видим, полученные значения для
сходятся с приведёнными в последовательности
A126198. Значит, я нигде не наглюкал.