#!/usr/bin/python

import sys, math

test_notes = int(sys.argv[1]) or 49
final_notes = int(sys.argv[2]) or 64

basehz = 3579545/32.0

results = []

for a440 in range(4300,4500):
    error = 0
    for note in range(1,test_notes):
        notehz = a440 / 10.0 * math.pow(2.0, (note - 49) / 12.0);
        period = int(round(basehz / notehz))
        tonehz = basehz / period
        error += abs(notehz-tonehz)
        #print a440,note,notehz,notehz-tonehz,period
        #if a440 == 4405:
        #  print '%d,%f,%d' % (note,tonehz-notehz,period)
    results.append((error, a440))

results.sort()
best_error, best_a440 = results[0]
best_a440 /= 10.0
print '//', best_a440, best_error, test_notes

print "const int note_table[%d] = {" % final_notes
for note in range(0,final_notes):
    notehz = best_a440 * math.pow(2.0, (note - 49) / 12.0);
    period = int(round(basehz / notehz))
    print '%d,' % period,
print "};"