A387425 Partial sums of A387413.

0, 1, 2, 2, 4, 5, 7, 10, 13, 13, 16, 19, 22, 25, 28, 31, 33, 33, 36, 40, 44, 47, 51, 55, 56, 58, 61, 64, 66, 69, 73, 78, 79, 84, 87, 92, 96, 99, 103, 106, 108, 110, 113, 118, 121, 124, 129, 131, 135, 140, 145, 150, 154, 159, 164, 169, 169, 174, 177, 182, 187, 191, 196, 201, 207, 213, 216, 221, 225, 231, 235, 240, 243

Offset

1,3

Links

Antti Karttunen, Table of n, a(n) for n = 1..65537

Index entries for sequences related to binary expansion of n.

Index entries for sequences related to prime indices in the factorization of n.

Prog

(PARI)
up_to = 65537;
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A387413(n) = { my(a=binary(n), b=binary(A003961(n)), i=1); while(i<=#a, if(a[i]!=b[i], return(#a-(i-1))); i++); (0); };
A387425list(up_to) = { my(v=vector(up_to)); v[1] = 0; for(n=2, up_to, v[n] = v[n-1]+A387413(n)); (v); };
v387425 = A387425list(up_to);
A387425(n) = v387425[n];

Crossrefs

Cf. A003961, A387413.

Cf. also A387424.

Sequence in context: A026930 A211862 A286948 * A211863 A333190 A098859

Adjacent sequences: A387422 A387423 A387424 * A387426 A387427 A387428

Keyword

nonn

Author

Antti Karttunen, Sep 01 2025

Status

approved