135 lines
3.4 KiB
Nim
135 lines
3.4 KiB
Nim
import math, options, sequtils, random, sets
|
|
import combinators, game, fixedseq
|
|
|
|
|
|
type
|
|
ScoreSet* = array[Color, int]
|
|
|
|
ScoreSpread = object
|
|
lo: array[Color, float]
|
|
hi: array[Color, float]
|
|
|
|
LegResults* = tuple[scores: ScoreSet, endStates: HashSet[Board]]
|
|
|
|
|
|
proc update*(scores: var ScoreSet, toAdd: ScoreSet) =
|
|
for i, s in toAdd:
|
|
scores[i] += s
|
|
|
|
|
|
proc projectLeg(b: Board): LegResults =
|
|
var scores: ScoreSet
|
|
var endStates: HashSet[Board]
|
|
|
|
var diceRemaining: ColorStack
|
|
diceRemaining.initFixedSeq
|
|
for i, c in b.diceRolled:
|
|
if not c: diceRemaining.add(i)
|
|
|
|
for future in possibleFutures(diceRemaining):
|
|
var prediction = b # make a copy
|
|
for dieRoll in future:
|
|
prediction.advance(dieRoll)
|
|
inc scores[prediction.leader.get]
|
|
# deduplicate results
|
|
endStates.incl(prediction)
|
|
|
|
result = (scores, endStates)
|
|
|
|
|
|
proc projectOutcomes(b: Board, maxDepth = 1): ScoreSet =
|
|
var outcomeStack = @[ [b].toHashSet ]
|
|
for depth in 1..maxDepth:
|
|
echo "simulating ", outcomeStack[^1].len, " possible legs."
|
|
var endStates: HashSet[Board]
|
|
|
|
for o in outcomeStack[^1]:
|
|
var o = o # make it mutable
|
|
if outcomeStack.len > 1:
|
|
o.resetDice # o was describina an end-of-leg state, so dice were exhausted
|
|
|
|
let projection = o.projectLeg
|
|
result.update(projection[0])
|
|
endStates.incl(projection[1])
|
|
stdout.write("simulated: " & $result.sum & "\r")
|
|
|
|
outcomeStack.add(endStates)
|
|
echo "\nDistinct end states: ", outcomeStack.mapIt(it.len).sum
|
|
|
|
|
|
proc randomGame(b: Board, r: var Rand): Color =
|
|
var projection = b
|
|
while true:
|
|
for roll in randomFuture(projection.diceRemaining, r):
|
|
projection.advance(roll)
|
|
if projection.gameOver:
|
|
return projection.leader.get
|
|
projection.resetDice
|
|
|
|
|
|
proc randomGames(b: Board, count: SomeInteger): ScoreSet =
|
|
randomize()
|
|
var r = initRand(rand(int64))
|
|
for i in 1 .. count:
|
|
let winner = b.randomGame(r)
|
|
inc result[winner]
|
|
# if i mod 100_000 == 0 or i == count - 1:
|
|
# stdout.write("simulating " & count & "random games: " & $i & "\r")
|
|
# echo ""
|
|
|
|
|
|
proc randomSpread(b: Board, nTests: SomeInteger, nSamples: SomeInteger): ScoreSpread =
|
|
for s in result.lo.mitems:
|
|
s = 1
|
|
|
|
for i in 0 ..< nTests:
|
|
let scores = b.randomGames(nSamples)
|
|
let total = scores.sum
|
|
for color, score in scores:
|
|
let pct = score / total
|
|
if pct < result.lo[color]:
|
|
result.lo[color] = pct
|
|
if pct > result.hi[color]:
|
|
result.hi[color] = pct
|
|
|
|
|
|
var b: Board
|
|
b.init
|
|
|
|
randomize()
|
|
var r = initRand(rand(int64))
|
|
|
|
var rolls: array[5, tuple[c: Color, p: int]]
|
|
for i, roll in randomFuture(b.diceRemaining, r):
|
|
rolls[i] = (roll[0], roll[1] + 1)
|
|
|
|
b.setState(rolls, @[])
|
|
echo "Starting state:"
|
|
b.display(1, 5)
|
|
|
|
# block outer:
|
|
# while true:
|
|
# for roll in randomFuture(b.diceRemaining, r):
|
|
# b.advance(roll)
|
|
# if b.gameOver:
|
|
# echo "last roll: ", roll
|
|
# break outer
|
|
# b.resetDice
|
|
|
|
# echo "winner: ", b.leader.get
|
|
# b.display(min(b.camels) - 1, 16)
|
|
|
|
# let scores = b.randomGames(10_000_000)
|
|
# let total = scores.sum
|
|
# for i, c in scores:
|
|
# echo Color(i), ": ", (100 * c / total).round(2), "% (", c, " / ", total, ")"
|
|
|
|
for i in 0..3:
|
|
let samples = 10 ^ (i + 3)
|
|
echo "Simulating ", samples, " games 100 times"
|
|
let spread = b.randomSpread(100, samples)
|
|
for color, lo in spread.lo:
|
|
stdout.write($color & ": " & $lo.round(4) & "-" & $spread.hi[color].round(4) & ", ")
|
|
stdout.flushFile
|
|
echo "\n"
|