172 lines
3.7 KiB
C

/*
* Copyright (c) Camden Dixie O'Brien
* SPDX-License-Identifier: AGPL-3.0-only
*/
#include "puzz.h"
#include <assert.h>
#include <stdbool.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <sys/time.h>
#define NRUNS 10000
#define MAX_ADJ 8
typedef enum {
FOUND_MINE,
FOUND_SAFE,
FOUND_NOTHING,
} update_res_t;
typedef struct {
int mines, unknown;
puzz_t field;
} state_t;
static void setadj(puzz_t field, int x, int y, uint8_t from, uint8_t to)
{
for (int yp = y - 1; yp < y + 2; ++yp) {
for (int xp = x - 1; xp < x + 2; ++xp) {
if (xp < 0 || xp >= WIDTH || yp < 0 || yp >= HEIGHT
|| (xp == x && yp == y))
continue;
field[xp][yp] = field[xp][yp] == from ? to : field[xp][yp];
}
}
}
static void getadj(puzz_t field, int *x, int *y, uint8_t val)
{
for (int yp = *y - 1; yp < *y + 2; ++yp) {
for (int xp = *x - 1; xp < *x + 2; ++xp) {
if (xp < 0 || xp >= WIDTH || yp < 0 || yp >= HEIGHT
|| (xp == *x && yp == *y))
continue;
if (field[xp][yp] == val) {
*x = xp;
*y = yp;
return;
}
}
}
assert(false);
}
static update_res_t update(state_t *state, int *x_out, int *y_out)
{
state->mines = state->unknown = 0;
for (int y = 0; y < HEIGHT; ++y) {
for (int x = 0; x < WIDTH; ++x) {
state->mines += state->field[x][y] == MINE ? 1 : 0;
state->unknown += state->field[x][y] == UNKNOWN ? 1 : 0;
}
}
for (int y = 0; y < HEIGHT; ++y) {
for (int x = 0; x < WIDTH; ++x) {
if (state->field[x][y] == UNKNOWN || state->field[x][y] == MINE
|| state->field[x][y] == SAFE)
continue;
const int mines = countadj(state->field, x, y, MINE);
const int unknowns = countadj(state->field, x, y, UNKNOWN);
if (unknowns == 0)
continue;
if (mines + unknowns == state->field[x][y]) {
setadj(state->field, x, y, UNKNOWN, MINE);
return FOUND_MINE;
}
if (mines == state->field[x][y]) {
getadj(state->field, &x, &y, UNKNOWN);
*x_out = x;
*y_out = y;
return FOUND_SAFE;
}
}
}
return FOUND_NOTHING;
}
static status_t solve(int *turns_out)
{
state_t state = { .mines = 0, .unknown = WIDTH * HEIGHT };
memset(state.field, UNKNOWN, sizeof(puzz_t));
int x = rand() % WIDTH;
int y = rand() % HEIGHT;
status_t status;
int turns = 0;
do {
++turns;
if (state.field[x][y] != MINE
&& (status = probe(x, y, state.field)) == DEAD)
break;
update_res_t update_res;
do
update_res = update(&state, &x, &y);
while (update_res == FOUND_MINE);
if (check(state.field) != OK) {
printf("Incorrect inference! State:\n");
print(state.field);
printsoln();
return INCORRECT;
}
if (update_res == FOUND_NOTHING) {
x = rand() % WIDTH;
y = rand() % HEIGHT;
}
} while (state.mines < NMINES || state.unknown > 0);
*turns_out = turns;
return status;
}
int main(void)
{
struct timeval tv;
if (gettimeofday(&tv, NULL) != 0) {
perror("Failed to get time");
exit(1);
}
srand(tv.tv_usec);
int nsolved = 0, nfirst = 0, nincorrect = 0, turns;
for (int i = 0; i < NRUNS; ++i) {
gen();
switch (solve(&turns)) {
case DEAD:
if (turns == 1)
++nfirst;
break;
case OK:
++nsolved;
break;
case INCORRECT:
++nincorrect;
break;
}
}
const double solved_prop = (double)nsolved / NRUNS;
const double incorrect_prop = (double)nincorrect / NRUNS;
const double first_prop = (double)nfirst / NRUNS;
const double solved_safe_start_prop = (double)nsolved / (NRUNS - nfirst);
printf("Solved %d (%0.1f%%)\n", nsolved, 100 * solved_prop);
printf("%d incorrect (%0.1f%%)\n", nincorrect, 100 * incorrect_prop);
printf("%d died on first turn (%0.1f%%)\n", nfirst, 100 * first_prop);
printf("%0.1f%% solved on safe start\n", 100 * solved_safe_start_prop);
return 0;
}