OFFSET
1,2
COMMENTS
a(p) = 2^p - 1 for prime p.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..1000
EXAMPLE
a(4) = 12, because 12 is the smallest number with a binary run length product of 4.
12 decimal = 1100 binary. Run lengths in binary are 2,2, and 2*2 = 4.
PROG
(PARI)
a(n)={
my(p=if(n==1, [], my(f=factor(n)); concat(vector(#f~, i, f[i, 1]*vector(f[i, 2], j, 1)))));
my(i=1, j=#p, b=0);
while(i<=j, if(bittest(b, 0),
if(p[j]>3||j==i||p[i+1]!=2, b<<=p[j]; j--, b<<=4; i+=2),
b++; b<<=p[i]; b--; i++));
b
} \\ Andrew Howroyd, Nov 09 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Weimholt, Nov 05 2009
EXTENSIONS
a(31) onward from Andrew Howroyd, Nov 09 2025
STATUS
approved
