# maxruntest 
# Fig. 38.2 
# input: 1024

import random, statistics

NUM_TRIALS = 100

n = int(input("n=? "))
maxRuns = []
for j in range(NUM_TRIALS): 
  maxRun = 0
  count = 0
  while count < n: 
    i = 0
    while (random.choice(['H','T']) == 'H') and (count < n): 
      i += 1
      count += 1
    if i > maxRun:
      maxRun = i
  print(f"{maxRun:2}", end=' ')
  if (j+1)%20 == 1:
    print()
  maxRuns.append(maxRun)
print()
print(f"Mean maxrun = {statistics.mean(maxRuns):.2f}")
print(f"Stdev = {statistics.stdev(maxRuns):.2f}")
