login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A319878
Odd numbers whose product of prime indices (A003963) is a square of a squarefree number (A062503).
1
1, 7, 9, 23, 25, 97, 121, 151, 161, 169, 175, 183, 185, 195, 207, 225, 227, 289, 541, 661, 679, 687, 781, 841, 847, 873, 957, 961, 1009, 1089, 1193, 1427, 1563, 1589, 1681, 1819, 1849, 1879, 1895, 2023, 2043, 2167, 2193, 2209, 2231, 2425, 2437, 2585, 2601
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. The multiset multisystem with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}. This sequence lists all MM-numbers of 2-regular (all vertex-degrees are 2) multiset partitions (no empty parts).
EXAMPLE
The sequence of multiset partitions whose MM-numbers belong to the sequence begins:
1: {}
7: {{1,1}}
9: {{1},{1}}
23: {{2,2}}
25: {{2},{2}}
97: {{3,3}}
121: {{3},{3}}
151: {{1,1,2,2}}
161: {{1,1},{2,2}}
169: {{1,2},{1,2}}
175: {{2},{2},{1,1}}
183: {{1},{1,2,2}}
185: {{2},{1,1,2}}
195: {{1},{2},{1,2}}
207: {{1},{1},{2,2}}
225: {{1},{1},{2},{2}}
227: {{4,4}}
289: {{4},{4}}
541: {{1,1,3,3}}
661: {{5,5}}
679: {{1,1},{3,3}}
687: {{1},{1,3,3}}
781: {{3},{1,1,3}}
841: {{1,3},{1,3}}
847: {{1,1},{3},{3}}
873: {{1},{1},{3,3}}
957: {{1},{3},{1,3}}
961: {{5},{5}}
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[1, 100, 2], Or[#==1, SameQ[##, 2]&@@Last/@FactorInteger[Times@@primeMS[#]]]&]
KEYWORD
nonn
AUTHOR
Gus Wiseman, Dec 17 2018
STATUS
approved