基本年金額 = 1.2 × (A値 + B値) × (1 + 0.05 × n/12)

def nps_deferral_breakeven(standard_monthly, deferral_years=5):
    """
    国民年金 繰下げ受給の損益分岐点を計算する。

    standard_monthly: 標準受給月額（万ウォン単位）
    deferral_years: 繰下げ年数（最大5年）
    """
    enhancement_rate = 0.072  # 年7.2%加算
    deferred_monthly = standard_monthly * (1 + enhancement_rate * deferral_years)

    # 繰下げ期間中に受け取れなかった総額
    foregone = standard_monthly * 12 * deferral_years

    # 繰下げによる月額増加分
    monthly_gain = deferred_monthly - standard_monthly

    # 損益分岐点（繰下げ後受給開始からの年数）
    breakeven_years = foregone / monthly_gain / 12

    print(f"標準受給月額: {standard_monthly}万ウォン")
    print(f"繰下げ後受給月額: {deferred_monthly:.1f}万ウォン")
    print(f"繰下げ期間中の機会損失: {foregone}万ウォン")
    print(f"月額増加分: {monthly_gain:.1f}万ウォン")
    print(f"損益分岐点: 受給開始から{breakeven_years:.1f}年後")
    print(f"=> 標準受給年齢から{deferral_years + breakeven_years:.1f}年以上生存すれば有利")

    return breakeven_years

# 例：月額120万ウォン、5年繰下げ
nps_deferral_breakeven(standard_monthly=120, deferral_years=5)

import numpy as np

def dc_simulation(
    current_age,
    retirement_age,
    monthly_contribution,   # 万ウォン単位
    current_balance,        # 万ウォン単位
    annual_return,
    inflation=0.025
):
    """
    DC型退職年金の積立シミュレーション。

    Returns: (名目将来価値, 実質将来価値) — 万ウォン単位
    """
    years = retirement_age - current_age
    monthly_rate = annual_return / 12
    n = years * 12

    # 現在残高の将来価値
    fv_balance = current_balance * (1 + annual_return) ** years

    # 毎月の拠出金の将来価値（年金終価式）
    if monthly_rate > 0:
        fv_contributions = monthly_contribution * (
            ((1 + monthly_rate) ** n - 1) / monthly_rate
        )
    else:
        fv_contributions = monthly_contribution * n

    nominal = fv_balance + fv_contributions
    real = nominal / (1 + inflation) ** years

    print(f"=== DC型退職年金 積立シミュレーション ===")
    print(f"現在{current_age}歳 → 退職{retirement_age}歳 ({years}年間)")
    print(f"毎月拠出: {monthly_contribution}万ウォン")
    print(f"現在残高: {current_balance}万ウォン")
    print(f"想定年利: {annual_return*100:.1f}%")
    print()
    print(f"退職時 積立見込み額（名目）: {nominal:,.0f}万ウォン")
    print(f"退職時 積立見込み額（実質）: {real:,.0f}万ウォン")

    return nominal, real

# 例：35歳、60歳定年、月30万ウォン拠出、現在2,000万ウォン
dc_simulation(
    current_age=35,
    retirement_age=60,
    monthly_contribution=30,
    current_balance=2000,
    annual_return=0.06
)

def irp_tax_credit(annual_contribution, annual_wage):
    """
    IRP + 年金貯蓄の税額控除額を計算する。

    annual_contribution: 年間拠出額（万ウォン、最大900）
    annual_wage: 年収（万ウォン）
    """
    max_base = 900
    eligible = min(annual_contribution, max_base)

    if annual_wage <= 5500:
        rate = 0.165
        bracket = "低所得区分（5,500万ウォン以下）"
    else:
        rate = 0.132
        bracket = "標準区分（5,500万ウォン超）"

    credit = eligible * rate

    print(f"年収: {annual_wage * 10:,}ウォン → {bracket}")
    print(f"拠出額: {annual_contribution * 10:,}ウォン")
    print(f"控除対象額: {eligible * 10:,}ウォン")
    print(f"税額控除率: {rate*100:.1f}%")
    print(f"税額控除額: {credit * 10:,}ウォン")
    print(f"実質利回り換算: {credit/annual_contribution*100:.1f}%")

    return credit

# 例1：年収4,000万ウォン、900万ウォン拠出
print("[ 低所得区分の例 ]")
irp_tax_credit(900, 4000)
print()

# 例2：年収7,000万ウォン、900万ウォン拠出
print("[ 標準区分の例 ]")
irp_tax_credit(900, 7000)

年金貯蓄ファンド  → 年600万ウォンまで（税額控除）
IRP             → 年300万ウォン追加（合計900万ウォン上限）

def isa_transfer_strategy(
    isa_maturity_amount,        # 万ウォン単位
    existing_pension_contrib,   # 万ウォン単位
    annual_wage                 # 万ウォン単位
):
    """
    ISA満期資金の年金口座移換による税額控除の計算。
    """
    # ISA移換追加控除基準（移換額の10%、上限300万ウォン）
    isa_extra_base = min(isa_maturity_amount * 0.10, 300)

    # 通常の年金控除基準（上限900万ウォン）
    pension_base = min(existing_pension_contrib, 900)

    if annual_wage <= 5500:
        rate = 0.165
        bracket = "16.5%（低所得区分）"
    else:
        rate = 0.132
        bracket = "13.2%（標準区分）"

    pension_credit = pension_base * rate
    isa_extra_credit = isa_extra_base * rate
    total_credit = pension_credit + isa_extra_credit

    print(f"=== ISA移換戦略 税額控除シミュレーション ===")
    print(f"ISA満期額: {isa_maturity_amount * 10:,}ウォン")
    print(f"ISA追加控除基準: {isa_extra_base * 10:,}ウォン")
    print(f"通常年金控除基準: {pension_base * 10:,}ウォン")
    print(f"適用税額控除率: {bracket}")
    print()
    print(f"年金貯蓄+IRP 税額控除: {pension_credit * 10:,}ウォン")
    print(f"ISA移換 追加税額控除: {isa_extra_credit * 10:,}ウォン")
    print(f"合計 税額控除額: {total_credit * 10:,}ウォン")

# 例：ISA満期3,000万ウォン、年金拠出900万ウォン、年収5,000万ウォン
isa_transfer_strategy(3000, 900, 5000)

def compounding_demo(
    monthly_investment=50,  # 万ウォン/月
    annual_return=0.07,
    years=30
):
    """
    長期積立による複利効果のデモンストレーション。
    """
    monthly_rate = annual_return / 12
    n = years * 12

    fv = monthly_investment * ((1 + monthly_rate) ** n - 1) / monthly_rate
    total_invested = monthly_investment * n
    gain = fv - total_invested

    print(f"毎月投資額: {monthly_investment * 10:,}ウォン")
    print(f"積立期間: {years}年")
    print(f"想定年利: {annual_return*100:.0f}%")
    print(f"元本合計: {total_invested * 10:,}ウォン")
    print(f"最終残高: {fv * 10:,.0f}ウォン")
    print(f"運用益: {gain * 10:,.0f}ウォン")
    print(f"元本倍率: {fv/total_invested:.1f}倍")

# 25歳開始（35年間）vs 35歳開始（25年間）の比較
print("--- 25歳から始めた場合（35年間）---")
compounding_demo(monthly_investment=50, annual_return=0.07, years=35)
print()
print("--- 35歳から始めた場合（25年間）---")
compounding_demo(monthly_investment=50, annual_return=0.07, years=25)

def pension_withdrawal_tax(
    annual_withdrawal,  # 万ウォン単位
    subscriber_age
):
    """
    韓国 私的年金受給時の税金計算。

    年間1,500万ウォン以下：年齢別優遇分離課税
    年間1,500万ウォン超：総合所得税に合算
    """
    THRESHOLD = 1500  # 1,500万ウォン（万ウォン単位）

    # 年齢別分離課税率
    if subscriber_age < 70:
        sep_rate = 0.055   # 5.5%
    elif subscriber_age < 80:
        sep_rate = 0.044   # 4.4%
    else:
        sep_rate = 0.033   # 3.3%

    if annual_withdrawal <= THRESHOLD:
        tax = annual_withdrawal * sep_rate
        tax_method = f"分離課税 {sep_rate*100:.1f}%"
    else:
        base_tax = THRESHOLD * sep_rate
        # 超過分は総合所得税に合算（簡易推計）
        excess_tax = (annual_withdrawal - THRESHOLD) * 0.165
        tax = base_tax + excess_tax
        tax_method = "総合所得税合算（簡易推計）"

    after_tax = annual_withdrawal - tax
    effective_rate = tax / annual_withdrawal * 100

    print(f"年間受給額: {annual_withdrawal * 10:,}ウォン")
    print(f"受給時年齢: {subscriber_age}歳")
    print(f"課税方式: {tax_method}")
    print(f"税額: {tax * 10:,.0f}ウォン")
    print(f"税引後受給額: {after_tax * 10:,.0f}ウォン")
    print(f"実効税率: {effective_rate:.2f}%")

# 例1：72歳、年1,200万ウォン受給
print("[ 72歳・年1,200万ウォン受給 ]")
pension_withdrawal_tax(annual_withdrawal=1200, subscriber_age=72)
print()

# 例2：65歳、年2,000万ウォン受給
print("[ 65歳・年2,000万ウォン受給 ]")
pension_withdrawal_tax(annual_withdrawal=2000, subscriber_age=65)

def full_retirement_simulation(
    # 年齢設定
    current_age=35,
    retirement_age=60,
    life_expectancy=90,

    # 現在の口座残高（万ウォン）
    nps_monthly_estimate=80,      # 国民年金 受給見込み月額（現在価値）
    dc_balance=1000,              # DC型退職年金 現在残高
    dc_monthly_contrib=25,        # DC型 毎月拠出額
    pension_savings_balance=500,  # 年金貯蓄+IRP 現在残高
    pension_monthly_contrib=50,   # 年金貯蓄+IRP 毎月拠出額

    # 運用・物価上昇率の仮定
    dc_return=0.05,
    pension_return=0.06,
    inflation=0.025,

    # 老後に必要な月額生活費（現在価値、万ウォン）
    monthly_living_cost=250
):
    """
    3階建て年金制度を統合した老後資金シミュレーション。
    """
    acc_years = retirement_age - current_age

    # DC型 退職時積立金（名目）
    mr_dc = dc_return / 12
    n_dc = acc_years * 12
    dc_at_retirement = (
        dc_balance * (1 + dc_return) ** acc_years +
        dc_monthly_contrib * ((1 + mr_dc) ** n_dc - 1) / mr_dc
    )

    # 年金貯蓄+IRP 退職時積立金（名目）
    mr_ps = pension_return / 12
    n_ps = acc_years * 12
    ps_at_retirement = (
        pension_savings_balance * (1 + pension_return) ** acc_years +
        pension_monthly_contrib * ((1 + mr_ps) ** n_ps - 1) / mr_ps
    )

    total_assets = dc_at_retirement + ps_at_retirement

    # 国民年金（名目将来価値）
    nps_future_monthly = nps_monthly_estimate * (1 + inflation) ** acc_years
    nps_annual = nps_future_monthly * 12

    # 老後の必要生活費（名目）
    future_monthly_cost = monthly_living_cost * (1 + inflation) ** acc_years
    future_annual_cost = future_monthly_cost * 12

    # 4%ルールによる取り崩し可能額
    annual_dc_withdrawal = dc_at_retirement * 0.04
    annual_ps_withdrawal = ps_at_retirement * 0.04
    total_annual_income = annual_dc_withdrawal + annual_ps_withdrawal + nps_annual

    income_gap = future_annual_cost - total_annual_income

    print("=" * 55)
    print("       総合 老後資金 シミュレーション結果")
    print("=" * 55)
    print(f"現在{current_age}歳 → 退職{retirement_age}歳 → 想定寿命{life_expectancy}歳")
    print()
    print("[退職時の予想積立金]")
    print(f"  DC型退職年金:      {dc_at_retirement:>12,.0f}万ウォン")
    print(f"  年金貯蓄+IRP:      {ps_at_retirement:>12,.0f}万ウォン")
    print(f"  運用資産 合計:     {total_assets:>12,.0f}万ウォン")
    print()
    print("[退職後の年間収入（4%ルール適用）]")
    print(f"  DC取り崩し:        {annual_dc_withdrawal:>12,.0f}万ウォン/年")
    print(f"  年金貯蓄取り崩し:  {annual_ps_withdrawal:>12,.0f}万ウォン/年")
    print(f"  国民年金:          {nps_annual:>12,.0f}万ウォン/年")
    print(f"  年収合計:          {total_annual_income:>12,.0f}万ウォン/年")
    print()
    print("[退職後の必要生活費（名目）]")
    print(f"  月額（名目）:      {future_monthly_cost:>12,.0f}万ウォン")
    print(f"  年額（名目）:      {future_annual_cost:>12,.0f}万ウォン")
    print()
    if income_gap > 0:
        print(f"[警告] 年間収入の不足額: {income_gap:,.0f}万ウォン")
        print("  => 積立額の増加または生活費の見直しが必要です。")
    else:
        print(f"[良好] 年間収入の余裕額: {-income_gap:,.0f}万ウォン")
        print("  => 現状の計画で老後の目標を達成できる見通しです。")

# 35歳会社員のベースシナリオ実行
full_retirement_simulation()