OFFSET
1,2
COMMENTS
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 the prime shadow A181819(n) to be the product of primes indexed by the exponents in the prime factorization of n. For example, 90 = prime(1)*prime(2)^2*prime(3) has prime shadow prime(1)*prime(2)*prime(1) = 12.
EXAMPLE
The terms together with their prime indices begin:
1: {} 78: {1,2,6} 158: {1,22}
3: {2} 83: {23} 165: {2,3,5}
5: {3} 85: {3,7} 174: {1,2,10}
11: {5} 86: {1,14} 177: {2,17}
15: {2,3} 93: {2,11} 179: {41}
17: {7} 94: {1,15} 187: {5,7}
26: {1,6} 109: {29} 191: {43}
31: {11} 123: {2,13} 196: {1,1,4,4}
33: {2,5} 126: {1,2,2,4} 201: {2,19}
41: {13} 127: {31} 202: {1,26}
51: {2,7} 130: {1,3,6} 205: {3,13}
55: {3,5} 146: {1,21} 211: {47}
58: {1,10} 148: {1,1,12} 241: {53}
59: {17} 155: {3,11} 244: {1,1,18}
67: {19} 157: {37} 249: {2,23}
For example, 126 is in the sequence because its prime indices {1,2,2,4} have shadows {1,2,2,3}, with product 12, which is also the prime shadow of 126.
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
red[n_]:=If[n==1, 1, Times@@Prime/@Last/@FactorInteger[n]];
Select[Range[100], Times@@red/@primeMS[#]==red[#]&]
CROSSREFS
The prime terms are A006450.
This is a ranking of the partitions counted by A353396.
A003963 gives product of prime indices.
A324850 lists numbers divisible by the product of their prime indices.
Numbers divisible by their prime shadow:
- counted by A325702
- listed by A325755
- co-recursive version A325756
- nonprime recursive version A353389
KEYWORD
nonn
AUTHOR
Gus Wiseman, May 17 2022
STATUS
approved