অসম জীবনকালের অল্টারনেটিভ তুলনা

import math

# P/A ও A/P ফ্যাক্টর -- M3-এ প্রতিষ্ঠিত
def pa_factor(i, n):
    return ((1 + i) ** n - 1) / (i * (1 + i) ** n)

def ap_factor(i, n):
    return (i * (1 + i) ** n) / ((1 + i) ** n - 1)

def present_worth(cost, revenue, i, n):
    return -cost + revenue * pa_factor(i, n)

def annual_worth(cost, revenue, i, n):
    return present_worth(cost, revenue, i, n) * ap_factor(i, n)

i = 0.10
cost_X, revenue_X, life_X = 20000, 9000, 3
cost_Y, revenue_Y, life_Y = 28000, 9500, 4

# --- পদ্ধতি ১: সরাসরি AW, প্রতিটি অল্টারনেটিভের নিজস্ব লাইফে (জীবনকাল-নিরপেক্ষ) ---
AW_X = annual_worth(cost_X, revenue_X, i, life_X)
AW_Y = annual_worth(cost_Y, revenue_Y, i, life_Y)

print("পদ্ধতি ১ -- সরাসরি AW (নিজস্ব লাইফে, রিপিটেবিলিটি ধরে নিয়ে)")
print(f"AW_X (৩ বছরের লাইফে) = ${AW_X:,.2f}")
print(f"AW_Y (৪ বছরের লাইফে) = ${AW_Y:,.2f}")
winner_aw = "X" if AW_X > AW_Y else "Y"
print(f"AW পদ্ধতিতে বিজয়ী: অল্টারনেটিভ {winner_aw}")

# --- পদ্ধতি ২: LCM পর্যন্ত ক্যাশ ফ্লো পুনরাবৃত্তি করে PW তুলনা ---
lcm = life_X * life_Y // math.gcd(life_X, life_Y)

def lcm_present_worth(cost, revenue, life, i, horizon):
    cycles = horizon // life
    pw = 0.0
    for c in range(cycles):
        start = c * life
        pw += -cost / (1 + i) ** start          # প্রতি চক্রের শুরুতে প্রাথমিক খরচ (পুনর্বিনিয়োগ)
        for t in range(1, life + 1):
            pw += revenue / (1 + i) ** (start + t)  # প্রতি বছরের বার্ষিক আয়
    return pw

PW_X_lcm = lcm_present_worth(cost_X, revenue_X, life_X, i, lcm)
PW_Y_lcm = lcm_present_worth(cost_Y, revenue_Y, life_Y, i, lcm)

print(f"\nপদ্ধতি ২ -- LCM({life_X}, {life_Y}) = {lcm} বছর পর্যন্ত ক্যাশ ফ্লো রিপিট করে PW")
print(f"অল্টারনেটিভ X {lcm // life_X} বার পুনরাবৃত্ত হয়ে {lcm} বছরে PW_X = ${PW_X_lcm:,.2f}")
print(f"অল্টারনেটিভ Y {lcm // life_Y} বার পুনরাবৃত্ত হয়ে {lcm} বছরে PW_Y = ${PW_Y_lcm:,.2f}")
winner_lcm = "X" if PW_X_lcm > PW_Y_lcm else "Y"
print(f"LCM-PW পদ্ধতিতে বিজয়ী: অল্টারনেটিভ {winner_lcm}")

print(f"\nদুই পদ্ধতির র‍্যাঙ্কিং মিলেছে? {winner_aw == winner_lcm}")

# ক্রস-চেক: LCM-PW সবসময় AW x P/A(i, lcm)-এর সমান হওয়া উচিত (যেহেতু AW লাইফ-নিরপেক্ষ সমতুল্য মূল্য)
pa_lcm = pa_factor(i, lcm)
print(f"\nক্রস-চেক: AW_X x P/A(i,{lcm}) = {AW_X:,.2f} x {pa_lcm:.6f} = ${AW_X * pa_lcm:,.2f} (তুলনা করুন PW_X উপরের সাথে)")
print(f"ক্রস-চেক: AW_Y x P/A(i,{lcm}) = {AW_Y:,.2f} x {pa_lcm:.6f} = ${AW_Y * pa_lcm:,.2f} (তুলনা করুন PW_Y উপরের সাথে)")