OFFSET
1,1
COMMENTS
Numbers n such that A324848(n) = 1.
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.
If k is in A324846, then k*prime(k) is in the sequence. - Robert Israel, Mar 22 2019
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The sequence of terms together with their prime indices begins:
2: {1}
10: {1,3}
14: {1,4}
15: {2,3}
22: {1,5}
26: {1,6}
34: {1,7}
38: {1,8}
45: {2,2,3}
46: {1,9}
50: {1,3,3}
55: {3,5}
58: {1,10}
62: {1,11}
70: {1,3,4}
74: {1,12}
82: {1,13}
86: {1,14}
94: {1,15}
98: {1,4,4}
MAPLE
filter:= proc(n) local F;
F:= select(t -> n mod numtheory:-pi(t[1])=0, ifactors(n)[2]);
nops(F)=1 and F[1][2]=1
end proc:
select(filter, [$2..1000]); # Robert Israel, Mar 22 2019
MATHEMATICA
Select[Range[100], Total[Cases[If[#==1, {}, FactorInteger[#]], {p_, k_}:>k/; Divisible[#, PrimePi[p]]]]==1&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 21 2019
STATUS
approved