@stevenvo780/st-lang 4.5.0 → 4.5.2
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/README.md +78 -37
- package/dist/bdd/index.d.ts +4 -0
- package/dist/bdd/index.d.ts.map +1 -0
- package/dist/bdd/index.js +29 -0
- package/dist/bdd/index.js.map +1 -0
- package/dist/bdd/manager.d.ts +122 -0
- package/dist/bdd/manager.d.ts.map +1 -0
- package/dist/bdd/manager.js +524 -0
- package/dist/bdd/manager.js.map +1 -0
- package/dist/bdd/types.d.ts +20 -0
- package/dist/bdd/types.d.ts.map +1 -0
- package/dist/bdd/types.js +28 -0
- package/dist/bdd/types.js.map +1 -0
- package/dist/game-theory/common-games.d.ts +50 -0
- package/dist/game-theory/common-games.d.ts.map +1 -0
- package/dist/game-theory/common-games.js +112 -0
- package/dist/game-theory/common-games.js.map +1 -0
- package/dist/game-theory/dominance.d.ts +26 -0
- package/dist/game-theory/dominance.d.ts.map +1 -0
- package/dist/game-theory/dominance.js +203 -0
- package/dist/game-theory/dominance.js.map +1 -0
- package/dist/game-theory/index.d.ts +8 -0
- package/dist/game-theory/index.d.ts.map +1 -0
- package/dist/game-theory/index.js +47 -0
- package/dist/game-theory/index.js.map +1 -0
- package/dist/game-theory/lemke-howson.d.ts +3 -0
- package/dist/game-theory/lemke-howson.d.ts.map +1 -0
- package/dist/game-theory/lemke-howson.js +269 -0
- package/dist/game-theory/lemke-howson.js.map +1 -0
- package/dist/game-theory/linalg.d.ts +7 -0
- package/dist/game-theory/linalg.d.ts.map +1 -0
- package/dist/game-theory/linalg.js +69 -0
- package/dist/game-theory/linalg.js.map +1 -0
- package/dist/game-theory/pure-nash.d.ts +3 -0
- package/dist/game-theory/pure-nash.d.ts.map +1 -0
- package/dist/game-theory/pure-nash.js +86 -0
- package/dist/game-theory/pure-nash.js.map +1 -0
- package/dist/game-theory/support-enumeration.d.ts +3 -0
- package/dist/game-theory/support-enumeration.d.ts.map +1 -0
- package/dist/game-theory/support-enumeration.js +214 -0
- package/dist/game-theory/support-enumeration.js.map +1 -0
- package/dist/game-theory/types.d.ts +57 -0
- package/dist/game-theory/types.d.ts.map +1 -0
- package/dist/game-theory/types.js +191 -0
- package/dist/game-theory/types.js.map +1 -0
- package/dist/index.d.ts +19 -1
- package/dist/index.d.ts.map +1 -1
- package/dist/index.js +75 -1
- package/dist/index.js.map +1 -1
- package/dist/information-theory/index.d.ts +29 -0
- package/dist/information-theory/index.d.ts.map +1 -0
- package/dist/information-theory/index.js +323 -0
- package/dist/information-theory/index.js.map +1 -0
- package/dist/namespaces/logic.d.ts +42 -0
- package/dist/namespaces/logic.d.ts.map +1 -0
- package/dist/namespaces/logic.js +112 -0
- package/dist/namespaces/logic.js.map +1 -0
- package/dist/namespaces/proof-systems.d.ts +18 -0
- package/dist/namespaces/proof-systems.d.ts.map +1 -0
- package/dist/namespaces/proof-systems.js +59 -0
- package/dist/namespaces/proof-systems.js.map +1 -0
- package/dist/namespaces/reasoning.d.ts +26 -0
- package/dist/namespaces/reasoning.d.ts.map +1 -0
- package/dist/namespaces/reasoning.js +72 -0
- package/dist/namespaces/reasoning.js.map +1 -0
- package/dist/namespaces/semantics.d.ts +21 -0
- package/dist/namespaces/semantics.d.ts.map +1 -0
- package/dist/namespaces/semantics.js +64 -0
- package/dist/namespaces/semantics.js.map +1 -0
- package/dist/namespaces/solvers.d.ts +22 -0
- package/dist/namespaces/solvers.d.ts.map +1 -0
- package/dist/namespaces/solvers.js +66 -0
- package/dist/namespaces/solvers.js.map +1 -0
- package/dist/namespaces/type-theory.d.ts +22 -0
- package/dist/namespaces/type-theory.d.ts.map +1 -0
- package/dist/namespaces/type-theory.js +65 -0
- package/dist/namespaces/type-theory.js.map +1 -0
- package/dist/optimization/branch-and-bound.d.ts +11 -0
- package/dist/optimization/branch-and-bound.d.ts.map +1 -0
- package/dist/optimization/branch-and-bound.js +311 -0
- package/dist/optimization/branch-and-bound.js.map +1 -0
- package/dist/optimization/index.d.ts +5 -0
- package/dist/optimization/index.d.ts.map +1 -0
- package/dist/optimization/index.js +20 -0
- package/dist/optimization/index.js.map +1 -0
- package/dist/optimization/simplex.d.ts +8 -0
- package/dist/optimization/simplex.d.ts.map +1 -0
- package/dist/optimization/simplex.js +408 -0
- package/dist/optimization/simplex.js.map +1 -0
- package/dist/optimization/standard-form.d.ts +11 -0
- package/dist/optimization/standard-form.d.ts.map +1 -0
- package/dist/optimization/standard-form.js +112 -0
- package/dist/optimization/standard-form.js.map +1 -0
- package/dist/optimization/types.d.ts +44 -0
- package/dist/optimization/types.d.ts.map +1 -0
- package/dist/optimization/types.js +15 -0
- package/dist/optimization/types.js.map +1 -0
- package/dist/profiles/hol/connectives.d.ts +46 -0
- package/dist/profiles/hol/connectives.d.ts.map +1 -0
- package/dist/profiles/hol/connectives.js +104 -0
- package/dist/profiles/hol/connectives.js.map +1 -0
- package/dist/profiles/hol/index.d.ts +7 -0
- package/dist/profiles/hol/index.d.ts.map +1 -0
- package/dist/profiles/hol/index.js +77 -0
- package/dist/profiles/hol/index.js.map +1 -0
- package/dist/profiles/hol/rules.d.ts +31 -0
- package/dist/profiles/hol/rules.d.ts.map +1 -0
- package/dist/profiles/hol/rules.js +258 -0
- package/dist/profiles/hol/rules.js.map +1 -0
- package/dist/profiles/hol/term.d.ts +52 -0
- package/dist/profiles/hol/term.d.ts.map +1 -0
- package/dist/profiles/hol/term.js +345 -0
- package/dist/profiles/hol/term.js.map +1 -0
- package/dist/profiles/hol/type-system.d.ts +38 -0
- package/dist/profiles/hol/type-system.d.ts.map +1 -0
- package/dist/profiles/hol/type-system.js +133 -0
- package/dist/profiles/hol/type-system.js.map +1 -0
- package/dist/profiles/hol/types.d.ts +59 -0
- package/dist/profiles/hol/types.d.ts.map +1 -0
- package/dist/profiles/hol/types.js +26 -0
- package/dist/profiles/hol/types.js.map +1 -0
- package/dist/smt-lib/ast.d.ts +138 -0
- package/dist/smt-lib/ast.d.ts.map +1 -0
- package/dist/smt-lib/ast.js +73 -0
- package/dist/smt-lib/ast.js.map +1 -0
- package/dist/smt-lib/emitter.d.ts +12 -0
- package/dist/smt-lib/emitter.d.ts.map +1 -0
- package/dist/smt-lib/emitter.js +174 -0
- package/dist/smt-lib/emitter.js.map +1 -0
- package/dist/smt-lib/index.d.ts +6 -0
- package/dist/smt-lib/index.d.ts.map +1 -0
- package/dist/smt-lib/index.js +31 -0
- package/dist/smt-lib/index.js.map +1 -0
- package/dist/smt-lib/parser.d.ts +13 -0
- package/dist/smt-lib/parser.d.ts.map +1 -0
- package/dist/smt-lib/parser.js +614 -0
- package/dist/smt-lib/parser.js.map +1 -0
- package/dist/smt-lib/tokenizer.d.ts +16 -0
- package/dist/smt-lib/tokenizer.d.ts.map +1 -0
- package/dist/smt-lib/tokenizer.js +234 -0
- package/dist/smt-lib/tokenizer.js.map +1 -0
- package/dist/tests/api/namespaces.test.d.ts +9 -0
- package/dist/tests/api/namespaces.test.d.ts.map +1 -0
- package/dist/tests/api/namespaces.test.js +218 -0
- package/dist/tests/api/namespaces.test.js.map +1 -0
- package/dist/tests/bdd/manager.test.d.ts +2 -0
- package/dist/tests/bdd/manager.test.d.ts.map +1 -0
- package/dist/tests/bdd/manager.test.js +472 -0
- package/dist/tests/bdd/manager.test.js.map +1 -0
- package/dist/tests/game-theory/game-theory.test.d.ts +2 -0
- package/dist/tests/game-theory/game-theory.test.d.ts.map +1 -0
- package/dist/tests/game-theory/game-theory.test.js +252 -0
- package/dist/tests/game-theory/game-theory.test.js.map +1 -0
- package/dist/tests/hol/hol.test.d.ts +2 -0
- package/dist/tests/hol/hol.test.d.ts.map +1 -0
- package/dist/tests/hol/hol.test.js +340 -0
- package/dist/tests/hol/hol.test.js.map +1 -0
- package/dist/tests/information-theory/information-theory.test.d.ts +2 -0
- package/dist/tests/information-theory/information-theory.test.d.ts.map +1 -0
- package/dist/tests/information-theory/information-theory.test.js +291 -0
- package/dist/tests/information-theory/information-theory.test.js.map +1 -0
- package/dist/tests/optimization/ilp.test.d.ts +2 -0
- package/dist/tests/optimization/ilp.test.d.ts.map +1 -0
- package/dist/tests/optimization/ilp.test.js +204 -0
- package/dist/tests/optimization/ilp.test.js.map +1 -0
- package/dist/tests/optimization/lp.test.d.ts +2 -0
- package/dist/tests/optimization/lp.test.d.ts.map +1 -0
- package/dist/tests/optimization/lp.test.js +224 -0
- package/dist/tests/optimization/lp.test.js.map +1 -0
- package/dist/tests/properties/agm.property.test.d.ts +2 -0
- package/dist/tests/properties/agm.property.test.d.ts.map +1 -0
- package/dist/tests/properties/agm.property.test.js +75 -0
- package/dist/tests/properties/agm.property.test.js.map +1 -0
- package/dist/tests/properties/anti-unification.property.test.d.ts +2 -0
- package/dist/tests/properties/anti-unification.property.test.d.ts.map +1 -0
- package/dist/tests/properties/anti-unification.property.test.js +43 -0
- package/dist/tests/properties/anti-unification.property.test.js.map +1 -0
- package/dist/tests/properties/argumentation.property.test.d.ts +2 -0
- package/dist/tests/properties/argumentation.property.test.d.ts.map +1 -0
- package/dist/tests/properties/argumentation.property.test.js +51 -0
- package/dist/tests/properties/argumentation.property.test.js.map +1 -0
- package/dist/tests/properties/bayesian.property.test.d.ts +2 -0
- package/dist/tests/properties/bayesian.property.test.d.ts.map +1 -0
- package/dist/tests/properties/bayesian.property.test.js +32 -0
- package/dist/tests/properties/bayesian.property.test.js.map +1 -0
- package/dist/tests/properties/bisimulation.property.test.d.ts +2 -0
- package/dist/tests/properties/bisimulation.property.test.d.ts.map +1 -0
- package/dist/tests/properties/bisimulation.property.test.js +50 -0
- package/dist/tests/properties/bisimulation.property.test.js.map +1 -0
- package/dist/tests/properties/cdcl.property.test.d.ts +2 -0
- package/dist/tests/properties/cdcl.property.test.d.ts.map +1 -0
- package/dist/tests/properties/cdcl.property.test.js +70 -0
- package/dist/tests/properties/cdcl.property.test.js.map +1 -0
- package/dist/tests/properties/coinduction.property.test.d.ts +2 -0
- package/dist/tests/properties/coinduction.property.test.d.ts.map +1 -0
- package/dist/tests/properties/coinduction.property.test.js +34 -0
- package/dist/tests/properties/coinduction.property.test.js.map +1 -0
- package/dist/tests/properties/constructive-reals.property.test.d.ts +2 -0
- package/dist/tests/properties/constructive-reals.property.test.d.ts.map +1 -0
- package/dist/tests/properties/constructive-reals.property.test.js +59 -0
- package/dist/tests/properties/constructive-reals.property.test.js.map +1 -0
- package/dist/tests/properties/csp.property.test.d.ts +2 -0
- package/dist/tests/properties/csp.property.test.d.ts.map +1 -0
- package/dist/tests/properties/csp.property.test.js +58 -0
- package/dist/tests/properties/csp.property.test.js.map +1 -0
- package/dist/tests/properties/generators.d.ts +78 -0
- package/dist/tests/properties/generators.d.ts.map +1 -0
- package/dist/tests/properties/generators.js +348 -0
- package/dist/tests/properties/generators.js.map +1 -0
- package/dist/tests/properties/ho-unify.property.test.d.ts +2 -0
- package/dist/tests/properties/ho-unify.property.test.d.ts.map +1 -0
- package/dist/tests/properties/ho-unify.property.test.js +46 -0
- package/dist/tests/properties/ho-unify.property.test.js.map +1 -0
- package/dist/tests/properties/hyperreal.property.test.d.ts +2 -0
- package/dist/tests/properties/hyperreal.property.test.d.ts.map +1 -0
- package/dist/tests/properties/hyperreal.property.test.js +33 -0
- package/dist/tests/properties/hyperreal.property.test.js.map +1 -0
- package/dist/tests/properties/intuit-nj.property.test.d.ts +2 -0
- package/dist/tests/properties/intuit-nj.property.test.d.ts.map +1 -0
- package/dist/tests/properties/intuit-nj.property.test.js +57 -0
- package/dist/tests/properties/intuit-nj.property.test.js.map +1 -0
- package/dist/tests/properties/lambda-calc.property.test.d.ts +2 -0
- package/dist/tests/properties/lambda-calc.property.test.d.ts.map +1 -0
- package/dist/tests/properties/lambda-calc.property.test.js +35 -0
- package/dist/tests/properties/lambda-calc.property.test.js.map +1 -0
- package/dist/tests/properties/mln.property.test.d.ts +2 -0
- package/dist/tests/properties/mln.property.test.d.ts.map +1 -0
- package/dist/tests/properties/mln.property.test.js +41 -0
- package/dist/tests/properties/mln.property.test.js.map +1 -0
- package/dist/tests/properties/mltt.property.test.d.ts +2 -0
- package/dist/tests/properties/mltt.property.test.d.ts.map +1 -0
- package/dist/tests/properties/mltt.property.test.js +33 -0
- package/dist/tests/properties/mltt.property.test.js.map +1 -0
- package/dist/tests/properties/nbe.property.test.d.ts +2 -0
- package/dist/tests/properties/nbe.property.test.d.ts.map +1 -0
- package/dist/tests/properties/nbe.property.test.js +44 -0
- package/dist/tests/properties/nbe.property.test.js.map +1 -0
- package/dist/tests/properties/planning.property.test.d.ts +2 -0
- package/dist/tests/properties/planning.property.test.d.ts.map +1 -0
- package/dist/tests/properties/planning.property.test.js +51 -0
- package/dist/tests/properties/planning.property.test.js.map +1 -0
- package/dist/tests/properties/profile-bridge.property.test.d.ts +2 -0
- package/dist/tests/properties/profile-bridge.property.test.d.ts.map +1 -0
- package/dist/tests/properties/profile-bridge.property.test.js +71 -0
- package/dist/tests/properties/profile-bridge.property.test.js.map +1 -0
- package/dist/tests/properties/refinement-types.property.test.d.ts +2 -0
- package/dist/tests/properties/refinement-types.property.test.d.ts.map +1 -0
- package/dist/tests/properties/refinement-types.property.test.js +37 -0
- package/dist/tests/properties/refinement-types.property.test.js.map +1 -0
- package/dist/tests/properties/sequent-g3.property.test.d.ts +2 -0
- package/dist/tests/properties/sequent-g3.property.test.d.ts.map +1 -0
- package/dist/tests/properties/sequent-g3.property.test.js +34 -0
- package/dist/tests/properties/sequent-g3.property.test.js.map +1 -0
- package/dist/tests/properties/symbolic-diff.property.test.d.ts +2 -0
- package/dist/tests/properties/symbolic-diff.property.test.d.ts.map +1 -0
- package/dist/tests/properties/symbolic-diff.property.test.js +49 -0
- package/dist/tests/properties/symbolic-diff.property.test.js.map +1 -0
- package/dist/tests/properties/system-f.property.test.d.ts +2 -0
- package/dist/tests/properties/system-f.property.test.d.ts.map +1 -0
- package/dist/tests/properties/system-f.property.test.js +37 -0
- package/dist/tests/properties/system-f.property.test.js.map +1 -0
- package/dist/tests/properties/theorem-cache.property.test.d.ts +2 -0
- package/dist/tests/properties/theorem-cache.property.test.d.ts.map +1 -0
- package/dist/tests/properties/theorem-cache.property.test.js +38 -0
- package/dist/tests/properties/theorem-cache.property.test.js.map +1 -0
- package/dist/tests/properties/trs.property.test.d.ts +2 -0
- package/dist/tests/properties/trs.property.test.d.ts.map +1 -0
- package/dist/tests/properties/trs.property.test.js +34 -0
- package/dist/tests/properties/trs.property.test.js.map +1 -0
- package/dist/tests/properties/unification.property.test.d.ts +2 -0
- package/dist/tests/properties/unification.property.test.d.ts.map +1 -0
- package/dist/tests/properties/unification.property.test.js +42 -0
- package/dist/tests/properties/unification.property.test.js.map +1 -0
- package/dist/tests/smt-lib/emitter.test.d.ts +2 -0
- package/dist/tests/smt-lib/emitter.test.d.ts.map +1 -0
- package/dist/tests/smt-lib/emitter.test.js +155 -0
- package/dist/tests/smt-lib/emitter.test.js.map +1 -0
- package/dist/tests/smt-lib/parser.test.d.ts +2 -0
- package/dist/tests/smt-lib/parser.test.d.ts.map +1 -0
- package/dist/tests/smt-lib/parser.test.js +164 -0
- package/dist/tests/smt-lib/parser.test.js.map +1 -0
- package/dist/tests/smt-lib/tokenizer.test.d.ts +2 -0
- package/dist/tests/smt-lib/tokenizer.test.d.ts.map +1 -0
- package/dist/tests/smt-lib/tokenizer.test.js +59 -0
- package/dist/tests/smt-lib/tokenizer.test.js.map +1 -0
- package/package.json +3 -1
|
@@ -0,0 +1,269 @@
|
|
|
1
|
+
"use strict";
|
|
2
|
+
// ============================================================
|
|
3
|
+
// Lemke-Howson — algoritmo de pivot para Nash de 2-player games
|
|
4
|
+
// ============================================================
|
|
5
|
+
//
|
|
6
|
+
// Versión clásica con tableros. Idea:
|
|
7
|
+
//
|
|
8
|
+
// - Normalizamos pagos para que sean estrictamente positivos
|
|
9
|
+
// (sumando un constante grande). Esto no cambia los Nash.
|
|
10
|
+
// - Trabajamos con dos sistemas en variables x ∈ R^n, y ∈ R^m
|
|
11
|
+
// más slacks r ∈ R^n, s ∈ R^m:
|
|
12
|
+
// B^T x + r = 1 con x ≥ 0, r ≥ 0
|
|
13
|
+
// A y + s = 1 con y ≥ 0, s ≥ 0
|
|
14
|
+
// donde A es la matriz de pagos del jugador 1 (NxM) y B la
|
|
15
|
+
// de pagos del jugador 2 (NxM).
|
|
16
|
+
// - Cada variable tiene una "etiqueta" en {1..n+m}. En un
|
|
17
|
+
// vértice completamente etiquetado (Nash), para cada label l
|
|
18
|
+
// o bien la variable es no-básica (= 0) o su complemento lo es.
|
|
19
|
+
// - Empezar en el vértice trivial (x=0, y=0; r=s=1) — todas las
|
|
20
|
+
// etiquetas están "cubiertas" por las variables no-básicas.
|
|
21
|
+
// Soltamos una etiqueta `startLabel` y pivotamos hasta volver
|
|
22
|
+
// a un vértice completamente etiquetado.
|
|
23
|
+
//
|
|
24
|
+
// Este algoritmo encuentra UN equilibrio (no necesariamente todos).
|
|
25
|
+
// Cambiar `startLabel` puede llevar a equilibria distintos.
|
|
26
|
+
//
|
|
27
|
+
// Implementación pedagógica: O(n+m) pivotes en práctica.
|
|
28
|
+
// El soporte que ya teníamos (support enumeration) sirve para juegos
|
|
29
|
+
// chicos. Lemke-Howson es alternativa para tamaños medios o cuando
|
|
30
|
+
// uno quiere "algún" Nash rápido.
|
|
31
|
+
Object.defineProperty(exports, "__esModule", { value: true });
|
|
32
|
+
exports.lemkeHowson = lemkeHowson;
|
|
33
|
+
const types_1 = require("./types");
|
|
34
|
+
const TOL = 1e-9;
|
|
35
|
+
function lemkeHowson(game, startLabel = 0) {
|
|
36
|
+
if (game.players !== 2) {
|
|
37
|
+
throw new Error('lemkeHowson supports only 2-player games');
|
|
38
|
+
}
|
|
39
|
+
const n = game.strategies[0];
|
|
40
|
+
const m = game.strategies[1];
|
|
41
|
+
if (startLabel < 0 || startLabel >= n + m) {
|
|
42
|
+
throw new Error(`startLabel ${startLabel} out of range [0, ${n + m})`);
|
|
43
|
+
}
|
|
44
|
+
const a1 = game.payoffs[0];
|
|
45
|
+
const a2 = game.payoffs[1];
|
|
46
|
+
// Construir matrices A (jugador 1) y B (jugador 2)
|
|
47
|
+
const A = [];
|
|
48
|
+
const B = [];
|
|
49
|
+
for (let i = 0; i < n; i++) {
|
|
50
|
+
const rowA = [];
|
|
51
|
+
const rowB = [];
|
|
52
|
+
for (let j = 0; j < m; j++) {
|
|
53
|
+
rowA.push(a1[i * m + j]);
|
|
54
|
+
rowB.push(a2[i * m + j]);
|
|
55
|
+
}
|
|
56
|
+
A.push(rowA);
|
|
57
|
+
B.push(rowB);
|
|
58
|
+
}
|
|
59
|
+
// Shift para que sean positivos
|
|
60
|
+
let minVal = Infinity;
|
|
61
|
+
for (let i = 0; i < n; i++) {
|
|
62
|
+
for (let j = 0; j < m; j++) {
|
|
63
|
+
const av = A[i][j];
|
|
64
|
+
const bv = B[i][j];
|
|
65
|
+
if (av < minVal)
|
|
66
|
+
minVal = av;
|
|
67
|
+
if (bv < minVal)
|
|
68
|
+
minVal = bv;
|
|
69
|
+
}
|
|
70
|
+
}
|
|
71
|
+
const shift = minVal <= 0 ? 1 - minVal : 0;
|
|
72
|
+
if (shift !== 0) {
|
|
73
|
+
for (let i = 0; i < n; i++) {
|
|
74
|
+
for (let j = 0; j < m; j++) {
|
|
75
|
+
A[i][j] = A[i][j] + shift;
|
|
76
|
+
B[i][j] = B[i][j] + shift;
|
|
77
|
+
}
|
|
78
|
+
}
|
|
79
|
+
}
|
|
80
|
+
const labelOf = (v) => {
|
|
81
|
+
if (v.kind === 0 || v.kind === 1)
|
|
82
|
+
return v.index;
|
|
83
|
+
return n + v.index;
|
|
84
|
+
};
|
|
85
|
+
const complementOf = (v) => {
|
|
86
|
+
if (v.kind === 0)
|
|
87
|
+
return { kind: 1, index: v.index };
|
|
88
|
+
if (v.kind === 1)
|
|
89
|
+
return { kind: 0, index: v.index };
|
|
90
|
+
if (v.kind === 2)
|
|
91
|
+
return { kind: 3, index: v.index };
|
|
92
|
+
return { kind: 2, index: v.index };
|
|
93
|
+
};
|
|
94
|
+
// System 1: B^T para que las filas sean ecuaciones indexadas por j.
|
|
95
|
+
const S1 = { mat: [], rhs: [] };
|
|
96
|
+
for (let j = 0; j < m; j++) {
|
|
97
|
+
const row = [];
|
|
98
|
+
for (let i = 0; i < n; i++)
|
|
99
|
+
row.push(B[i][j]);
|
|
100
|
+
S1.mat.push(row);
|
|
101
|
+
S1.rhs.push(1);
|
|
102
|
+
}
|
|
103
|
+
const S1_basic = Array.from({ length: m }, (_, j) => ({ kind: 3, index: j }));
|
|
104
|
+
const S1_nonbasic = Array.from({ length: n }, (_, i) => ({ kind: 0, index: i }));
|
|
105
|
+
const S2 = { mat: [], rhs: [] };
|
|
106
|
+
for (let i = 0; i < n; i++) {
|
|
107
|
+
const row = [];
|
|
108
|
+
for (let j = 0; j < m; j++)
|
|
109
|
+
row.push(A[i][j]);
|
|
110
|
+
S2.mat.push(row);
|
|
111
|
+
S2.rhs.push(1);
|
|
112
|
+
}
|
|
113
|
+
const S2_basic = Array.from({ length: n }, (_, i) => ({ kind: 1, index: i }));
|
|
114
|
+
const S2_nonbasic = Array.from({ length: m }, (_, j) => ({ kind: 2, index: j }));
|
|
115
|
+
// Pivot: pasar una no-básica `enter` a básica en el sistema que le
|
|
116
|
+
// corresponde, y sacar de la base la que pierda con la regla del
|
|
117
|
+
// cociente mínimo. Devuelve la variable que sale (Var) o null si
|
|
118
|
+
// el sistema es ilimitado (no debería pasar con shift > 0).
|
|
119
|
+
function pivotIn(sys, enter) {
|
|
120
|
+
const S = sys === 'S1' ? S1 : S2;
|
|
121
|
+
const nonbasic = sys === 'S1' ? S1_nonbasic : S2_nonbasic;
|
|
122
|
+
const basic = sys === 'S1' ? S1_basic : S2_basic;
|
|
123
|
+
// columna de `enter` en `nonbasic`
|
|
124
|
+
let col = -1;
|
|
125
|
+
for (let c = 0; c < nonbasic.length; c++) {
|
|
126
|
+
const v = nonbasic[c];
|
|
127
|
+
if (v.kind === enter.kind && v.index === enter.index) {
|
|
128
|
+
col = c;
|
|
129
|
+
break;
|
|
130
|
+
}
|
|
131
|
+
}
|
|
132
|
+
if (col < 0)
|
|
133
|
+
return null;
|
|
134
|
+
// Min ratio: filas con coef > 0 de la columna entrante
|
|
135
|
+
let bestRow = -1;
|
|
136
|
+
let bestRatio = Infinity;
|
|
137
|
+
for (let r = 0; r < S.mat.length; r++) {
|
|
138
|
+
const coef = S.mat[r][col];
|
|
139
|
+
if (coef <= TOL)
|
|
140
|
+
continue;
|
|
141
|
+
const ratio = S.rhs[r] / coef;
|
|
142
|
+
if (ratio < bestRatio - TOL) {
|
|
143
|
+
bestRatio = ratio;
|
|
144
|
+
bestRow = r;
|
|
145
|
+
}
|
|
146
|
+
}
|
|
147
|
+
if (bestRow < 0)
|
|
148
|
+
return null; // unbounded
|
|
149
|
+
// Pivot Gauss-Jordan en (bestRow, col)
|
|
150
|
+
const pivot = S.mat[bestRow][col];
|
|
151
|
+
const pivRow = S.mat[bestRow];
|
|
152
|
+
for (let c = 0; c < pivRow.length; c++)
|
|
153
|
+
pivRow[c] = pivRow[c] / pivot;
|
|
154
|
+
S.rhs[bestRow] = S.rhs[bestRow] / pivot;
|
|
155
|
+
for (let r = 0; r < S.mat.length; r++) {
|
|
156
|
+
if (r === bestRow)
|
|
157
|
+
continue;
|
|
158
|
+
const rowR = S.mat[r];
|
|
159
|
+
const factor = rowR[col];
|
|
160
|
+
if (factor === 0)
|
|
161
|
+
continue;
|
|
162
|
+
for (let c = 0; c < rowR.length; c++) {
|
|
163
|
+
rowR[c] = rowR[c] - factor * pivRow[c];
|
|
164
|
+
}
|
|
165
|
+
S.rhs[r] = S.rhs[r] - factor * S.rhs[bestRow];
|
|
166
|
+
}
|
|
167
|
+
const leaving = basic[bestRow];
|
|
168
|
+
basic[bestRow] = enter;
|
|
169
|
+
nonbasic[col] = leaving;
|
|
170
|
+
return leaving;
|
|
171
|
+
}
|
|
172
|
+
// Selección de pivote inicial: la variable cuyo label es `startLabel`.
|
|
173
|
+
// Esa variable está en `nonbasic` de uno de los dos sistemas.
|
|
174
|
+
let entering;
|
|
175
|
+
if (startLabel < n) {
|
|
176
|
+
entering = { kind: 0, index: startLabel }; // x_{startLabel}
|
|
177
|
+
}
|
|
178
|
+
else {
|
|
179
|
+
entering = { kind: 2, index: startLabel - n }; // y_{startLabel-n}
|
|
180
|
+
}
|
|
181
|
+
// Pivotear hasta que la complementaria del start vuelva a estar en
|
|
182
|
+
// la base, o hasta exceder un budget para evitar loops.
|
|
183
|
+
const maxIter = 4 * (n + m) * (n + m) + 50;
|
|
184
|
+
let currentSystem = entering.kind === 0 ? 'S1' : 'S2';
|
|
185
|
+
for (let iter = 0; iter < maxIter; iter++) {
|
|
186
|
+
const leaving = pivotIn(currentSystem, entering);
|
|
187
|
+
if (leaving === null)
|
|
188
|
+
return null; // unbounded → no Nash recovered
|
|
189
|
+
// Si la saliente es la complementaria de la variable que originalmente
|
|
190
|
+
// soltamos (startLabel), terminamos.
|
|
191
|
+
if (labelOf(leaving) === startLabel)
|
|
192
|
+
break;
|
|
193
|
+
// Caso contrario: la nueva entrante es la complementaria de la saliente,
|
|
194
|
+
// en el sistema opuesto.
|
|
195
|
+
entering = complementOf(leaving);
|
|
196
|
+
currentSystem = entering.kind === 0 || entering.kind === 1 ? 'S1' : 'S2';
|
|
197
|
+
// Pero entrar como x o como r requiere que la variable sea no-básica
|
|
198
|
+
// en su sistema. x_i / r_i son complementarias y ambas viven asociadas
|
|
199
|
+
// a S1. Igual y/s respecto a S2.
|
|
200
|
+
if (currentSystem === 'S1') {
|
|
201
|
+
// las no-básicas de S1 son x_*, las básicas son s_* o r_* o lo que
|
|
202
|
+
// hayamos pivoteado. Si la variable a entrar es r_i, debe entrar a
|
|
203
|
+
// S1 si está no-básica allí.
|
|
204
|
+
currentSystem = 'S1';
|
|
205
|
+
}
|
|
206
|
+
else {
|
|
207
|
+
currentSystem = 'S2';
|
|
208
|
+
}
|
|
209
|
+
// Si la entrante no está actualmente no-básica en S{1,2}, intercambiar
|
|
210
|
+
// sistemas.
|
|
211
|
+
const candidateSet = currentSystem === 'S1' ? S1_nonbasic : S2_nonbasic;
|
|
212
|
+
let found = false;
|
|
213
|
+
for (const v of candidateSet) {
|
|
214
|
+
if (v.kind === entering.kind && v.index === entering.index) {
|
|
215
|
+
found = true;
|
|
216
|
+
break;
|
|
217
|
+
}
|
|
218
|
+
}
|
|
219
|
+
if (!found) {
|
|
220
|
+
currentSystem = currentSystem === 'S1' ? 'S2' : 'S1';
|
|
221
|
+
}
|
|
222
|
+
}
|
|
223
|
+
// Extraer x e y de las básicas
|
|
224
|
+
const xVec = new Array(n).fill(0);
|
|
225
|
+
const yVec = new Array(m).fill(0);
|
|
226
|
+
for (let r = 0; r < S1_basic.length; r++) {
|
|
227
|
+
const v = S1_basic[r];
|
|
228
|
+
if (v.kind === 0)
|
|
229
|
+
xVec[v.index] = S1.rhs[r]; // x_i
|
|
230
|
+
}
|
|
231
|
+
for (let r = 0; r < S2_basic.length; r++) {
|
|
232
|
+
const v = S2_basic[r];
|
|
233
|
+
if (v.kind === 2)
|
|
234
|
+
yVec[v.index] = S2.rhs[r]; // y_j
|
|
235
|
+
}
|
|
236
|
+
// Si todo terminó en (0, 0) volvimos al inicial → no encontró Nash.
|
|
237
|
+
const sumX = xVec.reduce((s, v) => s + v, 0);
|
|
238
|
+
const sumY = yVec.reduce((s, v) => s + v, 0);
|
|
239
|
+
if (sumX < TOL || sumY < TOL)
|
|
240
|
+
return null;
|
|
241
|
+
// Normalizar a distribuciones de probabilidad
|
|
242
|
+
const p = xVec.map((v) => v / sumX);
|
|
243
|
+
const q = yVec.map((v) => v / sumY);
|
|
244
|
+
// Clean up negativos por ruido numérico
|
|
245
|
+
for (let i = 0; i < p.length; i++)
|
|
246
|
+
if (p[i] < 0 && p[i] > -TOL)
|
|
247
|
+
p[i] = 0;
|
|
248
|
+
for (let j = 0; j < q.length; j++)
|
|
249
|
+
if (q[j] < 0 && q[j] > -TOL)
|
|
250
|
+
q[j] = 0;
|
|
251
|
+
const dists = [p, q];
|
|
252
|
+
const payoffs = [
|
|
253
|
+
(0, types_1.expectedPayoffFromDistributions)(game, 0, dists),
|
|
254
|
+
(0, types_1.expectedPayoffFromDistributions)(game, 1, dists),
|
|
255
|
+
];
|
|
256
|
+
const supportP = p.filter((x) => x > TOL).length;
|
|
257
|
+
const supportQ = q.filter((x) => x > TOL).length;
|
|
258
|
+
const isPure = supportP === 1 && supportQ === 1;
|
|
259
|
+
return {
|
|
260
|
+
strategies: [
|
|
261
|
+
{ player: 0, distribution: p },
|
|
262
|
+
{ player: 1, distribution: q },
|
|
263
|
+
],
|
|
264
|
+
payoffs,
|
|
265
|
+
isPure,
|
|
266
|
+
isStrict: false,
|
|
267
|
+
};
|
|
268
|
+
}
|
|
269
|
+
//# sourceMappingURL=lemke-howson.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"lemke-howson.js","sourceRoot":"","sources":["../../src/game-theory/lemke-howson.ts"],"names":[],"mappings":";AAAA,+DAA+D;AAC/D,gEAAgE;AAChE,+DAA+D;AAC/D,EAAE;AACF,sCAAsC;AACtC,EAAE;AACF,+DAA+D;AAC/D,8DAA8D;AAC9D,gEAAgE;AAChE,mCAAmC;AACnC,2CAA2C;AAC3C,2CAA2C;AAC3C,+DAA+D;AAC/D,oCAAoC;AACpC,4DAA4D;AAC5D,iEAAiE;AACjE,oEAAoE;AACpE,kEAAkE;AAClE,gEAAgE;AAChE,kEAAkE;AAClE,6CAA6C;AAC7C,EAAE;AACF,oEAAoE;AACpE,4DAA4D;AAC5D,EAAE;AACF,yDAAyD;AACzD,qEAAqE;AACrE,mEAAmE;AACnE,kCAAkC;;AAOlC,kCAmQC;AAvQD,mCAA0D;AAE1D,MAAM,GAAG,GAAG,IAAI,CAAC;AAEjB,SAAgB,WAAW,CAAC,IAAoB,EAAE,UAAU,GAAG,CAAC;IAC9D,IAAI,IAAI,CAAC,OAAO,KAAK,CAAC,EAAE,CAAC;QACvB,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,CAAC;IACD,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IAC7B,IAAI,UAAU,GAAG,CAAC,IAAI,UAAU,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QAC1C,MAAM,IAAI,KAAK,CAAC,cAAc,UAAU,qBAAqB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACzE,CAAC;IACD,MAAM,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;IAC3B,MAAM,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;IAE3B,mDAAmD;IACnD,MAAM,CAAC,GAAe,EAAE,CAAC;IACzB,MAAM,CAAC,GAAe,EAAE,CAAC;IACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,IAAI,GAAa,EAAE,CAAC;QAC1B,MAAM,IAAI,GAAa,EAAE,CAAC;QAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACzB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC3B,CAAC;QACD,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACb,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACf,CAAC;IAED,gCAAgC;IAChC,IAAI,MAAM,GAAG,QAAQ,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnB,IAAI,EAAE,GAAG,MAAM;gBAAE,MAAM,GAAG,EAAE,CAAC;YAC7B,IAAI,EAAE,GAAG,MAAM;gBAAE,MAAM,GAAG,EAAE,CAAC;QAC/B,CAAC;IACH,CAAC;IACD,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3C,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;gBAC1B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;YAC5B,CAAC;QACH,CAAC;IACH,CAAC;IAoCD,MAAM,OAAO,GAAG,CAAC,CAAM,EAAU,EAAE;QACjC,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;YAAE,OAAO,CAAC,CAAC,KAAK,CAAC;QACjD,OAAO,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC;IACrB,CAAC,CAAC;IACF,MAAM,YAAY,GAAG,CAAC,CAAM,EAAO,EAAE;QACnC,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;YAAE,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,KAAK,EAAE,CAAC;QACrD,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;YAAE,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,KAAK,EAAE,CAAC;QACrD,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;YAAE,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,KAAK,EAAE,CAAC;QACrD,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,KAAK,EAAE,CAAC;IACrC,CAAC,CAAC;IAEF,oEAAoE;IACpE,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,EAAgB,EAAE,GAAG,EAAE,EAAc,EAAE,CAAC;IAC1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAa,EAAE,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9C,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACjB,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACjB,CAAC;IACD,MAAM,QAAQ,GAAU,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAQ,CAAC,CAAC;IAC5F,MAAM,WAAW,GAAU,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAQ,CAAC,CAAC;IAE/F,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,EAAgB,EAAE,GAAG,EAAE,EAAc,EAAE,CAAC;IAC1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAa,EAAE,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9C,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACjB,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACjB,CAAC;IACD,MAAM,QAAQ,GAAU,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAQ,CAAC,CAAC;IAC5F,MAAM,WAAW,GAAU,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAQ,CAAC,CAAC;IAE/F,mEAAmE;IACnE,iEAAiE;IACjE,iEAAiE;IACjE,4DAA4D;IAC5D,SAAS,OAAO,CAAC,GAAgB,EAAE,KAAU;QAC3C,MAAM,CAAC,GAAG,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACjC,MAAM,QAAQ,GAAG,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,WAAW,CAAC;QAC1D,MAAM,KAAK,GAAG,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC;QAEjD,mCAAmC;QACnC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACzC,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;YACtB,IAAI,CAAC,CAAC,IAAI,KAAK,KAAK,CAAC,IAAI,IAAI,CAAC,CAAC,KAAK,KAAK,KAAK,CAAC,KAAK,EAAE,CAAC;gBACrD,GAAG,GAAG,CAAC,CAAC;gBACR,MAAM;YACR,CAAC;QACH,CAAC;QACD,IAAI,GAAG,GAAG,CAAC;YAAE,OAAO,IAAI,CAAC;QAEzB,uDAAuD;QACvD,IAAI,OAAO,GAAG,CAAC,CAAC,CAAC;QACjB,IAAI,SAAS,GAAG,QAAQ,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACtC,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAC3B,IAAI,IAAI,IAAI,GAAG;gBAAE,SAAS;YAC1B,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;YAC9B,IAAI,KAAK,GAAG,SAAS,GAAG,GAAG,EAAE,CAAC;gBAC5B,SAAS,GAAG,KAAK,CAAC;gBAClB,OAAO,GAAG,CAAC,CAAC;YACd,CAAC;QACH,CAAC;QACD,IAAI,OAAO,GAAG,CAAC;YAAE,OAAO,IAAI,CAAC,CAAC,YAAY;QAE1C,uCAAuC;QACvC,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC;QAClC,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtE,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,KAAK,CAAC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACtC,IAAI,CAAC,KAAK,OAAO;gBAAE,SAAS;YAC5B,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACtB,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;YACzB,IAAI,MAAM,KAAK,CAAC;gBAAE,SAAS;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBACrC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACzC,CAAC;YACD,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAChD,CAAC;QAED,MAAM,OAAO,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC;QAC/B,KAAK,CAAC,OAAO,CAAC,GAAG,KAAK,CAAC;QACvB,QAAQ,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC;QACxB,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,uEAAuE;IACvE,8DAA8D;IAC9D,IAAI,QAAa,CAAC;IAClB,IAAI,UAAU,GAAG,CAAC,EAAE,CAAC;QACnB,QAAQ,GAAG,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,UAAU,EAAE,CAAC,CAAC,iBAAiB;IAC9D,CAAC;SAAM,CAAC;QACN,QAAQ,GAAG,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,UAAU,GAAG,CAAC,EAAE,CAAC,CAAC,mBAAmB;IACpE,CAAC;IAED,mEAAmE;IACnE,wDAAwD;IACxD,MAAM,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;IAC3C,IAAI,aAAa,GAAgB,QAAQ,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC;IACnE,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,OAAO,EAAE,IAAI,EAAE,EAAE,CAAC;QAC1C,MAAM,OAAO,GAAG,OAAO,CAAC,aAAa,EAAE,QAAQ,CAAC,CAAC;QACjD,IAAI,OAAO,KAAK,IAAI;YAAE,OAAO,IAAI,CAAC,CAAC,gCAAgC;QACnE,uEAAuE;QACvE,qCAAqC;QACrC,IAAI,OAAO,CAAC,OAAO,CAAC,KAAK,UAAU;YAAE,MAAM;QAC3C,yEAAyE;QACzE,yBAAyB;QACzB,QAAQ,GAAG,YAAY,CAAC,OAAO,CAAC,CAAC;QACjC,aAAa,GAAG,QAAQ,CAAC,IAAI,KAAK,CAAC,IAAI,QAAQ,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC;QACzE,qEAAqE;QACrE,uEAAuE;QACvE,iCAAiC;QACjC,IAAI,aAAa,KAAK,IAAI,EAAE,CAAC;YAC3B,mEAAmE;YACnE,mEAAmE;YACnE,6BAA6B;YAC7B,aAAa,GAAG,IAAI,CAAC;QACvB,CAAC;aAAM,CAAC;YACN,aAAa,GAAG,IAAI,CAAC;QACvB,CAAC;QACD,uEAAuE;QACvE,YAAY;QACZ,MAAM,YAAY,GAAG,aAAa,KAAK,IAAI,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,WAAW,CAAC;QACxE,IAAI,KAAK,GAAG,KAAK,CAAC;QAClB,KAAK,MAAM,CAAC,IAAI,YAAY,EAAE,CAAC;YAC7B,IAAI,CAAC,CAAC,IAAI,KAAK,QAAQ,CAAC,IAAI,IAAI,CAAC,CAAC,KAAK,KAAK,QAAQ,CAAC,KAAK,EAAE,CAAC;gBAC3D,KAAK,GAAG,IAAI,CAAC;gBACb,MAAM;YACR,CAAC;QACH,CAAC;QACD,IAAI,CAAC,KAAK,EAAE,CAAC;YACX,aAAa,GAAG,aAAa,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC;QACvD,CAAC;IACH,CAAC;IAED,+BAA+B;IAC/B,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACtB,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;YAAE,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM;IACrD,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACtB,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;YAAE,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM;IACrD,CAAC;IAED,oEAAoE;IACpE,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7C,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7C,IAAI,IAAI,GAAG,GAAG,IAAI,IAAI,GAAG,GAAG;QAAE,OAAO,IAAI,CAAC;IAE1C,8CAA8C;IAC9C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IACpC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IAEpC,wCAAwC;IACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG;YAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACzE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG;YAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAEzE,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrB,MAAM,OAAO,GAAG;QACd,IAAA,uCAA+B,EAAC,IAAI,EAAE,CAAC,EAAE,KAAK,CAAC;QAC/C,IAAA,uCAA+B,EAAC,IAAI,EAAE,CAAC,EAAE,KAAK,CAAC;KAChD,CAAC;IACF,MAAM,QAAQ,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACjD,MAAM,QAAQ,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACjD,MAAM,MAAM,GAAG,QAAQ,KAAK,CAAC,IAAI,QAAQ,KAAK,CAAC,CAAC;IAChD,OAAO;QACL,UAAU,EAAE;YACV,EAAE,MAAM,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE;YAC9B,EAAE,MAAM,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE;SAC/B;QACD,OAAO;QACP,MAAM;QACN,QAAQ,EAAE,KAAK;KAChB,CAAC;AACJ,CAAC"}
|
|
@@ -0,0 +1,7 @@
|
|
|
1
|
+
/**
|
|
2
|
+
* Resuelve A x = b por eliminación gaussiana con pivot parcial.
|
|
3
|
+
* Devuelve null si el sistema es singular (no único o sin solución).
|
|
4
|
+
* Diseñado para matrices chicas (n ≤ 20).
|
|
5
|
+
*/
|
|
6
|
+
export declare function solveLinear(A: number[][], b: number[], tol?: number): number[] | null;
|
|
7
|
+
//# sourceMappingURL=linalg.d.ts.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"linalg.d.ts","sourceRoot":"","sources":["../../src/game-theory/linalg.ts"],"names":[],"mappings":"AAMA;;;;GAIG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,SAAQ,GAAG,MAAM,EAAE,GAAG,IAAI,CAmDpF"}
|
|
@@ -0,0 +1,69 @@
|
|
|
1
|
+
"use strict";
|
|
2
|
+
// ============================================================
|
|
3
|
+
// Mini álgebra lineal (sin dependencias) — para resolver los
|
|
4
|
+
// sistemas indiferentes del support enumeration y los pivots de
|
|
5
|
+
// Lemke-Howson.
|
|
6
|
+
// ============================================================
|
|
7
|
+
Object.defineProperty(exports, "__esModule", { value: true });
|
|
8
|
+
exports.solveLinear = solveLinear;
|
|
9
|
+
/**
|
|
10
|
+
* Resuelve A x = b por eliminación gaussiana con pivot parcial.
|
|
11
|
+
* Devuelve null si el sistema es singular (no único o sin solución).
|
|
12
|
+
* Diseñado para matrices chicas (n ≤ 20).
|
|
13
|
+
*/
|
|
14
|
+
function solveLinear(A, b, tol = 1e-10) {
|
|
15
|
+
const n = A.length;
|
|
16
|
+
if (n === 0)
|
|
17
|
+
return [];
|
|
18
|
+
for (const row of A) {
|
|
19
|
+
if (row.length !== n)
|
|
20
|
+
throw new Error('solveLinear: matrix must be square');
|
|
21
|
+
}
|
|
22
|
+
if (b.length !== n)
|
|
23
|
+
throw new Error('solveLinear: b length mismatch');
|
|
24
|
+
// Matriz aumentada
|
|
25
|
+
const M = A.map((row, i) => {
|
|
26
|
+
const r = row.slice();
|
|
27
|
+
r.push(b[i]);
|
|
28
|
+
return r;
|
|
29
|
+
});
|
|
30
|
+
for (let col = 0; col < n; col++) {
|
|
31
|
+
// pivot parcial
|
|
32
|
+
let pivotRow = col;
|
|
33
|
+
let pivotVal = Math.abs(M[col][col]);
|
|
34
|
+
for (let r = col + 1; r < n; r++) {
|
|
35
|
+
const v = Math.abs(M[r][col]);
|
|
36
|
+
if (v > pivotVal) {
|
|
37
|
+
pivotVal = v;
|
|
38
|
+
pivotRow = r;
|
|
39
|
+
}
|
|
40
|
+
}
|
|
41
|
+
if (pivotVal < tol)
|
|
42
|
+
return null; // singular
|
|
43
|
+
if (pivotRow !== col) {
|
|
44
|
+
const tmp = M[col];
|
|
45
|
+
M[col] = M[pivotRow];
|
|
46
|
+
M[pivotRow] = tmp;
|
|
47
|
+
}
|
|
48
|
+
const pivot = M[col][col];
|
|
49
|
+
const pivRow = M[col];
|
|
50
|
+
for (let r = 0; r < n; r++) {
|
|
51
|
+
if (r === col)
|
|
52
|
+
continue;
|
|
53
|
+
const rowR = M[r];
|
|
54
|
+
const factor = rowR[col] / pivot;
|
|
55
|
+
if (factor === 0)
|
|
56
|
+
continue;
|
|
57
|
+
for (let c = col; c <= n; c++) {
|
|
58
|
+
rowR[c] = rowR[c] - factor * pivRow[c];
|
|
59
|
+
}
|
|
60
|
+
}
|
|
61
|
+
}
|
|
62
|
+
const x = new Array(n);
|
|
63
|
+
for (let i = 0; i < n; i++) {
|
|
64
|
+
const row = M[i];
|
|
65
|
+
x[i] = row[n] / row[i];
|
|
66
|
+
}
|
|
67
|
+
return x;
|
|
68
|
+
}
|
|
69
|
+
//# sourceMappingURL=linalg.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"linalg.js","sourceRoot":"","sources":["../../src/game-theory/linalg.ts"],"names":[],"mappings":";AAAA,+DAA+D;AAC/D,6DAA6D;AAC7D,gEAAgE;AAChE,gBAAgB;AAChB,+DAA+D;;AAO/D,kCAmDC;AAxDD;;;;GAIG;AACH,SAAgB,WAAW,CAAC,CAAa,EAAE,CAAW,EAAE,GAAG,GAAG,KAAK;IACjE,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,EAAE,CAAC;IACvB,KAAK,MAAM,GAAG,IAAI,CAAC,EAAE,CAAC;QACpB,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;IAC9E,CAAC;IACD,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;IAEtE,mBAAmB;IACnB,MAAM,CAAC,GAAe,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE;QACrC,MAAM,CAAC,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QACtB,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACb,OAAO,CAAC,CAAC;IACX,CAAC,CAAC,CAAC;IAEH,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,gBAAgB;QAChB,IAAI,QAAQ,GAAG,GAAG,CAAC;QACnB,IAAI,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC9B,IAAI,CAAC,GAAG,QAAQ,EAAE,CAAC;gBACjB,QAAQ,GAAG,CAAC,CAAC;gBACb,QAAQ,GAAG,CAAC,CAAC;YACf,CAAC;QACH,CAAC;QACD,IAAI,QAAQ,GAAG,GAAG;YAAE,OAAO,IAAI,CAAC,CAAC,WAAW;QAC5C,IAAI,QAAQ,KAAK,GAAG,EAAE,CAAC;YACrB,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YACnB,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC;YACrB,CAAC,CAAC,QAAQ,CAAC,GAAG,GAAG,CAAC;QACpB,CAAC;QACD,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC1B,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,KAAK,GAAG;gBAAE,SAAS;YACxB,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAClB,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;YACjC,IAAI,MAAM,KAAK,CAAC;gBAAE,SAAS;YAC3B,KAAK,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACzC,CAAC;QACH,CAAC;IACH,CAAC;IAED,MAAM,CAAC,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC;IAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjB,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"pure-nash.d.ts","sourceRoot":"","sources":["../../src/game-theory/pure-nash.ts"],"names":[],"mappings":"AAQA,OAAO,KAAK,EAAE,eAAe,EAAE,cAAc,EAAE,MAAM,SAAS,CAAC;AAS/D,wBAAgB,YAAY,CAAC,IAAI,EAAE,cAAc,GAAG,eAAe,EAAE,CAYpE"}
|
|
@@ -0,0 +1,86 @@
|
|
|
1
|
+
"use strict";
|
|
2
|
+
// ============================================================
|
|
3
|
+
// Pure-strategy Nash equilibria
|
|
4
|
+
// ============================================================
|
|
5
|
+
//
|
|
6
|
+
// Un perfil puro (s_1, ..., s_n) es Nash sii cada s_i es best
|
|
7
|
+
// response del jugador i contra (s_{-i}). Lo verificamos
|
|
8
|
+
// enumerando todos los perfiles y comprobando la condición.
|
|
9
|
+
Object.defineProperty(exports, "__esModule", { value: true });
|
|
10
|
+
exports.findPureNash = findPureNash;
|
|
11
|
+
const types_1 = require("./types");
|
|
12
|
+
function findPureNash(game) {
|
|
13
|
+
const sizes = game.strategies;
|
|
14
|
+
const total = (0, types_1.totalProfiles)(sizes);
|
|
15
|
+
const equilibria = [];
|
|
16
|
+
for (let idx = 0; idx < total; idx++) {
|
|
17
|
+
const profile = (0, types_1.indexToProfile)(idx, sizes);
|
|
18
|
+
if (isProfileNash(game, profile)) {
|
|
19
|
+
equilibria.push(buildPureEquilibrium(game, profile));
|
|
20
|
+
}
|
|
21
|
+
}
|
|
22
|
+
return equilibria;
|
|
23
|
+
}
|
|
24
|
+
function isProfileNash(game, profile) {
|
|
25
|
+
const sizes = game.strategies;
|
|
26
|
+
for (let p = 0; p < game.players; p++) {
|
|
27
|
+
const mySize = sizes[p];
|
|
28
|
+
const baseIdx = (0, types_1.profileToIndex)(profile, sizes);
|
|
29
|
+
const row = game.payoffs[p];
|
|
30
|
+
if (!row)
|
|
31
|
+
throw new Error(`payoff missing for player ${p}`);
|
|
32
|
+
const curU = row[baseIdx];
|
|
33
|
+
for (let alt = 0; alt < mySize; alt++) {
|
|
34
|
+
if (alt === profile[p])
|
|
35
|
+
continue;
|
|
36
|
+
const altProfile = profile.slice();
|
|
37
|
+
altProfile[p] = alt;
|
|
38
|
+
const altIdx = (0, types_1.profileToIndex)(altProfile, sizes);
|
|
39
|
+
const altU = row[altIdx];
|
|
40
|
+
if (altU > curU + 1e-12)
|
|
41
|
+
return false;
|
|
42
|
+
}
|
|
43
|
+
}
|
|
44
|
+
return true;
|
|
45
|
+
}
|
|
46
|
+
function isStrictPureNash(game, profile) {
|
|
47
|
+
const sizes = game.strategies;
|
|
48
|
+
for (let p = 0; p < game.players; p++) {
|
|
49
|
+
const mySize = sizes[p];
|
|
50
|
+
const baseIdx = (0, types_1.profileToIndex)(profile, sizes);
|
|
51
|
+
const row = game.payoffs[p];
|
|
52
|
+
if (!row)
|
|
53
|
+
throw new Error(`payoff missing for player ${p}`);
|
|
54
|
+
const curU = row[baseIdx];
|
|
55
|
+
for (let alt = 0; alt < mySize; alt++) {
|
|
56
|
+
if (alt === profile[p])
|
|
57
|
+
continue;
|
|
58
|
+
const altProfile = profile.slice();
|
|
59
|
+
altProfile[p] = alt;
|
|
60
|
+
const altIdx = (0, types_1.profileToIndex)(altProfile, sizes);
|
|
61
|
+
const altU = row[altIdx];
|
|
62
|
+
if (altU >= curU - 1e-12)
|
|
63
|
+
return false;
|
|
64
|
+
}
|
|
65
|
+
}
|
|
66
|
+
return true;
|
|
67
|
+
}
|
|
68
|
+
function buildPureEquilibrium(game, profile) {
|
|
69
|
+
const sizes = game.strategies;
|
|
70
|
+
const strategies = profile.map((s, p) => ({
|
|
71
|
+
player: p,
|
|
72
|
+
distribution: (0, types_1.pureDistribution)(s, sizes[p]),
|
|
73
|
+
}));
|
|
74
|
+
const dists = strategies.map((s) => s.distribution);
|
|
75
|
+
const payoffs = new Array(game.players);
|
|
76
|
+
for (let p = 0; p < game.players; p++) {
|
|
77
|
+
payoffs[p] = (0, types_1.expectedPayoffFromDistributions)(game, p, dists);
|
|
78
|
+
}
|
|
79
|
+
return {
|
|
80
|
+
strategies,
|
|
81
|
+
payoffs,
|
|
82
|
+
isPure: true,
|
|
83
|
+
isStrict: isStrictPureNash(game, profile),
|
|
84
|
+
};
|
|
85
|
+
}
|
|
86
|
+
//# sourceMappingURL=pure-nash.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"pure-nash.js","sourceRoot":"","sources":["../../src/game-theory/pure-nash.ts"],"names":[],"mappings":";AAAA,+DAA+D;AAC/D,gCAAgC;AAChC,+DAA+D;AAC/D,EAAE;AACF,8DAA8D;AAC9D,yDAAyD;AACzD,4DAA4D;;AAW5D,oCAYC;AApBD,mCAMiB;AAEjB,SAAgB,YAAY,CAAC,IAAoB;IAC/C,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC;IAC9B,MAAM,KAAK,GAAG,IAAA,qBAAa,EAAC,KAAK,CAAC,CAAC;IACnC,MAAM,UAAU,GAAsB,EAAE,CAAC;IAEzC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,KAAK,EAAE,GAAG,EAAE,EAAE,CAAC;QACrC,MAAM,OAAO,GAAG,IAAA,sBAAc,EAAC,GAAG,EAAE,KAAK,CAAC,CAAC;QAC3C,IAAI,aAAa,CAAC,IAAI,EAAE,OAAO,CAAC,EAAE,CAAC;YACjC,UAAU,CAAC,IAAI,CAAC,oBAAoB,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC;QACvD,CAAC;IACH,CAAC;IACD,OAAO,UAAU,CAAC;AACpB,CAAC;AAED,SAAS,aAAa,CAAC,IAAoB,EAAE,OAAiB;IAC5D,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC;IAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,MAAM,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,OAAO,GAAG,IAAA,sBAAc,EAAC,OAAO,EAAE,KAAK,CAAC,CAAC;QAC/C,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;QAC5B,IAAI,CAAC,GAAG;YAAE,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,EAAE,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,GAAG,CAAC,OAAO,CAAC,CAAC;QAC1B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC;YACtC,IAAI,GAAG,KAAK,OAAO,CAAC,CAAC,CAAC;gBAAE,SAAS;YACjC,MAAM,UAAU,GAAG,OAAO,CAAC,KAAK,EAAE,CAAC;YACnC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YACpB,MAAM,MAAM,GAAG,IAAA,sBAAc,EAAC,UAAU,EAAE,KAAK,CAAC,CAAC;YACjD,MAAM,IAAI,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC;YACzB,IAAI,IAAI,GAAG,IAAI,GAAG,KAAK;gBAAE,OAAO,KAAK,CAAC;QACxC,CAAC;IACH,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,gBAAgB,CAAC,IAAoB,EAAE,OAAiB;IAC/D,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC;IAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,MAAM,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,OAAO,GAAG,IAAA,sBAAc,EAAC,OAAO,EAAE,KAAK,CAAC,CAAC;QAC/C,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;QAC5B,IAAI,CAAC,GAAG;YAAE,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,EAAE,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,GAAG,CAAC,OAAO,CAAC,CAAC;QAC1B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC;YACtC,IAAI,GAAG,KAAK,OAAO,CAAC,CAAC,CAAC;gBAAE,SAAS;YACjC,MAAM,UAAU,GAAG,OAAO,CAAC,KAAK,EAAE,CAAC;YACnC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YACpB,MAAM,MAAM,GAAG,IAAA,sBAAc,EAAC,UAAU,EAAE,KAAK,CAAC,CAAC;YACjD,MAAM,IAAI,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC;YACzB,IAAI,IAAI,IAAI,IAAI,GAAG,KAAK;gBAAE,OAAO,KAAK,CAAC;QACzC,CAAC;IACH,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,oBAAoB,CAAC,IAAoB,EAAE,OAAiB;IACnE,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC;IAC9B,MAAM,UAAU,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACxC,MAAM,EAAE,CAAC;QACT,YAAY,EAAE,IAAA,wBAAgB,EAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;KAC5C,CAAC,CAAC,CAAC;IACJ,MAAM,KAAK,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC;IACpD,MAAM,OAAO,GAAG,IAAI,KAAK,CAAS,IAAI,CAAC,OAAO,CAAC,CAAC;IAChD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,OAAO,CAAC,CAAC,CAAC,GAAG,IAAA,uCAA+B,EAAC,IAAI,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IAC/D,CAAC;IACD,OAAO;QACL,UAAU;QACV,OAAO;QACP,MAAM,EAAE,IAAI;QACZ,QAAQ,EAAE,gBAAgB,CAAC,IAAI,EAAE,OAAO,CAAC;KAC1C,CAAC;AACJ,CAAC"}
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"support-enumeration.d.ts","sourceRoot":"","sources":["../../src/game-theory/support-enumeration.ts"],"names":[],"mappings":"AAsBA,OAAO,KAAK,EAAE,eAAe,EAAE,cAAc,EAAE,MAAM,SAAS,CAAC;AAM/D,wBAAgB,gBAAgB,CAAC,IAAI,EAAE,cAAc,EAAE,OAAO,CAAC,EAAE,MAAM,GAAG,eAAe,EAAE,CA8B1F"}
|