A335549 - OEIS

A335549

Number of normal patterns matched by the multiset of prime indices of n in weakly increasing order.

1, 2, 2, 3, 2, 3, 2, 4, 3, 3, 2, 5, 2, 3, 3, 5, 2, 5, 2, 5, 3, 3, 2, 7, 3, 3, 4, 5, 2, 4, 2, 6, 3, 3, 3, 7, 2, 3, 3, 7, 2, 4, 2, 5, 5, 3, 2, 9, 3, 5, 3, 5, 2, 7, 3, 7, 3, 3, 2, 7, 2, 3, 5, 7, 3, 4, 2, 5, 3, 4, 2, 10, 2, 3, 5, 5, 3, 4, 2, 9, 5, 3, 2, 7, 3, 3, 3

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

COMMENTS

First differs from A181796 at a(90) = 8 A181796(90) = 7.

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

We define a (normal) pattern to be a finite sequence covering an initial interval of positive integers. Patterns are counted by A000670 and ranked by A333217. A sequence S is said to match a pattern P if there is a not necessarily contiguous subsequence of S whose parts have the same relative order as P. For example, (3,1,1,3) matches (1,1,2), (2,1,1), and (2,1,2), but avoids (1,2,1), (1,2,2), and (2,2,1).

LINKS

Table of n, a(n) for n=1..87.

Wikipedia, Permutation pattern

Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.

EXAMPLE

The Heinz number of (1,2,2,3) is 90 and it matches 8 patterns: (), (1), (11), (12), (112), (122), (123), (1223); so a(90) = 8.

MATHEMATICA

primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];

mstype[q_]:=q/.Table[Union[q][[i]]->i, {i, Length[Union[q]]}];

Table[Length[Union[mstype/@Subsets[primeMS[n]]]], {n, 100}]

CROSSREFS

The version for standard compositions instead of prime indices is A335454.