OFFSET
1,1
COMMENTS
The asymptotic density of this sequence is 1 - 2/zeta(4) + Product_{p prime} (1 - 2/p^4) = 0.001856185541538432217... - Amiram Eldar, Feb 16 2021
Below 9508685764, it suffices to check for n such that either n or n+1 is divisible by p^4 for some p <= 19. - Charles R Greathouse IV, Jul 17 2024
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
80 is a term as 80 and 81 both are divisible by a fourth power, 2^4 and 3^4 respectively.
MATHEMATICA
Select[ Range[2, 25000], Max[ Transpose[ FactorInteger[ # ]] [[2]]] > 3 && Max[ Transpose[ FactorInteger[ # + 1]] [[2]]] > 3 &]
PROG
(PARI) has(n)=vecmax(factor(n)[, 2])>3
is(n)=has(n+1)&&has(n) \\ Charles R Greathouse IV, Dec 19 2018
(PARI) list(lim)=my(v=List(), x=1); forfactored(n=81, lim\1+1, if(vecmax(n[2][, 2])>3, if(x, listput(v, n[1]-1), x=1), x=0)); Vec(v) \\ Charles R Greathouse IV, Dec 19 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Mar 04 2002
EXTENSIONS
a(0) = 0 removed by Charles R Greathouse IV, Dec 19 2018
STATUS
approved