A356229 Number of maximal gapless submultisets of the prime indices of 2n.
1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 2, 2, 1, 2, 2, 1, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 1, 2, 2, 2, 2, 2, 1, 3, 2, 2, 2, 2, 1, 2, 2, 2, 1, 3, 2, 2, 2, 2, 2, 2, 1, 2, 2, 1, 2, 2, 2, 2, 2, 1, 2, 2, 2, 3, 2, 2, 2, 2, 1, 3, 2, 2, 2, 3, 1, 2, 2, 2, 2, 2, 2, 2, 2, 1
Offset: 1
Keywords
Examples
The prime indices of 2*9282 are {1,1,2,4,6,7}, with maximal gapless submultisets {1,1,2}, {4}, {6,7}, so a(9282) = 3.
Links
- Antti Karttunen, Table of n, a(n) for n = 1..100000
Crossrefs
Programs
-
Mathematica
primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]; Table[Length[Split[primeMS[2n],#1>=#2-1&]],{n,100}]
-
PARI
A287170(n) = { my(f=factor(n)); if(#f~==0, return (0), return(#f~ - sum(i=1, #f~-1, if (primepi(f[i, 1])+1 == primepi(f[i+1, 1]), 1, 0)))); }; A356229(n) = A287170(2*n); \\ Antti Karttunen, Jan 19 2025
Formula
a(n) = A287170(2n).
Extensions
Data section extended to a(105) by Antti Karttunen, Jan 19 2025
Comments