(Python)

def main():
    pass

if __name__ == '__main__':
    main()
figs = {}
max = 18000000000000
current = 1
increment = 1000000
row = 6
col = 3
inc = 3
colinc = 3
name = "m_densities_big"
fnum = 0
total = 0
file = open(str(name)+str(fnum)+".txt","w")
while current <= max:
   if row > current-increment:
      value = row
      colinc = 3
   else:
      value = current-increment-(current-increment-row)%col+col
      colinc = int((current-increment-row)/col)+4
   while value<=current:
      if value in figs:
         figs[value] += [inc,colinc]
      else:
         figs[value] = [inc,colinc]
      colinc += 1
      value += col
   col = row
   inc += 1
   row += inc
   if row >current:
      for key in sorted(figs.keys()):
       total += len(figs[key])/2
       mult = int(len(figs[key])/2)
       if len(figs[key])/2 >=1:
         file.write (str(key)+" "+str(mult)+" "+str(figs[key])+"\n")
      file.close()
      fnum += 1
      file = open(str(name)+str(fnum)+".txt","w")
      row = 6
      col = 3
      inc = 3
      current += increment
      figs = {}