diff --git a/reasoning_core/tasks/generated/matrix_induction.py b/reasoning_core/tasks/generated/matrix_induction.py
new file mode 100644
index 0000000..0bee7f7
--- /dev/null
+++ b/reasoning_core/tasks/generated/matrix_induction.py
@@ -0,0 +1,164 @@
+import random
+import re
+from dataclasses import dataclass
+
+from reasoning_core.template import Config, Entry, Task, edict, stochastic_rounding as sround
+from ._base import GeneratedMixin
+
+
+ATTRIBUTES = {
+    "shape": ("circle", "triangle", "square", "diamond"),
+    "count": ("1", "2", "3", "4"),
+    "orientation": ("north", "east", "south", "west"),
+    "fill": ("open", "filled"),
+}
+
+
+def _rules(k):
+    rules = [
+        ("left", lambda a, b, k=k: a),
+        ("right", lambda a, b, k=k: b),
+        ("min", lambda a, b, k=k: min(a, b)),
+        ("max", lambda a, b, k=k: max(a, b)),
+    ]
+    rules += [(f"add+{c}", lambda a, b, c=c, k=k: (a + b + c) % k) for c in range(k)]
+    if k & (k - 1) == 0:
+        rules += [(f"xor+{c}", lambda a, b, c=c, k=k: a ^ b ^ c) for c in range(k)]
+    return rules
+
+
+def _satisfies(grid, rule):
+    return all(rule(row[0], row[1]) == row[2] for row in grid) and all(
+        rule(grid[0][j], grid[1][j]) == grid[2][j] for j in range(3)
+    )
+
+
+def _complete_values(partial, missing, k):
+    values = set()
+    survivors = {}
+    for value in range(k):
+        grid = [row[:] for row in partial]
+        grid[missing[0]][missing[1]] = value
+        names = [name for name, rule in _rules(k) if _satisfies(grid, rule)]
+        if names:
+            values.add(value)
+            survivors[value] = names
+    return values, survivors
+
+
+def _make_grid(rule, k):
+    grid = [[None] * 3 for _ in range(3)]
+    for i in range(2):
+        for j in range(2):
+            grid[i][j] = random.randrange(k)
+    for i in range(2):
+        grid[i][2] = rule(grid[i][0], grid[i][1])
+    for j in range(2):
+        grid[2][j] = rule(grid[0][j], grid[1][j])
+    grid[2][2] = rule(grid[2][0], grid[2][1])
+    return grid if rule(grid[0][2], grid[1][2]) == grid[2][2] else None
+
+
+def _cell_text(values, attrs):
+    return " ".join(f"{name}={ATTRIBUTES[name][values[name]]}" for name in attrs)
+
+
+@dataclass
+class MatrixInductionConfig(Config):
+    n_attributes: int = 1
+    hard_position_p: float = 0.15
+    hard_rule_p: float = 0.35
+
+    def apply_difficulty(self, level):
+        self.n_attributes = sround(self.n_attributes + 0.6 * level)
+        self.hard_position_p = min(0.9, self.hard_position_p + 0.12 * level)
+        self.hard_rule_p = min(0.9, self.hard_rule_p + 0.1 * level)
+
+
+class MatrixInduction(GeneratedMixin, Task):
+    summary = "Infer a missing multi-attribute matrix cell under a certified finite rule family."
+    config_cls = MatrixInductionConfig
+
+    def generate_entry(self):
+        n_attributes = min(4, max(1, int(self.config.n_attributes)))
+        attrs = random.sample(list(ATTRIBUTES), n_attributes)
+
+        for _ in range(300):
+            missing = random.choice([(0, 0), (0, 1), (1, 0), (1, 1)]) if random.random() < self.config.hard_position_p else (2, 2)
+            grids = {}
+            hidden_rules = {}
+            survivor_counts = {}
+            ok = True
+
+            for attr in attrs:
+                k = len(ATTRIBUTES[attr])
+                bank = _rules(k)
+                hard = [x for x in bank if x[0].startswith(("add", "xor"))]
+                name, rule = random.choice(hard if hard and random.random() < self.config.hard_rule_p else bank)
+                grid = _make_grid(rule, k)
+                if grid is None or len({x for row in grid for x in row}) < 2:
+                    ok = False
+                    break
+                partial = [row[:] for row in grid]
+                partial[missing[0]][missing[1]] = None
+                values, survivors = _complete_values(partial, missing, k)
+                if values != {grid[missing[0]][missing[1]]}:
+                    ok = False
+                    break
+                grids[attr] = grid
+                hidden_rules[attr] = name
+                survivor_counts[attr] = len(survivors[grid[missing[0]][missing[1]]])
+
+            if not ok:
+                continue
+
+            cells = []
+            for i in range(3):
+                row = []
+                for j in range(3):
+                    if (i, j) == missing:
+                        row.append("?")
+                    else:
+                        row.append(_cell_text({a: grids[a][i][j] for a in attrs}, attrs))
+                cells.append(row)
+
+            gold_values = {a: grids[a][missing[0]][missing[1]] for a in attrs}
+            answer = _cell_text(gold_values, attrs)
+            return Entry(
+                metadata=edict(
+                    attributes=attrs,
+                    cells=cells,
+                    missing=missing,
+                    hidden_rules=hidden_rules,
+                    survivor_counts=survivor_counts,
+                ),
+                answer=answer,
+            )
+        raise RuntimeError("MatrixInduction: could not build a uniquely completable instance")
+
+    def render_prompt(self, m):
+        domains = "\n".join(f"- {a}: " + ", ".join(ATTRIBUTES[a]) for a in m.attributes)
+        rows = "\n".join(" | ".join(row) for row in m.cells)
+        return (
+            "Complete the missing cell of the 3x3 matrix. Each attribute is independent and uses one fixed rule "
+            "for every row and every column. Encode each listed domain by indices 0,1,... in the shown order. "
+            "For a row or column with encoded values a,b,c, c is obtained from a,b by one of: left=a; right=b; "
+            "min=min(a,b); max=max(a,b); add+t=(a+b+t) mod k for some t; xor+t=a xor b xor t for some t "
+            "(xor is used only for power-of-two domain sizes). Different attributes may use different rules.\n"
+            f"Domains:\n{domains}\nMatrix:\n{rows}\n"
+            "The answer is the missing cell written with exactly the displayed attribute names as name=value pairs."
+        )
+
+    def score_answer(self, answer, entry):
+        def parse(text):
+            out = {}
+            for name in entry.metadata["attributes"]:
+                m = re.search(rf"\b{name}\s*=\s*([A-Za-z0-9_-]+)", str(text), re.I)
+                if not m:
+                    return None
+                out[name] = m.group(1).lower()
+            return out
+        return float(parse(answer) == parse(entry.answer))
+
+    def balancing_key(self, problem):
+        return tuple(problem.metadata.hidden_rules.values())
diff --git a/reasoning_core/tasks/generated/rule_switching.py b/reasoning_core/tasks/generated/rule_switching.py
new file mode 100644
index 0000000..1e87a34
--- /dev/null
+++ b/reasoning_core/tasks/generated/rule_switching.py
@@ -0,0 +1,162 @@
+import random
+from dataclasses import dataclass
+
+from reasoning_core.template import Config, Entry, Task, edict, stochastic_rounding as sround
+from ._base import GeneratedMixin
+
+
+PRIMITIVES = {
+    "rotate-left": (1, 2, 0),
+    "rotate-right": (2, 0, 1),
+    "swap-first-two": (1, 0, 2),
+    "swap-last-two": (0, 2, 1),
+    "swap-outer": (2, 1, 0),
+}
+OPCODES = ("X", "Y", "Z")
+
+
+def _apply(state, regs, perm):
+    old = [state[r] for r in regs]
+    for dst, src in enumerate(perm):
+        state[regs[dst]] = old[src]
+
+
+def _lineage(target, operations):
+    reg = target
+    relevant = []
+    for index in range(len(operations) - 1, -1, -1):
+        op = operations[index]
+        regs = op["regs"]
+        if reg not in regs:
+            continue
+        dst = regs.index(reg)
+        previous = regs[op["perm"][dst]]
+        if previous != reg:
+            relevant.append(index)
+        reg = previous
+    relevant.reverse()
+    return relevant
+
+
+def _has_switch_interference(operations, relevant):
+    by_opcode = {}
+    for i in relevant:
+        op = operations[i]
+        by_opcode.setdefault(op["opcode"], []).append((op["mode"], op["primitive"]))
+    for uses in by_opcode.values():
+        if len({mode for mode, _ in uses}) >= 2 and len({primitive for _, primitive in uses}) >= 2:
+            return True
+    return False
+
+
+@dataclass
+class RuleSwitchingConfig(Config):
+    n_registers: int = 5
+    n_steps: int = 7
+    n_modes: int = 2
+    min_dependency_depth: int = 2
+
+    def apply_difficulty(self, level):
+        self.n_registers = sround(self.n_registers + 0.5 * level)
+        self.n_steps = sround(self.n_steps + 1.8 * level)
+        self.n_modes = sround(self.n_modes + 0.3 * level)
+        self.min_dependency_depth = sround(self.min_dependency_depth + 0.6 * level)
+
+
+class RuleSwitching(GeneratedMixin, Task):
+    summary = "Track symbolic state while identical opcodes change meaning across active rule modes."
+    config_cls = RuleSwitchingConfig
+
+    def generate_entry(self):
+        n_registers = min(9, max(4, int(self.config.n_registers)))
+        n_steps = max(5, int(self.config.n_steps))
+        n_modes = min(4, max(2, int(self.config.n_modes)))
+        min_depth = min(n_steps - 1, max(2, int(self.config.min_dependency_depth)))
+        registers = [f"r{i + 1}" for i in range(n_registers)]
+        values = [chr(ord("A") + i) for i in range(n_registers)]
+
+        for _ in range(400):
+            mappings = []
+            for mode in range(n_modes):
+                for _ in range(50):
+                    chosen = random.sample(list(PRIMITIVES), len(OPCODES))
+                    mapping = dict(zip(OPCODES, chosen))
+                    if mode == 0 or mapping["X"] != mappings[-1]["X"]:
+                        mappings.append(mapping)
+                        break
+
+            state = dict(zip(registers, values))
+            current_mode = 0
+            operations = []
+            program = []
+            switch_every = max(1, n_steps // (n_modes + 1))
+
+            for step in range(n_steps):
+                if step and (step % switch_every == 0 or random.random() < 0.18):
+                    choices = [m for m in range(n_modes) if m != current_mode]
+                    current_mode = random.choice(choices)
+                    program.append({"kind": "mode", "mode": current_mode})
+                opcode = random.choices(OPCODES, weights=[3, 1, 1])[0]
+                regs = tuple(random.sample(registers, 3))
+                primitive = mappings[current_mode][opcode]
+                perm = PRIMITIVES[primitive]
+                operations.append({
+                    "mode": current_mode,
+                    "opcode": opcode,
+                    "regs": regs,
+                    "primitive": primitive,
+                    "perm": perm,
+                })
+                program.append({"kind": "op", "operation": len(operations) - 1})
+                _apply(state, regs, perm)
+
+            candidates = []
+            for target in registers:
+                relevant = _lineage(target, operations)
+                if len(relevant) >= min_depth and _has_switch_interference(operations, relevant):
+                    candidates.append((len(relevant), target, relevant))
+            if not candidates:
+                continue
+            depth, target, relevant = random.choice([x for x in candidates if x[0] == max(y[0] for y in candidates)])
+            return Entry(
+                metadata=edict(
+                    registers=registers,
+                    initial=dict(zip(registers, values)),
+                    mappings=mappings,
+                    program=program,
+                    operations=operations,
+                    target=target,
+                    dependency_depth=depth,
+                    relevant_operations=relevant,
+                ),
+                answer=state[target],
+            )
+        raise RuntimeError("RuleSwitching: could not build a sufficiently interfering trace")
+
+    def render_prompt(self, m):
+        mapping_lines = []
+        for i, mapping in enumerate(m.mappings, 1):
+            mapping_lines.append(f"M{i}: " + ", ".join(f"{op}={mapping[op]}" for op in OPCODES))
+        program_lines = []
+        for item in m.program:
+            if item["kind"] == "mode":
+                program_lines.append(f"mode M{item['mode'] + 1}")
+            else:
+                op = m.operations[item["operation"]]
+                program_lines.append(f"{op['opcode']} " + " ".join(op["regs"]))
+        initial = " ".join(f"{r}={m.initial[r]}" for r in m.registers)
+        return (
+            "Maintain the register values while executing the program. The current mode determines each opcode's "
+            "meaning. For an instruction on registers (a,b,c): rotate-left maps their values to (b,c,a); "
+            "rotate-right to (c,a,b); swap-first-two to (b,a,c); swap-last-two to (a,c,b); and swap-outer to (c,b,a). "
+            "Mode changes affect following instructions.\n"
+            f"Modes:\n" + "\n".join(mapping_lines) + "\n"
+            f"Initial state: {initial}\nStart mode: M1\nProgram:\n" + "\n".join(program_lines) + "\n"
+            f"What value is in {m.target} after the program? The answer is one value label."
+        )
+
+    def score_answer(self, answer, entry):
+        return float(str(answer).strip().upper() == str(entry.answer).strip().upper())
+
+    def balancing_key(self, problem):
+        return problem.answer
diff --git a/reasoning_core/tasks/generated/spatial_folding.py b/reasoning_core/tasks/generated/spatial_folding.py
new file mode 100644
index 0000000..06e3a98
--- /dev/null
+++ b/reasoning_core/tasks/generated/spatial_folding.py
@@ -0,0 +1,118 @@
+import random
+import re
+from dataclasses import dataclass
+
+from reasoning_core.template import Config, Entry, Task, edict, stochastic_rounding as sround
+from ._base import GeneratedMixin
+
+
+FOLD_TEXT = {
+    "left": "fold the left half over the right half",
+    "right": "fold the right half over the left half",
+    "top": "fold the top half over the bottom half",
+    "bottom": "fold the bottom half over the top half",
+}
+
+
+@dataclass
+class SpatialFoldingConfig(Config):
+    n_folds: int = 1
+    n_punches: int = 1
+
+    def apply_difficulty(self, level):
+        self.n_folds = sround(self.n_folds + 0.6 * level)
+        self.n_punches = sround(self.n_punches + 0.2 * level)
+
+
+def _fold_shape(height, width, direction):
+    if direction in {"left", "right"}:
+        return height, width // 2
+    return height // 2, width
+
+
+def _unfold(points, folds):
+    points = set(points)
+    for direction, height, width in reversed(folds):
+        expanded = set()
+        if direction == "left":
+            half = width // 2
+            for row, col in points:
+                expanded.add((row, half + col))
+                expanded.add((row, half - 1 - col))
+        elif direction == "right":
+            for row, col in points:
+                expanded.add((row, col))
+                expanded.add((row, width - 1 - col))
+        elif direction == "top":
+            half = height // 2
+            for row, col in points:
+                expanded.add((half + row, col))
+                expanded.add((half - 1 - row, col))
+        else:
+            for row, col in points:
+                expanded.add((row, col))
+                expanded.add((height - 1 - row, col))
+        points = expanded
+    return points
+
+
+def _format_points(points):
+    return "; ".join(f"{row + 1},{col + 1}" for row, col in sorted(points))
+
+
+class SpatialFolding(GeneratedMixin, Task):
+    summary = "Track hole positions while folding and unfolding a square grid."
+    config_cls = SpatialFoldingConfig
+
+    def generate_entry(self):
+        n_folds = min(4, max(1, int(self.config.n_folds)))
+        n_punches = min(2, max(1, int(self.config.n_punches)))
+        side = 2 ** max(2, n_folds)
+        height = width = side
+        folds = []
+        directions = []
+
+        for _ in range(n_folds):
+            candidates = []
+            if width % 2 == 0 and width > 1:
+                candidates += ["left", "right"]
+            if height % 2 == 0 and height > 1:
+                candidates += ["top", "bottom"]
+            direction = random.choice(candidates)
+            folds.append((direction, height, width))
+            directions.append(direction)
+            height, width = _fold_shape(height, width, direction)
+
+        positions = [(r, c) for r in range(height) for c in range(width)]
+        punches = random.sample(positions, min(n_punches, len(positions)))
+        holes = _unfold(punches, folds)
+        metadata = edict(
+            side=side,
+            folds=directions,
+            folded_shape=(height, width),
+            punches=punches,
+            holes=sorted(holes),
+        )
+        return Entry(metadata=metadata, answer=_format_points(holes))
+
+    def render_prompt(self, m):
+        fold_text = "\n".join(f"{i}. {FOLD_TEXT[d]}." for i, d in enumerate(m.folds, 1))
+        punches = "; ".join(f"{r + 1},{c + 1}" for r, c in m.punches)
+        return (
+            f"A {m.side}x{m.side} square sheet is divided into unit cells. Rows are numbered top-to-bottom "
+            "and columns left-to-right, starting at 1. After every fold, renumber the visible folded rectangle "
+            "from its new top-left corner.\n"
+            f"Folds, in order:\n{fold_text}\n"
+            f"After all folds the sheet is {m.folded_shape[0]}x{m.folded_shape[1]}. Punch holes through cells: {punches}.\n"
+            "Unfold the sheet completely. The answer is all punched cells as row,column pairs separated by semicolons, "
+            "in row-major order."
+        )
+
+    def score_answer(self, answer, entry):
+        pairs = re.findall(r"(\d+)\s*,\s*(\d+)", str(answer))
+        parsed = [(int(r), int(c)) for r, c in pairs]
+        gold = [(r + 1, c + 1) for r, c in entry.metadata["holes"]]
+        return float(len(parsed) == len(gold) and set(parsed) == set(gold))
+
+    def balancing_key(self, problem):
+        return len(problem.metadata.holes)