All pastes #443396 Raw Edit

Something

public text v1 · immutable
·published 2007-04-16 19:51 UTC
rendered paste body
     pop   profit      APs
       1   254616      218
      11   300776      218
      21   346936      218
      31   393096      218
      41   439256      436
      51   485416      436
      61   531576      436
      71   577736      436
      81   623896      654
      91   670056      654
     101   716216      654
     111   762376      654
     121   808536      872
     131   854696      872
     141   900856      872
     151   947016      872
     161   993176     1090
     171  1039336     1090
     181  1085496     1090
     191  1131656     1090
     201  1177816     1308
     211  1223976     1308
     221  1270136     1308
     231  1316296     1308
     241  1362456     1526
     251  1408616     1526
     261  1454776     1526
     271  1500936     1744
     281  1547096     1744
     291  1593256     1744
     301  1639416     1744
     311  1685576     1962
     321  1731736     1962
     331  1777896     1962
     341  1824056     1962
     351  1870216     2180
     361  1916376     2180
     371  1962536     2180
     381  2008696     2180
     391  2054856     2398
     401  2101016     2398
     411  2147176     2398
     421  2193336     2398
     431  2239496     2616
     441  2285656     2616
     451  2331816     2616
     461  2377976     2616
     471  2424136     2834
     481  2470296     2834
     491  2516456     2834

---
import math

round_trip_cost = 104*2+10
cargo_capacity = 203
food_price = 34
water_price = 23
energy_price = 100

subsidy = 250000

print '%8s %8s %8s'%('pop', 'profit', 'APs')

for pop in range(1, 500, 10):
  upkeep_cost = pop*(3*food_price+2*water_price-5.25*energy_price)
  profit = pop*200-upkeep_cost+subsidy/8.
  trips = math.ceil(pop*5.25/cargo_capacity)
  ap = round_trip_cost*trips

  print '%8d %8d %8d'%(pop, profit*8, ap)