Mesa 3.0.0a5__py3-none-any.whl → 3.0.0b1__py3-none-any.whl

examples/advanced/sugarscape_g1mt/Readme.md

@@ -0,0 +1,87 @@
+# Sugarscape Constant Growback Model with Traders
+
+## Summary
+
+This is Epstein & Axtell's Sugarscape model with Traders, a detailed description is in Chapter four of
+*Growing Artificial Societies: Social Science from the Bottom Up (1996)*. The model shows an emergent price equilibrium can happen via a decentralized dynamics.
+
+This code generally matches the code in the Complexity Explorer Tutorial, but in `.py` instead of `.ipynb` format.
+
+### Agents:
+
+- **Resource**:  Resource agents grow back at one unit of sugar and spice per time step up to a specified max amount and can be harvested and traded by the trader agents.
+  (if you do the interactive run, the color will be green if the resource agent has a bigger amount of sugar, or yellow if it has a bigger amount of spice)
+- **Traders**: Trader agents have the following attributes: (1) metabolism for sugar, (2) metabolism for spice, (3) vision,
+  (4) initial sugar endowment and (5) initial spice endowment. The traverse the landscape harvesting sugar and spice and
+trading with other agents. If they run out of sugar or spice then they are removed from the model. (red circle if you do the interactive run)
+
+The trader agents traverse the landscape according to rule **M**:
+- Look out as far as vision permits in the four principal lattice directions and identify the unoccupied site(s).
+- Considering only unoccupied sites find the nearest position that produces the most welfare using the Cobb-Douglas function.
+- Move to the new position
+- Collect all the resources (sugar and spice) at that location
+(Epstein and Axtell, 1996, p. 99)
+
+The traders trade according to rule **T**:
+- Agents and potential trade partner compute their marginal rates of substitution (MRS), if they are equal *end*.
+- Exchange resources, with spice flowing from the agent with the higher MRS to the agent with the lower MRS and sugar
+flowing the opposite direction.
+- The price (p) is calculated by taking the geometric mean of the agents' MRS.
+- If p > 1 then p units of spice are traded for 1 unit of sugar; if p < 1 then 1/p units of sugar for 1 unit of spice
+- The trade occurs if it will (a) make both agent better off (increases MRS) and (b) does not cause the agents' MRS to
+cross over one another otherwise *end*.
+- This process then repeats until an *end* condition is met.
+(Epstein and Axtell, 1996, p. 105)
+
+The model demonstrates several Mesa concepts and features:
+ - MultiGrid
+ - Multiple agent types (traders, sugar, spice)
+ - Dynamically removing agents from the grid and schedule when they die
+ - Data Collection at the model and agent level
+ - Batchrunner (i.e. parameter sweeps)
+
+## Installation
+
+To install the dependencies use pip and the requirements.txt in this directory. e.g.
+
+```
+  $ pip install -r requirements.txt
+```
+
+## How to Run
+
+To run the model a single instance of the model:
+
+```
+  $ python run.py -s
+```
+
+To run the model with BatchRunner:
+
+```
+  $ python run.py -b
+```
+
+To run the model interactively:
+
+```
+  $ mesa runserver
+```
+
+Then open your browser to [http://127.0.0.1:8521/](http://127.0.0.1:8521/) and press Reset, then Run.
+
+## Files
+
+* `sugarscape_g1mt/trader_agents.py`: Defines the Trader agent class.
+* `sugarscape_g1mt/resource_agents.py`: Defines the Resource agent class which contains an amount of sugar and spice.
+* `sugarscape_g1mt/model.py`: Manages the Sugarscape Constant Growback with Traders model.
+* `sugarscape_g1mt/sugar_map.txt`: Provides sugar and spice landscape in raster type format.
+* `server.py`: Sets up an interactive visualization server.
+* `run.py`: Runs Server, Single Run or Batch Run  with data collection and basic analysis.
+* `app.py`: Runs a visualization server via Solara (`solara run app.py`).
+* `tests.py`: Has tests to ensure that the model reproduces the results in shown in Growing Artificial Societies.
+
+## Additional Resources
+
+- [Growing Artificial Societies](https://mitpress.mit.edu/9780262550253/growing-artificial-societies/)
+- [Complexity Explorer Sugarscape with Traders Tutorial](https://www.complexityexplorer.org/courses/172-agent-based-models-with-python-an-introduction-to-mesa)

examples/advanced/sugarscape_g1mt/app.py

@@ -0,0 +1,61 @@
+import numpy as np
+import solara
+from matplotlib.figure import Figure
+from mesa.visualization import SolaraViz, make_plot_measure
+from sugarscape_g1mt.model import SugarscapeG1mt
+from sugarscape_g1mt.trader_agents import Trader
+
+
+def SpaceDrawer(model):
+    def portray(g):
+        layers = {
+            "sugar": [[np.nan for j in range(g.height)] for i in range(g.width)],
+            "spice": [[np.nan for j in range(g.height)] for i in range(g.width)],
+            "trader": {"x": [], "y": [], "c": "tab:red", "marker": "o", "s": 10},
+        }
+
+        for content, (i, j) in g.coord_iter():
+            for agent in content:
+                if isinstance(agent, Trader):
+                    layers["trader"]["x"].append(i)
+                    layers["trader"]["y"].append(j)
+                else:
+                    # Don't visualize resource with value <= 1.
+                    layers["sugar"][i][j] = (
+                        agent.sugar_amount if agent.sugar_amount > 1 else np.nan
+                    )
+                    layers["spice"][i][j] = (
+                        agent.spice_amount if agent.spice_amount > 1 else np.nan
+                    )
+        return layers
+
+    fig = Figure()
+    ax = fig.subplots()
+    out = portray(model.grid)
+    # Sugar
+    # Important note: imshow by default draws from upper left. You have to
+    # always explicitly specify origin="lower".
+    im = ax.imshow(out["sugar"], cmap="spring", origin="lower")
+    fig.colorbar(im, orientation="vertical")
+    # Spice
+    ax.imshow(out["spice"], cmap="winter", origin="lower")
+    # Trader
+    ax.scatter(**out["trader"])
+    ax.set_axis_off()
+    return solara.FigureMatplotlib(fig)
+
+
+model_params = {
+    "width": 50,
+    "height": 50,
+}
+
+model1 = SugarscapeG1mt(50, 50)
+
+page = SolaraViz(
+    model1,
+    components=[SpaceDrawer, make_plot_measure(["Trader", "Price"])],
+    name="Sugarscape {G1, M, T}",
+    play_interval=1500,
+)
+page  # noqa

examples/advanced/sugarscape_g1mt/requirements.txt

@@ -0,0 +1,6 @@
+jupyter
+mesa~=2.0
+numpy
+matplotlib
+networkx
+pandas

examples/advanced/sugarscape_g1mt/run.py

@@ -0,0 +1,105 @@
+import sys
+
+import matplotlib.pyplot as plt
+import mesa
+import networkx as nx
+import pandas as pd
+from sugarscape_g1mt.model import SugarscapeG1mt
+from sugarscape_g1mt.server import server
+
+
+# Analysis
+def assess_results(results, single_agent):
+    # Make dataframe of results
+    results_df = pd.DataFrame(results)
+    # Plot and show  mean price
+    plt.scatter(results_df["Step"], results_df["Price"], s=0.75)
+    plt.show()
+
+    if single_agent is not None:
+        plt.plot(results_df["Step"], results_df["Trader"])
+        plt.show()
+    else:
+        n = max(results_df["RunId"])
+        # Plot number of Traders
+        for i in range(n):
+            results_explore = results_df[results_df["RunId"] == i]
+            plt.plot(results_explore["Step"], results_explore["Trader"])
+        plt.show()
+
+    if single_agent is not None:
+        results_df = single_agent
+
+    # Show Trade Networks
+    #  create graph object
+    print("Making Network")
+    G = nx.Graph()
+    trade = results_df.dropna(subset=["Trade Network"])
+    # add agent keys to make initial node set
+    G.add_nodes_from(list(trade["AgentID"].unique()))
+
+    # create edge list
+    for idx, row in trade.iterrows():
+        if len(row["Trade Network"]) > 0:
+            for agent in row["Trade Network"]:
+                G.add_edge(row["AgentID"], agent)
+
+    # Get Basic Network Statistics
+    print(f"Node Connectivity {nx.node_connectivity(G)}")
+    print(f"Average Clustering {nx.average_clustering(G)}")
+    print(f"Global Efficiency {nx.global_efficiency(G)}")
+
+    # Plot histogram of degree distribution
+    degree_sequence = sorted((d for n, d in G.degree()), reverse=True)
+    degree_sequence = [d for n, d in G.degree()]
+    plt.hist(degree_sequence)
+    plt.show()
+
+    # Plot network
+    nx.draw(G)
+    plt.show()
+
+
+# Run the model
+def main():
+    args = sys.argv[1:]
+
+    if len(args) == 0:
+        server.launch()
+
+    elif args[0] == "-s":
+        print("Running Single Model")
+        model = SugarscapeG1mt()
+        model.run_model()
+        model_results = model.datacollector.get_model_vars_dataframe()
+        model_results["Step"] = model_results.index
+        agent_results = model.datacollector.get_agent_vars_dataframe()
+        agent_results = agent_results.reset_index()
+        assess_results(model_results, agent_results)
+
+    elif args[0] == "-b":
+        print("Conducting a Batch Run")
+        params = {
+            "width": 50,
+            "height": 50,
+            "vision_min": range(1, 4),
+            "metabolism_max": [2, 3, 4, 5],
+        }
+
+        results_batch = mesa.batch_run(
+            SugarscapeG1mt,
+            parameters=params,
+            iterations=1,
+            number_processes=1,
+            data_collection_period=1,
+            display_progress=True,
+        )
+
+        assess_results(results_batch, None)
+
+    else:
+        raise Exception("Option not found")
+
+
+if __name__ == "__main__":
+    main()

examples/advanced/sugarscape_g1mt/sugarscape_g1mt/model.py

@@ -0,0 +1,180 @@
+from pathlib import Path
+
+import mesa
+import numpy as np
+from mesa.experimental.cell_space import OrthogonalVonNeumannGrid
+
+from .resource_agents import Resource
+from .trader_agents import Trader
+
+
+# Helper Functions
+def flatten(list_of_lists):
+    """
+    helper function for model datacollector for trade price
+    collapses agent price list into one list
+    """
+    return [item for sublist in list_of_lists for item in sublist]
+
+
+def geometric_mean(list_of_prices):
+    """
+    find the geometric mean of a list of prices
+    """
+    return np.exp(np.log(list_of_prices).mean())
+
+
+def get_trade(agent):
+    """
+    For agent reporters in data collector
+
+    return list of trade partners and None for other agents
+    """
+    if isinstance(agent, Trader):
+        return agent.trade_partners
+    else:
+        return None
+
+
+class SugarscapeG1mt(mesa.Model):
+    """
+    Manager class to run Sugarscape with Traders
+    """
+
+    def __init__(
+        self,
+        width=50,
+        height=50,
+        initial_population=200,
+        endowment_min=25,
+        endowment_max=50,
+        metabolism_min=1,
+        metabolism_max=5,
+        vision_min=1,
+        vision_max=5,
+        enable_trade=True,
+    ):
+        super().__init__()
+        # Initiate width and height of sugarscape
+        self.width = width
+        self.height = height
+        # Initiate population attributes
+        self.initial_population = initial_population
+        self.endowment_min = endowment_min
+        self.endowment_max = endowment_max
+        self.metabolism_min = metabolism_min
+        self.metabolism_max = metabolism_max
+        self.vision_min = vision_min
+        self.vision_max = vision_max
+        self.enable_trade = enable_trade
+        self.running = True
+
+        # initiate mesa grid class
+        self.grid = OrthogonalVonNeumannGrid((self.width, self.height), torus=False)
+        # initiate datacollector
+        self.datacollector = mesa.DataCollector(
+            model_reporters={
+                "Trader": lambda m: len(m.agents_by_type[Trader]),
+                "Trade Volume": lambda m: sum(
+                    len(a.trade_partners) for a in m.agents_by_type[Trader]
+                ),
+                "Price": lambda m: geometric_mean(
+                    flatten([a.prices for a in m.agents_by_type[Trader]])
+                ),
+            },
+            agent_reporters={"Trade Network": lambda a: get_trade(a)},
+        )
+
+        # read in landscape file from supplmentary material
+        sugar_distribution = np.genfromtxt(Path(__file__).parent / "sugar-map.txt")
+        spice_distribution = np.flip(sugar_distribution, 1)
+
+        for cell in self.grid.all_cells:
+            max_sugar = sugar_distribution[cell.coordinate]
+            max_spice = spice_distribution[cell.coordinate]
+            Resource(self, max_sugar, max_spice, cell)
+
+        for _ in range(self.initial_population):
+            # get agent position
+            x = self.random.randrange(self.width)
+            y = self.random.randrange(self.height)
+            # see Growing Artificial Societies p. 108 for initialization
+            # give agents initial endowment
+            sugar = int(self.random.uniform(self.endowment_min, self.endowment_max + 1))
+            spice = int(self.random.uniform(self.endowment_min, self.endowment_max + 1))
+            # give agents initial metabolism
+            metabolism_sugar = int(
+                self.random.uniform(self.metabolism_min, self.metabolism_max + 1)
+            )
+            metabolism_spice = int(
+                self.random.uniform(self.metabolism_min, self.metabolism_max + 1)
+            )
+            # give agents vision
+            vision = int(self.random.uniform(self.vision_min, self.vision_max + 1))
+
+            cell = self.grid[(x, y)]
+            # create Trader object
+            Trader(
+                self,
+                cell,
+                sugar=sugar,
+                spice=spice,
+                metabolism_sugar=metabolism_sugar,
+                metabolism_spice=metabolism_spice,
+                vision=vision,
+            )
+
+    def step(self):
+        """
+        Unique step function that does staged activation of sugar and spice
+        and then randomly activates traders
+        """
+        # step Resource agents
+        self.agents_by_type[Resource].do("step")
+
+        # step trader agents
+        # to account for agent death and removal we need a separate data structure to
+        # iterate
+        trader_shuffle = self.agents_by_type[Trader].shuffle()
+
+        for agent in trader_shuffle:
+            agent.prices = []
+            agent.trade_partners = []
+            agent.move()
+            agent.eat()
+            agent.maybe_die()
+
+        if not self.enable_trade:
+            # If trade is not enabled, return early
+            self.datacollector.collect(self)
+            return
+
+        trader_shuffle = self.agents_by_type[Trader].shuffle()
+
+        for agent in trader_shuffle:
+            agent.trade_with_neighbors()
+
+        # collect model level data
+        self.datacollector.collect(self)
+        """
+        Mesa is working on updating datacollector agent reporter
+        so it can collect information on specific agents from
+        mesa.time.RandomActivationByType.
+
+        Please see issue #1419 at
+        https://github.com/projectmesa/mesa/issues/1419
+        (contributions welcome)
+
+        Below is one way to update agent_records to get specific Trader agent data
+        """
+        # Need to remove excess data
+        # Create local variable to store trade data
+        agent_trades = self.datacollector._agent_records[self.steps]
+        # Get rid of all None to reduce data storage needs
+        agent_trades = [agent for agent in agent_trades if agent[2] is not None]
+        # Reassign the dictionary value with lean trade data
+        self.datacollector._agent_records[self.steps] = agent_trades
+
+    def run_model(self, step_count=1000):
+        for i in range(step_count):
+            self.step()

examples/advanced/sugarscape_g1mt/sugarscape_g1mt/resource_agents.py

@@ -0,0 +1,26 @@
+from mesa.experimental.cell_space import FixedAgent
+
+
+class Resource(FixedAgent):
+    """
+    Resource:
+    - contains an amount of sugar and spice
+    - grows 1 amount of sugar at each turn
+    - grows 1 amount of spice at each turn
+    """
+
+    def __init__(self, model, max_sugar, max_spice, cell):
+        super().__init__(model)
+        self.sugar_amount = max_sugar
+        self.max_sugar = max_sugar
+        self.spice_amount = max_spice
+        self.max_spice = max_spice
+        self.cell = cell
+
+    def step(self):
+        """
+        Growth function, adds one unit of sugar and spice each step up to
+        max amount
+        """
+        self.sugar_amount = min([self.max_sugar, self.sugar_amount + 1])
+        self.spice_amount = min([self.max_spice, self.spice_amount + 1])

examples/advanced/sugarscape_g1mt/sugarscape_g1mt/server.py

@@ -0,0 +1,61 @@
+import mesa
+
+from .model import SugarscapeG1mt
+from .resource_agents import Resource
+from .trader_agents import Trader
+
+sugar_dic = {4: "#005C00", 3: "#008300", 2: "#00AA00", 1: "#00F800"}
+spice_dic = {4: "#acac00", 3: "#c5c500", 2: "#dfdf00", 1: "#f8f800"}
+
+
+def Agent_portrayal(agent):
+    if agent is None:
+        return
+
+    if isinstance(agent, Trader):
+        return {
+            "Shape": "circle",
+            "Filled": "true",
+            "r": 0.5,
+            "Layer": 0,
+            "Color": "#FF0A01",
+        }
+
+    elif isinstance(agent, Resource):
+        resource_type = "sugar" if agent.max_sugar > agent.max_spice else "spice"
+        if resource_type == "sugar":
+            color = (
+                sugar_dic[agent.sugar_amount] if agent.sugar_amount != 0 else "#D6F5D6"
+            )
+            layer = 1 if agent.sugar_amount > 2 else 0
+        else:
+            color = (
+                spice_dic[agent.spice_amount] if agent.spice_amount != 0 else "#D6F5D6"
+            )
+            layer = 1 if agent.spice_amount > 2 else 0
+        return {
+            "Color": color,
+            "Shape": "rect",
+            "Filled": "true",
+            "Layer": layer,
+            "w": 1,
+            "h": 1,
+        }
+
+    return {}
+
+
+canvas_element = mesa.visualization.CanvasGrid(Agent_portrayal, 50, 50, 500, 500)
+chart_element = mesa.visualization.ChartModule(
+    [{"Label": "Trader", "Color": "#AA0000"}]
+)
+chart_element2 = mesa.visualization.ChartModule(
+    [{"Label": "Price", "Color": "#000000"}]
+)
+
+server = mesa.visualization.ModularServer(
+    SugarscapeG1mt,
+    [canvas_element, chart_element, chart_element2],
+    "Sugarscape with Traders",
+)
+# server.launch()

examples/advanced/sugarscape_g1mt/sugarscape_g1mt/sugar-map.txt

@@ -0,0 +1,50 @@
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 2 2 2 2 2 2 2 2
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 3 4 4 4 4 3 3 3 3 3 3 3 2 2
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 4 4 3 3 3 3 3 3 2
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 2
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 3
+0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3
+0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3
+0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3
+0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3
+0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 3
+0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 2
+0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 3 3 3 3 3 3 2
+1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 4 4 4 4 3 3 3 3 3 3 3 2 2
+1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
+1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2
+1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2
+1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2
+1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2
+1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1
+1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1
+1 1 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1
+1 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1
+1 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1
+2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1
+2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1
+2 2 2 2 2 2 3 3 3 3 3 3 3 4 4 4 4 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 4 4 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 0 0 0
+2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 3 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0
+2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 3 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0
+2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0
+2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0
+2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0
+2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 2 2 2 2 2 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0
+2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 3 2 2 2 2 2 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0
+2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 3 3 3 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0
+2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 4 4 3 3 3 3 3 3 2 2 2 2 2 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0
+2 2 2 2 2 2 3 3 3 3 3 3 3 4 4 4 4 3 3 3 3 3 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0
+2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0
+2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0
+1 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0
+1 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0
+1 1 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0
+1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0