Gruul 15/18/4/24 6 mill Sphinx's tutelage

import random
from collections import Counter

def simulate_mill_realistic_gruul(deck_size=60, red=15, green=18, redgreen=4, lands=23, mill_per_trigger=6, num_simulations=100_000, verbose_simulations=3):
    full_mill_count = 0
    total_milled_list = []

    for sim in range(num_simulations):
        # Deck contains single- and multicolor cards
        deck = ['R'] * red + ['G'] * green + ['RG'] * redgreen + ['L'] * lands
        random.shuffle(deck)

        milled_total = 0
        mill_round = 1
        verbose = sim < verbose_simulations

        if verbose:
            print(f"\n--- Simulation #{sim + 1} ---")

        while True:
            to_mill = min(mill_per_trigger, len(deck))
            if to_mill == 0:
                break

            milled = [deck.pop() for _ in range(to_mill)]
            milled_total += to_mill

            # Tally color presence
            color_counts = {'R': 0, 'G': 0}
            for card in milled:
                if card == 'R':
                    color_counts['R'] += 1
                elif card == 'G':
                    color_counts['G'] += 1
                elif card == 'RG':
                    color_counts['R'] += 1
                    color_counts['G'] += 1
                # Lands are ignored

            has_color_pair = any(count >= 2 for count in color_counts.values())

            if verbose:
                mill_str = ','.join(milled)
                if has_color_pair:
                    print(f"Mill round {mill_round}: Reveals [{mill_str}] → Continue (color pair found)")
                else:
                    print(f"Mill round {mill_round}: Reveals [{mill_str}] → Stop (no color pair)")

            if not has_color_pair:
                break

            mill_round += 1

        total_milled_list.append(milled_total)
        if milled_total >= deck_size:
            full_mill_count += 1

    avg_milled = sum(total_milled_list) / num_simulations
    full_mill_percent = full_mill_count / num_simulations * 100

    print(f"\n=== Summary over {num_simulations} simulations ===")
    print(f"Average cards milled per draw: {avg_milled:.2f}")
    print(f"Full deck mill rate: {full_mill_percent:.2f}%")

simulate_mill_realistic_gruul()