A381836 k/25 is in this list if A053824(k) < A112765(k), i.e. if digitsum(k, 5) < valuation(k, 5).
1, 5, 10, 25, 30, 50, 75, 125, 130, 150, 175, 250, 275, 375, 500, 625, 630, 650, 675, 750, 775, 875, 1000, 1250, 1275, 1375, 1500, 1875, 2000, 2500, 3125, 3130, 3150, 3175, 3250, 3275, 3375, 3500, 3750, 3775, 3875, 4000, 4375, 4500, 5000, 5625, 6250, 6275, 6375
Offset: 1
Programs
-
Maple
aList := upto -> local k; [seq(k/25, k in select(n -> add(convert(n, base, 5)) < padic[ordp](n, 5), [seq(25..upto,25)]))]: aList(160000);
-
Mathematica
Select[Range[160000],DigitSum[#,5]
Stefano Spezia, Mar 08 2025 *) -
SageMath
def aList(upto, b): return [n/b^2 for n in srange(b^2, upto, b^2) if sum(n.digits(b)) < valuation(n, b)] print(aList(160000, 5))