roi.py 9.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244
  1. from __future__ import annotations
  2. import dataclasses
  3. import json
  4. import typing
  5. import cache
  6. def main() -> None:
  7. recipes: list[Recipe] = cache.get('https://api.prunplanner.org/data/recipes/')
  8. buildings: dict[str, Building] = {m['building_ticker']: m for m in cache.get('https://api.prunplanner.org/data/buildings/')}
  9. materials: dict[str, Material] = {m['ticker']: m for m in cache.get('https://api.prunplanner.org/data/materials/')}
  10. raw_prices: list[RawPrice] = cache.get('https://refined-prun.github.io/refined-prices/all.json')
  11. for cx in ['AI1', 'CI1', 'IC1', 'NC1']:
  12. profits = calc_for_cx(cx, recipes, buildings, materials, raw_prices)
  13. with open(f'www/roi_{cx.lower()}.json', 'w') as f:
  14. json.dump([dataclasses.asdict(p) for p in profits], f, indent='\t')
  15. def calc_for_cx(cx: str, recipes: typing.Collection[Recipe], buildings: typing.Mapping[str, Building],
  16. materials: typing.Mapping[str, Material], raw_prices: typing.Collection[RawPrice]) -> typing.Sequence[Profit]:
  17. prices: dict[str, Price] = {
  18. p['MaterialTicker']: Price(p['VWAP7D'], p['AverageTraded7D'], p['VWAP30D'], p['Bid'], p['Ask']) for p in raw_prices # pyright: ignore[reportArgumentType]
  19. if p['ExchangeCode'] == cx
  20. }
  21. habitation: typing.Mapping[Worker, str] = {
  22. 'pioneers': 'HB1',
  23. 'settlers': 'HB2',
  24. 'technicians': 'HB3',
  25. 'engineers': 'HB4',
  26. 'scientists': 'HB5',
  27. }
  28. hab_area_cost: dict[Worker, float] = {}
  29. hab_capex: dict[Worker, dict[str, float]] = {}
  30. for worker, hab in habitation.items():
  31. hab_area_cost[worker] = buildings[hab]['area_cost'] / 100
  32. base_capex = building_construction_cost(buildings[hab], prices)
  33. hab_capex[worker] = {k: v / 100 for k, v in base_capex.items()}
  34. profits: list[Profit] = []
  35. for recipe in recipes:
  36. if profit := calc_profit(recipe, buildings, hab_area_cost, hab_capex, materials, prices):
  37. profits.append(profit)
  38. profits.sort()
  39. return profits
  40. def get_metrics(amount: float, price: Price) -> dict[str, float]:
  41. v = price.vwap_7d or price.vwap_30d or 0.0
  42. b = price.bid if price.bid is not None else v
  43. a = price.ask if price.ask is not None else v
  44. return {'vwap': amount * v, 'bid': amount * b, 'ask': amount * a}
  45. def calc_profit(recipe: Recipe, buildings: typing.Mapping[str, Building], hab_area_cost: typing.Mapping[Worker, float],
  46. hab_capex: typing.Mapping[Worker, dict[str, float]], materials: typing.Mapping[str, Material],
  47. prices: typing.Mapping[str, Price]) -> Profit | None:
  48. if len(recipe['outputs']) == 0:
  49. return
  50. building = buildings[recipe['building_ticker']]
  51. area = building['area_cost'] + sum(hab_area_cost[worker] * building[worker] for worker in hab_area_cost)
  52. runs_per_day = 24 * 60 * 60 * 1000 / recipe['time_ms'] * 1.25 # assume CoGC
  53. if building['building_ticker'] in ('FRM', 'ORC'):
  54. runs_per_day *= 1.1212 # promitor's fertility
  55. outputs: list[MatPrice] = []
  56. revenue = {'vwap': 0.0, 'bid': 0.0, 'ask': 0.0}
  57. output_prices: dict[str, PriceNonNull] = {}
  58. for output in recipe['outputs']:
  59. price = prices[output['material_ticker']]
  60. if price.vwap_7d is None or price.average_traded_7d is None:
  61. return # skip recipes with thinly traded outputs
  62. output_prices[output['material_ticker']] = typing.cast(PriceNonNull, price)
  63. m = get_metrics(output['material_amount'] * runs_per_day, price)
  64. for k in revenue: revenue[k] += m[k]
  65. outputs.append(MatPrice(output['material_ticker'], output['material_amount'], price.vwap_7d, price.bid, price.ask))
  66. input_costs: list[MatPrice] = []
  67. opex = {'vwap': 0.0, 'bid': 0.0, 'ask': 0.0}
  68. for input in recipe['inputs']:
  69. price = prices[input['material_ticker']]
  70. if price.vwap_7d is None:
  71. return # skip recipes with thinly traded inputs
  72. m = get_metrics(input['material_amount'] * runs_per_day, price)
  73. for k in opex: opex[k] += m[k]
  74. input_costs.append(MatPrice(input['material_ticker'], input['material_amount'], price.vwap_7d, price.bid, price.ask))
  75. worker_consumable = building_daily_cost(building, prices)
  76. for k in opex: opex[k] += worker_consumable[k]
  77. capex = building_construction_cost(building, prices)
  78. for worker, hab_cost in hab_capex.items():
  79. workers = building[worker]
  80. if workers > 0:
  81. for k in capex: capex[k] += hab_cost[k] * workers
  82. lowest_liquidity = min(recipe['outputs'],
  83. key=lambda output: output['material_amount'] / output_prices[output['material_ticker']].average_traded_7d)
  84. output_per_day = lowest_liquidity['material_amount'] * runs_per_day
  85. average_traded_7d = output_prices[lowest_liquidity['material_ticker']].average_traded_7d
  86. # EXTREME DETAIL: We convert our metrics from 'per area' to 'per base' (1 base = 500 area).
  87. # We divide the building's area footprint by 500 to determine what fraction of a base it consumes.
  88. output_per_base = output_per_day / (area / 500)
  89. market_capacity_base = average_traded_7d / output_per_base
  90. logistics_per_base = max(
  91. sum(materials[input['material_ticker']]['weight'] * input['material_amount'] for input in recipe['inputs']),
  92. sum(materials[input['material_ticker']]['volume'] * input['material_amount'] for input in recipe['inputs']),
  93. sum(materials[output['material_ticker']]['weight'] * output['material_amount'] for output in recipe['outputs']),
  94. sum(materials[output['material_ticker']]['volume'] * output['material_amount'] for output in recipe['outputs']),
  95. ) * runs_per_day / (area / 500)
  96. return Profit(outputs, recipe['recipe_name'],
  97. expertise=building['expertise'],
  98. building=building['building_ticker'],
  99. area=area,
  100. capex=capex,
  101. opex=opex,
  102. revenue=revenue,
  103. input_costs=input_costs,
  104. runs_per_day=runs_per_day,
  105. logistics_per_base=logistics_per_base,
  106. output_per_day=output_per_day,
  107. average_traded_7d=average_traded_7d,
  108. market_capacity_base=market_capacity_base)
  109. def building_construction_cost(building: Building, prices: typing.Mapping[str, Price]) -> dict[str, float]:
  110. cost = {'vwap': 0.0, 'bid': 0.0, 'ask': 0.0}
  111. for bc in building['costs']:
  112. m = get_metrics(bc['material_amount'], prices[bc['material_ticker']])
  113. for k in cost: cost[k] += m[k]
  114. # https://handbook.apex.prosperousuniverse.com/wiki/building-costs/#rocky-planets
  115. mcg = get_metrics(building['area_cost'] * 4, prices['MCG'])
  116. for k in cost: cost[k] += mcg[k]
  117. return cost
  118. def building_daily_cost(building: Building, prices: typing.Mapping[str, Price]) -> dict[str, float]:
  119. consumption = {
  120. 'pioneers': [('COF', 0.5), ('DW', 4), ('RAT', 4), ('OVE', 0.5), ('PWO', 0.2)],
  121. 'settlers': [('DW', 5), ('RAT', 6), ('KOM', 1), ('EXO', 0.5), ('REP', 0.2), ('PT', 0.5)],
  122. 'technicians': [('DW', 7.5), ('RAT', 7), ('ALE', 1), ('MED', 0.5), ('SC', 0.1), ('HMS', 0.5), ('SCN', 0.1)],
  123. 'engineers': [('DW', 10), ('MED', 0.5), ('GIN', 1), ('FIM', 7), ('VG', 0.2), ('HSS', 0.2), ('PDA', 0.1)],
  124. 'scientists': [('DW', 10), ('MED', 0.5), ('WIN', 1), ('MEA', 7), ('NST', 0.1), ('LC', 0.2), ('WS', 0.05)],
  125. }
  126. cost = {'vwap': 0.0, 'bid': 0.0, 'ask': 0.0}
  127. for worker, mats in consumption.items():
  128. workers = building[worker]
  129. for mat, per_100 in mats:
  130. m = get_metrics(workers * per_100 / 100, prices[mat])
  131. for k in cost: cost[k] += m[k]
  132. return cost
  133. Worker = typing.Literal['pioneers', 'settlers', 'technicians', 'engineers', 'scientists']
  134. class Recipe(typing.TypedDict):
  135. recipe_name: str
  136. building_ticker: str
  137. inputs: list[RecipeMat]
  138. outputs: list[RecipeMat]
  139. time_ms: int
  140. class RecipeMat(typing.TypedDict):
  141. material_ticker: str
  142. material_amount: int
  143. class Building(typing.TypedDict):
  144. building_ticker: str
  145. expertise: str
  146. area_cost: int
  147. costs: list[BuildingMat]
  148. pioneers: int
  149. settlers: int
  150. technicians: int
  151. engineers: int
  152. scientists: int
  153. class BuildingMat(typing.TypedDict):
  154. material_ticker: str
  155. material_amount: int
  156. class Material(typing.TypedDict):
  157. ticker: str
  158. weight: float
  159. volume: float
  160. class RawPrice(typing.TypedDict):
  161. MaterialTicker: str
  162. ExchangeCode: str
  163. VWAP7D: float | None
  164. AverageTraded7D: float | None
  165. VWAP30D: float | None
  166. Bid: float | None
  167. Ask: float | None
  168. @dataclasses.dataclass(eq=False, frozen=True, slots=True)
  169. class Price:
  170. vwap_7d: float | None
  171. average_traded_7d: float | None
  172. vwap_30d: float | None
  173. bid: float | None
  174. ask: float | None
  175. @dataclasses.dataclass(eq=False, frozen=True, slots=True)
  176. class PriceNonNull:
  177. vwap_7d: float
  178. average_traded_7d: float
  179. @dataclasses.dataclass(eq=False, frozen=True, slots=True)
  180. class Profit:
  181. outputs: typing.Collection[MatPrice]
  182. recipe: str
  183. expertise: str
  184. building: str
  185. area: float
  186. capex: dict[str, float]
  187. opex: dict[str, float]
  188. revenue: dict[str, float]
  189. input_costs: typing.Collection[MatPrice]
  190. runs_per_day: float
  191. logistics_per_base: float # Renamed from area to base
  192. output_per_day: float
  193. average_traded_7d: float
  194. market_capacity_base: float # Renamed from area to base
  195. def __lt__(self, other: Profit) -> bool:
  196. p_a = self.revenue['vwap'] - self.opex['vwap']
  197. be_a = (self.capex['vwap'] + 3 * self.opex['vwap']) / p_a if p_a > 0 else 10000 - p_a
  198. p_b = other.revenue['vwap'] - other.opex['vwap']
  199. be_b = (other.capex['vwap'] + 3 * other.opex['vwap']) / p_b if p_b > 0 else 10000 - p_b
  200. return be_a < be_b
  201. @dataclasses.dataclass(eq=False, frozen=True, slots=True)
  202. class MatPrice:
  203. ticker: str
  204. amount: int
  205. vwap_7d: float
  206. bid: float | None
  207. ask: float | None
  208. if __name__ == '__main__':
  209. main()