OFFSET
0,2
COMMENTS
Graphically at large scales this sequence is vaguely self-similar, though in certain sections it acts in a rough manner, in particular in regions surrounding apparent cusps. See the program for a graph to zoom in on these sections.
See program for zoomable graph.
{I assume "program" refers to the Python program? - N. J. A. Sloane, Apr 09 2020)
LINKS
Nathaniel J. Strout, Table of n, a(n) for n = 0..116000
MATHEMATICA
a[0] = 1; a[n_] := a[n] = If[a[n - 1] != 1 && CoprimeQ[n, a[n - 1]], a[n - 1] - 1, a[n - 1] + GCD[a[n - 1], n] + 1]; Array[a, 101, 0] (* Amiram Eldar, Feb 24 2020 *)
nxt[{n_, a_}]:={n+1, If[CoprimeQ[n+1, a]&&a!=1, a-1, a+GCD[a, n+1]+1]}; NestList[nxt, {0, 1}, 70][[;; , 2]] (* Harvey P. Dale, Jun 08 2024 *)
PROG
(Python)
import math
import matplotlib.pyplot as plt
num = 10000
x = []
y = []
# y is the main sequence
def sequence():
a = 1
y.append(a)
for i in range(num):
if (a != 1) and (math.gcd(a, i+1) == 1):
a -= 1
else:
a += math.gcd(a, i+1)+1
x.append(i)
y.append(a)
x.append(num)
sequence()
# code only regarding the plot.
plt.xlim(0, num)
plt.ylim(0, num)
plt.plot(x, y)
plt.xlabel('x - axis')
plt.ylabel('y - axis')
plt.title('Plot of Sequence')
plt.show()
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Nathaniel J. Strout, Feb 24 2020
STATUS
approved