login
Nonprimes n with gpf(n) > gpf(n+1), where gpf=A006530 (greatest prime factor).
2

%I #11 Jul 03 2018 21:22:01

%S 14,15,26,34,35,38,39,44,49,51,55,62,63,65,69,74,76,80,86,87,94,95,99,

%T 104,111,116,118,119,122,123,124,129,134,142,143,146,152,153,155,158,

%U 159,161,164,174,183,185,186,188,194,195,202,203,206,207,209,214,215

%N Nonprimes n with gpf(n) > gpf(n+1), where gpf=A006530 (greatest prime factor).

%C Subsequence of A070087.

%H Robert Israel, <a href="/A087430/b087430.txt">Table of n, a(n) for n = 1..10000</a>

%p N:= 1000: # to get all terms < N

%p V:= Vector(N):

%p p:= 1:

%p do

%p p:= nextprime(p);

%p if p > N then break fi;

%p V[[seq(k*p,k=1..N/p)]]:= p

%p od:

%p select(t -> V[t] > V[t+1] and not isprime(t), [$1..N-1]); # _Robert Israel_, Jul 03 2018

%t With[{nn=250},Select[Complement[Range[nn],Prime[Range[PrimePi[ nn]]]], FactorInteger[ #][[-1,1]]>FactorInteger[#+1][[-1,1]]&]] (* _Harvey P. Dale_, Jul 14 2016 *)

%Y Cf. A087429.

%K nonn

%O 1,1

%A _Reinhard Zumkeller_, Sep 02 2003