OFFSET
1,2
COMMENTS
There are 3988536 terms below 2*10^9.
Up to a(3988356):
- a(1) is the only odd term.
- a(1) is the only term with 5 consecutive numbers where phi(k) <= phi(k+1).
- the only powers of 2 are a(1)=1 and a(2) = 2.
- of the residues of a(n) mod 10, 4 accounts for greater than 91%.
- if a(n) is divisible by 4, then phi(a(n)) is divisible by 4.
Numbers k such that A057000(k) >= 0 for 3 consecutive terms. - Michel Marcus, Mar 21 2018
EXAMPLE
14 is a term because phi(14) <= phi(15) <= phi(16) <= phi(17).
15 is not a term because phi(15) <= phi(16) <= phi(17) > phi(18).
MAPLE
Phi:= map(numtheory:-phi, [$1..20001]):
DPhi:= Phi[2..-1]-Phi[1..-2]:
C:= select(t -> DPhi[t]>=0, [$1..20000]):
C[select(t -> C[t+2]=C[t]+2, [$1..nops(C)-3])]; # Robert Israel, Mar 26 2018
MATHEMATICA
Drop[#, -2] & /@ Select[SplitBy[#, Last@ # >= 0 &], Length@ # > 2 && #[[1, -1]] >= 0 &][[All, All, 1]] &@ MapIndexed[{First@ #2, #1} &, Differences@ Array[EulerPhi, 2^14]] // Flatten (* Michael De Vlieger, Mar 26 2018 *)
PROG
(PARI) isok(n) = {my(v = vector(4, k, eulerphi(n+k-1))); (v[1] <= v[2]) && (v[2] <= v[3]) && (v[3] <= v[4]); } \\ Michel Marcus, Mar 21 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Torlach Rush, Feb 13 2018
STATUS
approved