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)