171 lines
3.6 KiB
C
171 lines
3.6 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)
|
|
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;
|
|
}
|