# Hofstadter rekursiv
# Plot der Funktion i --> hof(i)

import matplotlib.pyplot as plt
import numpy as np
import pylab as py
import sys
#print("Rekursionstiefe: ",sys.getrecursionlimit())
#z=int(input("gewuenschte Rekursionstiefe: "))
#sys.setrecursionlimit(z)
#print("aktuelle Rekursionstiefe: ",sys.getrecursionlimit())
#print()

def hof(n):
     if n == 1 or n == 2:
          return (1)
     else:
          return hof(n-hof(n-1)) + hof(n - hof(n-2))
          
n=int(input('n =  '))
print()
x = list(range(1, n+1, 1))
y=[]
for i in range(1,n+1):
     y.append(hof(i))
     
print ('hof(',1,') . . . hof(',n,')') 
print()
print(y)

plt.plot(x, y,'bo')     
plt.show()