কস্ট এস্টিমেটিং টেকনিক ও অ্যাকুরেসি

# ১) পার-ইউনিট এস্টিমেট ও ২) ফ্যাক্টর এস্টিমেট -- দুটো ভিন্ন নির্ভুলতার ধাপ

def per_unit_estimate(rate_per_unit, quantity):
    return rate_per_unit * quantity

def accuracy_range(estimate, low_pct, high_pct):
    low = estimate * (1 + low_pct)
    high = estimate * (1 + high_pct)
    return low, high

# ১) পার-ইউনিট এস্টিমেট -- অর্ডার-অফ-ম্যাগনিটিউড / "ballpark" (প্রজেক্টের একদম শুরুতে)
rate = 180           # $ প্রতি বর্গফুট
area = 12_000        # বর্গফুট
rough_estimate = per_unit_estimate(rate, area)
rough_low, rough_high = accuracy_range(rough_estimate, -0.30, 0.50)   # -30% থেকে +50%

print("== ১) পার-ইউনিট (রাফ) এস্টিমেট ==")
print(f"রেট = ${rate}/sq ft, এরিয়া = {area:,} sq ft")
print(f"এস্টিমেট = ${rough_estimate:,.0f}")
print(f"অ্যাকুরেসি রেঞ্জ (-30% থেকে +50%): ${rough_low:,.0f} -- ${rough_high:,.0f}")

# ২) ফ্যাক্টর (রেশিও) এস্টিমেট -- জানা ইকুইপমেন্ট কস্টের উপর শতকরা ফ্যাক্টর যোগ
equipment_cost = 200_000
installation_factor = 0.40
engineering_factor = 0.15
contingency_factor = 0.10

installation_cost = equipment_cost * installation_factor
engineering_cost = equipment_cost * engineering_factor
subtotal = equipment_cost + installation_cost + engineering_cost
contingency_cost = subtotal * contingency_factor
total_installed_cost = subtotal + contingency_cost

detailed_low, detailed_high = accuracy_range(total_installed_cost, -0.05, 0.15)  # -5% থেকে +15%

print("\n== ২) ফ্যাক্টর এস্টিমেট (ডিজাইনের পরের ধাপে, বেশি নিখুঁত) ==")
print(f"ইকুইপমেন্ট কস্ট              = ${equipment_cost:,.0f}")
print(f"+ ইনস্টলেশন (x{installation_factor}) = ${installation_cost:,.0f}")
print(f"+ ইঞ্জিনিয়ারিং (x{engineering_factor})    = ${engineering_cost:,.0f}")
print(f"= সাবটোটাল                   = ${subtotal:,.0f}")
print(f"+ কন্টিনজেন্সি (x{contingency_factor} of সাবটোটাল) = ${contingency_cost:,.0f}")
print(f"= মোট ইনস্টলড কস্ট            = ${total_installed_cost:,.0f}")
print(f"অ্যাকুরেসি রেঞ্জ (-5% থেকে +15%): ${detailed_low:,.0f} -- ${detailed_high:,.0f}")