OFFSET
1,2
COMMENTS
FORMULA
a(n) = A129597(n)/2 for n > 1.
EXAMPLE
The terms together with their prime indices begin:
1: {} 31: {11} 71: {20}
2: {1} 32: {1,1,1,1,1} 73: {21}
3: {2} 37: {12} 76: {1,1,8}
5: {3} 41: {13} 79: {22}
7: {4} 43: {14} 80: {1,1,1,1,3}
8: {1,1,1} 44: {1,1,5} 83: {23}
11: {5} 45: {2,2,3} 89: {24}
12: {1,1,2} 47: {15} 92: {1,1,9}
13: {6} 48: {1,1,1,1,2} 97: {25}
17: {7} 52: {1,1,6} 99: {2,2,5}
19: {8} 53: {16} 101: {26}
20: {1,1,3} 59: {17} 103: {27}
23: {9} 61: {18} 107: {28}
27: {2,2,2} 63: {2,2,4} 108: {1,1,2,2,2}
28: {1,1,4} 67: {19} 109: {29}
29: {10} 68: {1,1,7} 112: {1,1,1,1,4}
MAPLE
filter:= proc(n) issqr(n/max(numtheory:-factorset(n))) end proc:
filter(1):= true:
select(filter, [$1..200]); # Robert Israel, Nov 26 2022
MATHEMATICA
sqrQ[n_]:=IntegerQ[Sqrt[n]];
Select[Range[100], sqrQ[#*FactorInteger[#][[-1, 1]]]&]
PROG
(PARI) isok(m) = (m==1) || issquare(m/vecmax(factor(m)[, 1])); \\ Michel Marcus, Aug 12 2021
CROSSREFS
Removing 1 gives a subset of A026424.
The unsorted even version is A129597.
Except the first term, the even version is 2*a(n).
A000290 lists squares.
A001221 counts distinct prime factors.
A001222 counts all prime factors.
A006530 gives the greatest prime factor.
A061395 gives the greatest prime index.
A027193 counts partitions of odd length.
A344606 counts alternating permutations of prime indices.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 10 2021
STATUS
approved