OFFSET
1,1
COMMENTS
These form a subsemigroup and a subsequence of the sequence A289484.
Density: Only 4.3% of the integers between 1 and 400 are doubly A289484.divisible by at least 3 primes. If a term in the sequence is squarefree, it must be divisible by at least 4 primes. If a number n is in the sequence, then every multiple is also in it. Using Wolfram Alpha, about 48% of the integers between 10^40+1 to 10^40+62 were found to be doubly A289484.
EXAMPLE
60=2^2*3*5 is a term because 2^2 > 3 and 2^2*3 > 5.
315=3^2*5*7 is a term because 3^2 > 5 and 3^2*5 > 7.
MAPLE
isA291125 := proc(n)
local pset, p, pprodidx, pprod, nu, falls ;
pset := sort(convert(numtheory[factorset](n), list)) ;
pprod := 1;
falls := 0 ;
for pprodidx from 1 to nops(pset)-1 do
p := pset[pprodidx] ;
nu := padic[ordp](n, p) ;
pprod := pprod*p^nu ;
if pprod > pset[pprodidx+1] then
falls := falls+1 ;
if falls >= 2 then
return true;
end if;
end if;
end do:
return false ;
end proc:
for n from 1 to 3000 do
if isA291125(n) then
printf("%d, ", n) ;
end if;
end do: # R. J. Mathar, Oct 20 2017
PROG
(PARI) is(n, f=factor(n))=my(p=1, t, s); for(i=1, #f~, t=f[i, 1]^f[i, 2]; if(p>t, s++); p*=t); s>1 \\ Charles R Greathouse IV, Jun 10 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Richard Locke Peterson, Aug 17 2017
EXTENSIONS
New name from Charles R Greathouse IV, Jun 11 2020
STATUS
approved