cem 0.1.4 → 0.1.5

data/examples/aoc2018/day24.rb

@@ -0,0 +1,156 @@
+#
+# https://adventofcode.com/2018/day/24 part 1 and 2
+#
+# Does not use any cem functions
+#
+
+Group = Struct.new("Group", :immuneSystem, :units, :hitPoints, :weakness, :immunities, :damage, :damageType, :initiative) {
+
+  def damageCalc(other)
+    damage * units * (other.immunities.include?(damageType) ? 0 : 1) * (other.weakness.include?(damageType) ? 2 : 1)
+  end
+  
+  def effectivePow
+    damage * units
+  end
+}
+
+lines = File.readlines("inputs/day24_input.txt")
+
+groups = []
+immuneSystem = false
+
+lines.each { |line|
+  line.strip!
+
+  if line =~ /^Immune System:$/
+    immuneSystem = true
+  elsif line =~ /^Infection:$/
+    immuneSystem = false
+  elsif line =~ /^([+-]?\d+) units each with ([+-]?\d+) hit points\s*(.*?) with an attack that does ([+-]?\d+) (\w+) damage at initiative ([+-]?\d+)$/
+    units = $1.to_i
+    hitPoints = $2.to_i
+    immunities = []
+    weakness = []
+    w = $3
+    damage = $4.to_i
+    damageType = $5
+    initiative = $6.to_i
+    
+    w.gsub(/\(|\)/, "").split("; ").each { |s|
+      if s =~ /immune to (.*)/
+        immunities = $1.split(", ")
+      end
+      if s =~ /weak to (.*)/
+        weakness = $1.split(", ")
+      end
+    }
+    
+    g = Group.new(immuneSystem, units, hitPoints, weakness, immunities, damage, damageType, initiative)
+    groups << g
+    # puts g.inspect
+  elsif line == ""
+  
+  else
+    raise line.inspect
+  end    
+}
+
+def run(groups)
+
+  max_init = groups.map { |g| g.initiative }.max + 1
+
+  @max_init = max_init
+
+  round = 0
+  
+  while groups.count { |g| g.immuneSystem } != 0 && groups.count { |g| !g.immuneSystem } != 0
+
+    anyDamage = false  
+  
+    groups.sort_by! { |g| [-g.damage * g.units, -g.initiative] }
+    
+    # Target Selection
+    attackedGroups = []
+    attackPairs = []
+    groups.each { |g|
+      others = groups.select { |other| g.immuneSystem != other.immuneSystem && !attackedGroups.include?(other) }
+      
+      groupToAttack = others.sort_by { |other| [-g.damageCalc(other), -other.damage * other.units, -other.initiative] }.first
+      
+      if groupToAttack && g.damageCalc(groupToAttack) != 0
+        attackedGroups << groupToAttack
+        attackPairs << [g, groupToAttack]
+      end
+    }
+    
+    attackPairs.sort_by! { |x| -x[0].initiative } 
+    
+    # Attack phase
+    attackPairs.each { |x|
+      from, to = x
+      damage = from.damageCalc(to)
+      
+      killedUnits = damage / to.hitPoints
+      if killedUnits > 0
+        anyDamage = true
+      end
+      to.units -= killedUnits 
+          
+      if to.units <= 0
+        to.units = 0
+        groups = groups - [to]
+      end
+    }
+
+    groups.sort_by! { |g| -(g.immuneSystem ? 1 : 0) * 100000000000 - (g.damage * g.units * max_init + g.initiative)}
+    
+    # puts round
+    round += 1
+    # system("cls") if round % 1024 == 0    
+    # groups.each { |g| puts g.inspect } if round % 1024 == 0
+      
+    if !anyDamage 
+      groups.each { |g| puts g.inspect }
+      return 
+    end
+  end
+  
+  groups
+end  
+
+tempGroup = groups.map { |g| g.dup }
+
+puts "Part 1: #{run(tempGroup).sum{|g| g.units }}"
+
+def part2(groups)
+
+  minAdd = 1
+  
+  loop do  
+    
+    tempGroup = groups.map { |g| g.dup }
+  
+    tempGroup.each { |g| 
+      if g.immuneSystem
+        g.damage += minAdd
+      end
+    }
+  
+    res = run(tempGroup)
+    
+    if res == nil
+      puts "Boost: #{minAdd} - Got stuck"
+    else
+      puts "Boost: #{minAdd} - #{res[0].immuneSystem}"  
+    
+      if res[0].immuneSystem
+        return res
+      end
+    end
+  
+    minAdd += 1
+  end 
+end
+
+puts "Part 2: #{part2(groups).sum { |g| g.units }}"

data/examples/aoc2018/day25.rb

@@ -0,0 +1,52 @@
+
+#
+# https://adventofcode.com/2018/day/25 part 1 (there is only 1
+#
+# Uses Point4D, but not much to see here.
+#
+require 'cem'
+
+lines = File.readlines("inputs/day25_input.txt")
+
+input = []
+
+lines.each { |line|
+
+  if line =~ /^\s*([+-]?\d+),([+-]?\d+),([+-]?\d+),([+-]?\d+)\s*$/
+    v1 = $1.to_i
+    v2 = $2.to_i
+    v3 = $3.to_i
+    v4 = $4.to_i
+    input << Point4D.new(v1,v2,v3,v4)
+  else
+    raise line.inspect
+  end
+}
+
+constellations = []
+
+input.each { |p|
+  if c = constellations.find { |c| c.any? { |p2| p2.manhattan(p) <= 3 } }
+    c << p
+  else
+    constellations << [p]
+  end
+}
+
+loop do  
+  count = constellations.count { |c| c.size > 0 }
+  
+  constellations.each_with_index { |c,i|
+    constellations.each_with_index { |d,j|
+      next if j <= i 
+      if c.any? { |p| d.any? { |p2| p2.manhattan(p) <= 3 } }
+        c.push(*d)
+        d.clear
+      end
+    }
+  }
+  break if count == constellations.count { |c| c.size > 0 }
+end
+
+puts "Part 1: #{constellations.count { |c| c.size > 0 }}"
+

data/examples/aoc2018/day9.rb

@@ -0,0 +1,51 @@
+
+#
+# https://adventofcode.com/2018/day/9 part 1 and part 2
+#
+# Does not use any cem functions \o/, but did you know that 
+#  - Array#shift is supposed to be fast, but is not fast enough for this problem?
+# 
+
+players = 428 
+marbles = 70825
+
+def day9(players, marbles)
+
+  scores = [0] * players
+
+  nexta = 3
+  storage = [0, 2, 1] 
+  current = 1
+  front = 0
+
+  while nexta < marbles 
+
+    if nexta % 23 == 0
+      
+      current = (current + storage.size - 7) % storage.size
+      worth = nexta + storage.delete_at(current)
+      scores[nexta % players] += worth
+    
+    else
+    
+      # Avoid wrapping the index
+      while current + 2 > storage.size
+        front_value = storage[front]
+        front += 1
+        storage.push front_value  
+      end
+      
+      current += 2
+      storage.insert(current, nexta)
+      
+    end
+    
+    nexta += 1
+    
+  end
+  
+  return scores.max
+end  
+  
+puts "One Star : #{day9(players, marbles)}"
+puts "Two Stars: #{day9(players, marbles * 100)}"

data/examples/aoc2018/day9_circular.rb

@@ -0,0 +1,125 @@
+
+#
+# https://adventofcode.com/2018/day/9 part 1 and part 2
+#
+# Does not use any cem functions \o/ 
+# 
+# Unfinished attempt to develop a circular array
+# 
+
+players = 428 
+marbles = 7082500
+
+#players = 10 
+#marbles = 1618
+
+scores = [0] * 428
+
+nexta = 3
+lasta = 0
+storage = [1, 2, 1] 
+current = 1
+
+#
+#
+# Caveats:
+#   - When writing to indices greater than @size, the CircularArray does not expand, but the index wraps.
+#
+class CircularArray 
+
+  attr_accessor :data
+  attr_reader  :size, :start
+
+  def initialize(initialSize)
+	@data = Array.new(initialSize)
+	@size = 0
+  end
+  
+  def index(index)
+    return (@start + index) % @size
+  end
+  
+  def [](i)
+    return @data[index(i)]
+  end
+  
+  def []=(i, newValue)
+    @data[index(i)] = newValue
+  end
+  
+  def each
+    (0..@size).each do |x|
+      yield self[x]
+    end
+  end
+  
+  def each_with_index
+    (0..@size).each do |x|
+      yield self[x], x
+    end
+  end
+  
+  def <<(newValue)
+    if @size == @data.size
+      @data.insert(index(@size) + 1, [nil] * @size)
+    end	
+    
+	@data[index(@size)] = newValue
+	@size += 1	
+  end
+  
+  
+end
+
+
+while nexta < marbles 
+
+  if nexta / 1000 != lasta
+    lasta = nexta / 1000
+	puts lasta * 1000
+  end
+  
+  n = storage.size
+  if nexta % 23 == 0
+  
+    current = (current + n - 7) % n
+    worth = nexta + storage.delete_at(current)
+    scores[nexta % players] += worth
+    
+    if worth == marbles
+      break;
+    end
+  
+  else
+  
+    if nexta % 23 == 1 && current + 46 < n && nexta + 22 < marbles
+    
+      current += 2
+      
+      storage.insert(current, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
+      (0..21).each { |i|
+      
+        # puts storage[current..current + 40].inspect
+      
+        storage[current + i * 2] = nexta
+        nexta += 1
+        storage[current + i * 2 + 1] = storage[current + 20 + i]
+        
+
+      }
+      current += 42
+      nexta -= 1
+    
+    else
+      current = (current + 2) % n
+      storage.insert(current, nexta)
+    end
+    
+  end
+  
+  nexta += 1
+  
+end
+  
+  
+puts scores.max

data/examples/aoc2018/inputs/day10_input.txt

@@ -0,0 +1,395 @@
+position=<-41214, -10223> velocity=< 4,  1>
+position=< 41635, -10215> velocity=<-4,  1>
+position=<-51569,  20845> velocity=< 5, -2>
+position=<-30826,  20850> velocity=< 3, -2>
+position=< 31249, -30930> velocity=<-3,  3>
+position=<-30862, -20578> velocity=< 3,  2>
+position=<-20528, -30930> velocity=< 2,  3>
+position=< 52013, -51642> velocity=<-5,  5>
+position=<-20491,  41552> velocity=< 2, -4>
+position=<-41218,  31200> velocity=< 4, -3>
+position=< 31271, -10216> velocity=<-3,  1>
+position=< 52010,  51915> velocity=<-5, -5>
+position=< 41634, -30933> velocity=<-4,  3>
+position=<-30826,  41551> velocity=< 3, -4>
+position=<-20499,  31199> velocity=< 2, -3>
+position=< 20953,  41559> velocity=<-2, -4>
+position=<-51576, -41285> velocity=< 5,  4>
+position=<-20523,  51906> velocity=< 2, -5>
+position=< 10593,  51914> velocity=<-1, -5>
+position=<-51593,  31204> velocity=< 5, -3>
+position=< 20893, -30930> velocity=<-2,  3>
+position=< 20929,  20843> velocity=<-2, -2>
+position=< 51985,  20845> velocity=<-5, -2>
+position=<-51534, -10215> velocity=< 5,  1>
+position=<-41222,  51908> velocity=< 4, -5>
+position=<-30859,  51914> velocity=< 3, -5>
+position=<-30874,  51906> velocity=< 3, -5>
+position=<-10129, -10220> velocity=< 1,  1>
+position=<-41198,  31203> velocity=< 4, -3>
+position=<-41238,  20847> velocity=< 4, -2>
+position=<-10141,  20850> velocity=< 1, -2>
+position=< 41661,  51915> velocity=<-4, -5>
+position=<-41230,  31202> velocity=< 4, -3>
+position=< 10561, -20578> velocity=<-1,  2>
+position=<-30835,  20848> velocity=< 3, -2>
+position=<-20467, -30926> velocity=< 2,  3>
+position=<-51561,  31202> velocity=< 5, -3>
+position=< 20900, -41280> velocity=<-2,  4>
+position=< 10588,  31198> velocity=<-1, -3>
+position=<-30854,  41551> velocity=< 3, -4>
+position=< 31263,  10488> velocity=<-3, -1>
+position=<-10117,  10490> velocity=< 1, -1>
+position=<-51569,  20843> velocity=< 5, -2>
+position=< 20928, -51640> velocity=<-2,  5>
+position=< 10537, -10221> velocity=<-1,  1>
+position=<-41201, -41286> velocity=< 4,  4>
+position=<-10141, -20572> velocity=< 1,  2>
+position=< 20908, -30929> velocity=<-2,  3>
+position=<-20512,  51911> velocity=< 2, -5>
+position=<-30835, -10224> velocity=< 3,  1>
+position=< 41610, -10219> velocity=<-4,  1>
+position=< 10537,  20849> velocity=<-1, -2>
+position=< 41644, -20579> velocity=<-4,  2>
+position=< 51990, -41289> velocity=<-5,  4>
+position=< 10569, -20571> velocity=<-1,  2>
+position=< 20905,  41559> velocity=<-2, -4>
+position=< 20940,  20843> velocity=<-2, -2>
+position=<-20512,  20843> velocity=< 2, -2>
+position=<-41202,  51915> velocity=< 4, -5>
+position=<-10173,  51911> velocity=< 1, -5>
+position=<-41230, -10220> velocity=< 4,  1>
+position=< 31279,  31205> velocity=<-3, -3>
+position=<-41226,  10492> velocity=< 4, -1>
+position=<-51577, -41280> velocity=< 5,  4>
+position=<-10125,  20846> velocity=< 1, -2>
+position=< 41658,  20845> velocity=<-4, -2>
+position=< 31274,  41555> velocity=<-3, -4>
+position=<-10117,  51913> velocity=< 1, -5>
+position=< 31287,  51910> velocity=<-3, -5>
+position=< 10545, -20576> velocity=<-1,  2>
+position=<-20520,  51913> velocity=< 2, -5>
+position=<-20528, -10216> velocity=< 2,  1>
+position=<-41180, -20579> velocity=< 4,  2>
+position=< 10561, -41283> velocity=<-1,  4>
+position=< 31252,  51909> velocity=<-3, -5>
+position=<-30822, -30927> velocity=< 3,  3>
+position=< 20929, -41281> velocity=<-2,  4>
+position=<-51553,  10489> velocity=< 5, -1>
+position=<-30883, -41287> velocity=< 3,  4>
+position=<-30846,  10487> velocity=< 3, -1>
+position=< 10594,  10495> velocity=<-1, -1>
+position=< 10547,  31200> velocity=<-1, -3>
+position=< 41658, -10216> velocity=<-4,  1>
+position=<-51585,  41553> velocity=< 5, -4>
+position=<-41198, -20579> velocity=< 4,  2>
+position=<-10168, -10223> velocity=< 1,  1>
+position=<-30851, -20575> velocity=< 3,  2>
+position=< 10574, -41283> velocity=<-1,  4>
+position=< 31260,  20849> velocity=<-3, -2>
+position=< 41607, -51636> velocity=<-4,  5>
+position=< 20944,  10494> velocity=<-2, -1>
+position=< 31255, -10223> velocity=<-3,  1>
+position=<-20510,  51906> velocity=< 2, -5>
+position=<-30851, -41286> velocity=< 3,  4>
+position=<-41198, -51636> velocity=< 4,  5>
+position=<-41190,  41553> velocity=< 4, -4>
+position=<-10170, -30930> velocity=< 1,  3>
+position=<-41233,  20848> velocity=< 4, -2>
+position=< 31255,  31197> velocity=<-3, -3>
+position=<-30843,  10488> velocity=< 3, -1>
+position=< 20948,  10492> velocity=<-2, -1>
+position=<-41206, -51639> velocity=< 4,  5>
+position=<-10161, -30928> velocity=< 1,  3>
+position=< 31308,  20848> velocity=<-3, -2>
+position=<-10149, -30934> velocity=< 1,  3>
+position=<-10165,  20841> velocity=< 1, -2>
+position=< 20948, -51642> velocity=<-2,  5>
+position=<-41222, -51637> velocity=< 4,  5>
+position=<-30823,  51914> velocity=< 3, -5>
+position=< 41642, -20574> velocity=<-4,  2>
+position=<-20523, -41280> velocity=< 2,  4>
+position=< 10589,  41559> velocity=<-1, -4>
+position=< 10577, -41283> velocity=<-1,  4>
+position=< 20900, -30934> velocity=<-2,  3>
+position=<-30851, -51635> velocity=< 3,  5>
+position=< 31249,  10490> velocity=<-3, -1>
+position=< 51992, -41280> velocity=<-5,  4>
+position=<-30833, -20576> velocity=< 3,  2>
+position=<-20509, -41285> velocity=< 2,  4>
+position=< 41607,  51907> velocity=<-4, -5>
+position=<-20485, -51639> velocity=< 2,  5>
+position=<-20523,  31196> velocity=< 2, -3>
+position=< 10585,  41558> velocity=<-1, -4>
+position=< 41653, -30927> velocity=<-4,  3>
+position=< 20900, -41285> velocity=<-2,  4>
+position=< 31276, -20573> velocity=<-3,  2>
+position=<-10117,  41555> velocity=< 1, -4>
+position=<-10133,  51910> velocity=< 1, -5>
+position=< 31249,  10490> velocity=<-3, -1>
+position=< 31295,  10492> velocity=<-3, -1>
+position=<-30843, -20570> velocity=< 3,  2>
+position=<-30827, -20572> velocity=< 3,  2>
+position=< 10586, -51639> velocity=<-1,  5>
+position=< 10589, -10223> velocity=<-1,  1>
+position=< 20943,  20848> velocity=<-2, -2>
+position=<-10157, -10218> velocity=< 1,  1>
+position=< 10598, -41284> velocity=<-1,  4>
+position=< 41660, -30934> velocity=<-4,  3>
+position=< 31263,  20849> velocity=<-3, -2>
+position=< 31276,  20846> velocity=<-3, -2>
+position=<-41238,  20846> velocity=< 4, -2>
+position=<-10132, -41285> velocity=< 1,  4>
+position=< 20924,  51906> velocity=<-2, -5>
+position=<-20504,  51908> velocity=< 2, -5>
+position=< 52007, -20576> velocity=<-5,  2>
+position=< 51957,  41555> velocity=<-5, -4>
+position=< 31291, -10224> velocity=<-3,  1>
+position=<-10161, -10217> velocity=< 1,  1>
+position=<-30847, -20575> velocity=< 3,  2>
+position=< 41615, -41287> velocity=<-4,  4>
+position=< 41626,  20850> velocity=<-4, -2>
+position=< 31307, -10224> velocity=<-3,  1>
+position=<-20470,  51915> velocity=< 2, -5>
+position=<-51593,  31198> velocity=< 5, -3>
+position=<-30873,  31196> velocity=< 3, -3>
+position=<-20504, -10220> velocity=< 2,  1>
+position=<-10157, -20571> velocity=< 1,  2>
+position=<-10156,  41555> velocity=< 1, -4>
+position=< 51970, -30933> velocity=<-5,  3>
+position=< 41626,  41557> velocity=<-4, -4>
+position=< 52009, -20578> velocity=<-5,  2>
+position=< 51994, -20577> velocity=<-5,  2>
+position=< 41626,  10491> velocity=<-4, -1>
+position=<-41193,  10494> velocity=< 4, -1>
+position=<-51545, -51635> velocity=< 5,  5>
+position=< 31279, -20571> velocity=<-3,  2>
+position=< 41611,  31200> velocity=<-4, -3>
+position=<-41182,  41552> velocity=< 4, -4>
+position=< 52008,  10493> velocity=<-5, -1>
+position=< 31273,  51910> velocity=<-3, -5>
+position=<-20512, -10215> velocity=< 2,  1>
+position=<-30843, -51643> velocity=< 3,  5>
+position=<-30849,  10486> velocity=< 3, -1>
+position=<-20502,  10490> velocity=< 2, -1>
+position=<-20496, -51640> velocity=< 2,  5>
+position=< 41615,  10489> velocity=<-4, -1>
+position=<-20485, -41284> velocity=< 2,  4>
+position=<-20520,  20849> velocity=< 2, -2>
+position=< 51986, -30929> velocity=<-5,  3>
+position=<-51564,  31198> velocity=< 5, -3>
+position=<-30841,  10490> velocity=< 3, -1>
+position=<-51573,  20841> velocity=< 5, -2>
+position=< 10589, -41281> velocity=<-1,  4>
+position=< 20900,  31199> velocity=<-2, -3>
+position=<-41225,  20850> velocity=< 4, -2>
+position=< 20905, -30932> velocity=<-2,  3>
+position=< 41663, -30933> velocity=<-4,  3>
+position=<-30833, -51641> velocity=< 3,  5>
+position=<-41190, -30934> velocity=< 4,  3>
+position=< 31291, -51637> velocity=<-3,  5>
+position=<-10132, -10220> velocity=< 1,  1>
+position=< 41638, -10220> velocity=<-4,  1>
+position=<-10165,  10486> velocity=< 1, -1>
+position=<-20504,  41552> velocity=< 2, -4>
+position=< 20929,  41559> velocity=<-2, -4>
+position=<-10124,  31200> velocity=< 1, -3>
+position=< 31279,  41556> velocity=<-3, -4>
+position=< 20908, -51637> velocity=<-2,  5>
+position=<-30867,  31199> velocity=< 3, -3>
+position=< 41628, -51640> velocity=<-4,  5>
+position=< 52001,  41555> velocity=<-5, -4>
+position=<-10139, -10215> velocity=< 1,  1>
+position=< 51957, -10224> velocity=<-5,  1>
+position=< 52018,  51912> velocity=<-5, -5>
+position=< 20892, -20572> velocity=<-2,  2>
+position=<-20520,  31204> velocity=< 2, -3>
+position=<-51585,  10488> velocity=< 5, -1>
+position=< 31251,  51910> velocity=<-3, -5>
+position=<-41204,  10490> velocity=< 4, -1>
+position=<-51566,  20845> velocity=< 5, -2>
+position=< 41623,  41553> velocity=<-4, -4>
+position=<-41233, -51641> velocity=< 4,  5>
+position=<-41194,  31203> velocity=< 4, -3>
+position=<-10112, -20570> velocity=< 1,  2>
+position=< 10574,  10488> velocity=<-1, -1>
+position=< 41646, -20573> velocity=<-4,  2>
+position=<-51560,  31200> velocity=< 5, -3>
+position=<-41214, -41281> velocity=< 4,  4>
+position=< 31308, -51635> velocity=<-3,  5>
+position=<-10130,  10486> velocity=< 1, -1>
+position=<-51588, -41285> velocity=< 5,  4>
+position=< 41643,  10486> velocity=<-4, -1>
+position=<-51588,  31197> velocity=< 5, -3>
+position=< 31284, -10218> velocity=<-3,  1>
+position=< 31295,  31199> velocity=<-3, -3>
+position=<-20509,  41555> velocity=< 2, -4>
+position=<-51588,  41558> velocity=< 5, -4>
+position=<-20486,  31196> velocity=< 2, -3>
+position=< 10569, -41288> velocity=<-1,  4>
+position=< 10593,  51907> velocity=<-1, -5>
+position=< 20925,  31200> velocity=<-2, -3>
+position=<-20480, -30929> velocity=< 2,  3>
+position=<-30870,  20842> velocity=< 3, -2>
+position=<-30843,  31202> velocity=< 3, -3>
+position=< 41631, -41288> velocity=<-4,  4>
+position=<-30854,  41559> velocity=< 3, -4>
+position=<-30867, -20575> velocity=< 3,  2>
+position=<-30859,  20844> velocity=< 3, -2>
+position=< 41661, -51639> velocity=<-4,  5>
+position=< 31259, -20575> velocity=<-3,  2>
+position=< 31276,  51906> velocity=<-3, -5>
+position=< 41615,  10494> velocity=<-4, -1>
+position=<-41226, -51644> velocity=< 4,  5>
+position=< 20932,  31204> velocity=<-2, -3>
+position=< 31271,  51910> velocity=<-3, -5>
+position=< 10561,  51906> velocity=<-1, -5>
+position=< 41650, -10216> velocity=<-4,  1>
+position=<-51582, -20574> velocity=< 5,  2>
+position=< 41647,  20843> velocity=<-4, -2>
+position=<-41227,  31200> velocity=< 4, -3>
+position=< 20926, -41289> velocity=<-2,  4>
+position=< 51986,  51915> velocity=<-5, -5>
+position=< 41642, -51644> velocity=<-4,  5>
+position=< 20908, -30933> velocity=<-2,  3>
+position=< 10569,  51906> velocity=<-1, -5>
+position=<-10162,  41555> velocity=< 1, -4>
+position=< 31276,  51912> velocity=<-3, -5>
+position=< 41636,  20850> velocity=<-4, -2>
+position=< 51973, -51641> velocity=<-5,  5>
+position=<-51588, -10219> velocity=< 5,  1>
+position=< 20948, -41283> velocity=<-2,  4>
+position=<-30822, -30933> velocity=< 3,  3>
+position=<-30823,  51911> velocity=< 3, -5>
+position=<-30875,  51908> velocity=< 3, -5>
+position=< 31287, -10222> velocity=<-3,  1>
+position=< 31247, -10222> velocity=<-3,  1>
+position=<-10116, -41289> velocity=< 1,  4>
+position=< 41643,  10490> velocity=<-4, -1>
+position=< 20925,  31200> velocity=<-2, -3>
+position=<-30854,  20848> velocity=< 3, -2>
+position=<-30872,  20841> velocity=< 3, -2>
+position=< 41610,  41552> velocity=<-4, -4>
+position=<-41228,  20845> velocity=< 4, -2>
+position=< 31283,  10486> velocity=<-3, -1>
+position=<-20512,  41558> velocity=< 2, -4>
+position=< 41623,  41554> velocity=<-4, -4>
+position=< 51973, -30933> velocity=<-5,  3>
+position=< 31296,  20845> velocity=<-3, -2>
+position=< 41642, -10215> velocity=<-4,  1>
+position=<-20517, -41285> velocity=< 2,  4>
+position=<-51548,  31197> velocity=< 5, -3>
+position=< 20897, -51642> velocity=<-2,  5>
+position=< 41642,  51915> velocity=<-4, -5>
+position=< 20935,  41555> velocity=<-2, -4>
+position=< 20932,  41555> velocity=<-2, -4>
+position=< 41634, -20577> velocity=<-4,  2>
+position=<-51574,  10486> velocity=< 5, -1>
+position=< 52002,  10494> velocity=<-5, -1>
+position=< 51997,  41559> velocity=<-5, -4>
+position=< 41610, -10221> velocity=<-4,  1>
+position=< 20909, -41289> velocity=<-2,  4>
+position=< 31268,  10489> velocity=<-3, -1>
+position=<-51577,  10492> velocity=< 5, -1>
+position=<-51574, -20579> velocity=< 5,  2>
+position=< 31263, -51640> velocity=<-3,  5>
+position=<-51564,  10494> velocity=< 5, -1>
+position=< 41607,  20846> velocity=<-4, -2>
+position=<-10114, -51644> velocity=< 1,  5>
+position=< 31295,  20842> velocity=<-3, -2>
+position=< 31247,  51907> velocity=<-3, -5>
+position=<-10164,  41555> velocity=< 1, -4>
+position=<-30827, -51639> velocity=< 3,  5>
+position=<-41233,  10488> velocity=< 4, -1>
+position=< 52001,  41558> velocity=<-5, -4>
+position=< 51981,  41551> velocity=<-5, -4>
+position=<-20469,  10486> velocity=< 2, -1>
+position=< 41607,  10490> velocity=<-4, -1>
+position=< 20953, -20574> velocity=<-2,  2>
+position=< 41631,  41556> velocity=<-4, -4>
+position=< 20929, -41286> velocity=<-2,  4>
+position=<-51574, -41289> velocity=< 5,  4>
+position=<-30851,  20848> velocity=< 3, -2>
+position=< 51997, -41287> velocity=<-5,  4>
+position=<-41206,  31200> velocity=< 4, -3>
+position=<-10141, -51641> velocity=< 1,  5>
+position=< 31271,  10488> velocity=<-3, -1>
+position=<-20491, -30929> velocity=< 2,  3>
+position=< 10574,  51912> velocity=<-1, -5>
+position=< 51986, -30930> velocity=<-5,  3>
+position=<-51585,  20847> velocity=< 5, -2>
+position=< 31303, -10216> velocity=<-3,  1>
+position=< 10574, -20572> velocity=<-1,  2>
+position=< 41622,  31196> velocity=<-4, -3>
+position=< 52005, -51635> velocity=<-5,  5>
+position=< 52007,  20847> velocity=<-5, -2>
+position=< 31300,  10486> velocity=<-3, -1>
+position=<-41203,  51906> velocity=< 4, -5>
+position=< 31303, -30929> velocity=<-3,  3>
+position=<-30882,  10490> velocity=< 3, -1>
+position=<-10128,  41559> velocity=< 1, -4>
+position=< 31271, -10215> velocity=<-3,  1>
+position=<-20468,  20846> velocity=< 2, -2>
+position=<-20500, -20575> velocity=< 2,  2>
+position=< 31308, -20578> velocity=<-3,  2>
+position=<-10125, -20575> velocity=< 1,  2>
+position=< 10594,  31196> velocity=<-1, -3>
+position=< 41629,  31200> velocity=<-4, -3>
+position=<-10161, -10217> velocity=< 1,  1>
+position=< 31295,  41559> velocity=<-3, -4>
+position=<-51537, -30931> velocity=< 5,  3>
+position=<-51585,  10491> velocity=< 5, -1>
+position=<-41182,  31198> velocity=< 4, -3>
+position=<-41238,  20850> velocity=< 4, -2>
+position=<-51556,  10491> velocity=< 5, -1>
+position=< 51970,  31199> velocity=<-5, -3>
+position=< 41607,  31204> velocity=<-4, -3>
+position=< 10569,  31198> velocity=<-1, -3>
+position=<-20491,  31197> velocity=< 2, -3>
+position=<-30848,  31200> velocity=< 3, -3>
+position=<-20494, -10224> velocity=< 2,  1>
+position=<-10131,  20841> velocity=< 1, -2>
+position=<-10130, -51640> velocity=< 1,  5>
+position=<-10141,  31202> velocity=< 1, -3>
+position=<-41188, -51638> velocity=< 4,  5>
+position=< 20904, -30928> velocity=<-2,  3>
+position=< 20929, -51641> velocity=<-2,  5>
+position=<-10114, -20570> velocity=< 1,  2>
+position=<-30883, -10219> velocity=< 3,  1>
+position=< 41607,  10495> velocity=<-4, -1>
+position=<-30851,  51912> velocity=< 3, -5>
+position=< 20937,  20850> velocity=<-2, -2>
+position=<-51575, -30930> velocity=< 5,  3>
+position=< 10566, -30925> velocity=<-1,  3>
+position=< 51986, -20572> velocity=<-5,  2>
+position=<-51569, -51637> velocity=< 5,  5>
+position=<-10130, -20574> velocity=< 1,  2>
+position=<-10157, -41281> velocity=< 1,  4>
+position=<-41217,  31197> velocity=< 4, -3>
+position=<-30835,  41556> velocity=< 3, -4>
+position=<-10128, -20576> velocity=< 1,  2>
+position=< 20940,  10487> velocity=<-2, -1>
+position=< 52013,  41556> velocity=<-5, -4>
+position=<-51588, -20570> velocity=< 5,  2>
+position=<-10157, -51641> velocity=< 1,  5>
+position=<-10112,  20847> velocity=< 1, -2>
+position=<-51573, -10220> velocity=< 5,  1>
+position=< 20917,  51910> velocity=<-2, -5>
+position=< 51997, -30933> velocity=<-5,  3>
+position=< 41618,  20845> velocity=<-4, -2>
+position=<-41205, -10224> velocity=< 4,  1>
+position=<-30875,  41555> velocity=< 3, -4>
+position=<-51559,  20845> velocity=< 5, -2>
+position=< 31256,  51906> velocity=<-3, -5>
+position=< 10590,  41560> velocity=<-1, -4>
+position=< 41610, -41284> velocity=<-4,  4>
+position=<-30842,  31196> velocity=< 3, -3>
+position=< 10542, -10218> velocity=<-1,  1>
+position=<-41227, -41289> velocity=< 4,  4>
+position=< 10566,  10492> velocity=<-1, -1>
+position=< 10546, -20579> velocity=<-1,  2>
+position=< 10537, -30927> velocity=<-1,  3>
+position=<-51564,  51910> velocity=< 5, -5>
+position=<-41222, -41289> velocity=< 4,  4>
+position=< 10593, -10217> velocity=<-1,  1>
+position=<-20528, -20575> velocity=< 2,  2>