#shampo.py (shampoopy)
import math
import cmath

# ============================================================
#  Shampoo State Vector
# ============================================================

class ShampooState:
    def __init__(self, Sh=1.0, C=1.0, Ch=1.0, T=1.0, TD=1+0j):
        self.Sh = Sh      # shampoo amplitude
        self.C = C        # camphor oil density
        self.Ch = Ch      # chamomile volatility
        self.T = T        # tulip arbitrage index
        self.TD = TD      # time/date phasor (complex)


# ============================================================
#  ABCD Matrix Application
# ============================================================

def apply_ABCD(matrix, x, y):
    A, B = matrix[0]
    C, D = matrix[1]
    return A * x + B * y, C * x + D * y


# ============================================================
#  Quadrupole Shape Function
# ============================================================

def quad_shape(S):
    return (S.Sh - S.C + S.Ch) * 0.01


# ============================================================
#  Bit Extraction from Phasor
# ============================================================

def bit_from_phasor(TD):
    return 1 if TD.real > 0 else 0


# ============================================================
#  Instruction Execution
# ============================================================

def execute_instruction(instr, S, millionaire_threshold, time):
    itype = instr["type"]

    # -----------------------------
    # Supply-chain instructions
    # -----------------------------
    if itype == "SHIPYARD":
        amt = instr["amount"]
        S.C += amt
        S.Sh += amt * 0.1
        S.T  -= amt * 0.05
        return millionaire_threshold

    if itype == "CHAMOMILE":
        d = instr["delta"]
        S.Ch += d
        S.C  -= d * 0.05
        return millionaire_threshold

    if itype == "WAREHOUSE":
        S.T += instr["tulips"]
        return millionaire_threshold

    # -----------------------------
    # Input instructions
    # -----------------------------
    if itype == "CAMIN":
        v = float(input("Chamomile input: "))
        S.Ch += v
        S.C  -= v * 0.02
        return millionaire_threshold

    if itype == "CAMPIN":
        v = float(input("Camphor oil input: "))
        S.C += v
        S.Sh += v * 0.05
        return millionaire_threshold

    if itype == "TULIPIN":
        v = float(input("Tulip arbitrage input: "))
        S.T += v
        return millionaire_threshold

    # -----------------------------
    # Millionaire threshold
    # -----------------------------
    if itype == "MILLIONAIRE":
        return instr["threshold"]

    # -----------------------------
    # Output instruction
    # -----------------------------
    if itype == "ROTATE":
        ω = instr["omega"]
        S.TD *= cmath.exp(1j * ω * time)
        print(bit_from_phasor(S.TD), end="")
        return millionaire_threshold

    # -----------------------------
    # Quadrupole + ABCD propagation
    # -----------------------------
    if itype == "ATTACK":
        U = instr["U"]
        V = instr["V"]
        ω = instr["omega"]
        ABCD = instr["ABCD"]

        # millionaire-gated collapse
        if S.T <= millionaire_threshold:
            ABCD = [[1, 0], [0, 1]]
            V = 0

        Φ = (U + V * math.cos(ω * time)) * quad_shape(S)

        S.Sh += Φ
        S.C  -= Φ
        S.Ch += Φ * 0.5

        S.Sh, S.C = apply_ABCD(ABCD, S.Sh, S.C)
        S.Ch, S.T = apply_ABCD(ABCD, S.Ch, S.T)

        S.TD *= cmath.exp(1j * ω * time)
        return millionaire_threshold

    raise ValueError(f"Unknown instruction type: {itype}")


# ============================================================
#  Main Interpreter Loop (BATCH / EOF)
# ============================================================

def run_shampoo(program, initial_state=None, millionaire_threshold=5, time_step=1.0):
    S = initial_state or ShampooState()
    time = 0.0
    pc = 0
    batch_stack