from pulp import *
import numpy as np
import matplotlib.pyplot as plt

# ユーザー数
n_users = 3

# ユーザーごとの必要なデータレート (Mbps)
data_rates = np.array([5, 10, 8])

# ユーザーごとの信号対雑音比 (SNR) (dB)
snrs = np.array([10, 20, 15])

# 帯域幅 (Hz)
B = 2000000

# ノイズパワー (W)
N = 10**-13

# 問題の定義
prob = LpProblem("Bandwidth_Allocation", LpMinimize)

# 変数の定義
x = [LpVariable(f"x{i}", lowBound=0) for i in range(n_users)]

# 目的関数
prob += lpSum(x)

# 制約条件
for i in range(n_users):
    prob += x[i] * np.log2(1 + snrs[i]) >= data_rates[i]

# 最適化問題を解く
prob.solve()

# 結果を表示
print("最適な帯域幅割り当て: ", [value(xi) for xi in x])

# グラフ化
plt.bar(range(n_users), [value(xi) for xi in x])
plt.xlabel('User')
plt.ylabel('Bandwidth')
plt.show()

import cvxpy as cp
import numpy as np
import matplotlib.pyplot as plt

np.random.seed(1)
n = 50
d = 2
X = np.random.randn(n, d)
delta = np.random.uniform(-0.5, 0.5, size=(n,1))
Y = 2*(X[:,0:1] + delta > X[:,1:2]) - 1

w = cp.Variable((d, 1))
b = cp.Variable()
xi = cp.Variable((n, 1))
C = 1.0

objective = cp.Minimize(cp.norm(w, 2) + C*cp.sum(xi))
constraints = [cp.multiply(Y, X @ w + b) >= 1 - xi, xi >= 0]
prob = cp.Problem(objective, constraints)
prob.solve()

plt.figure(figsize=(8,8))
plt.scatter(X[:,0], X[:,1], c=Y.flatten())
slope = -w.value[0] / w.value[1]
intercept = -b.value / w.value[1]
x_vals = np.linspace(-2, 2, 100)
plt.plot(x_vals, slope * x_vals + intercept, 'k-')
plt.fill_between(x_vals, slope * x_vals + intercept - 1 / np.linalg.norm(w.value), slope * x_vals + intercept + 1 / np.linalg.norm(w.value), color='k', alpha=0.1)
plt.xlim([-2, 2])
plt.ylim([-2, 2])
plt.show()

import cvxpy as cp
import numpy as np
import matplotlib.pyplot as plt

# 資産のリターン
returns = np.array([0.1, 0.2, 0.15])

# 資産のリスク（標準偏差）
risks = np.array([0.05, 0.1, 0.08])

# 資産間の相関係数
correlation = np.array([
    [1, 0.2, 0.5],
    [0.2, 1, 0.3],
    [0.5, 0.3, 1]
])

# 共分散行列を計算
cov_matrix = np.outer(risks, risks) * correlation

# 投資の割合を表す変数
weights = cp.Variable(len(returns))

# リターンを最大化する問題を定義
objective = cp.Maximize(weights @ returns)

# 制約条件を定義
constraints = [cp.sum(weights) == 1,
               weights >= 0,
               cp.quad_form(weights, cov_matrix) <= 0.05**2]

# 問題を定義
problem = cp.Problem(objective, constraints)

# 問題を解く
problem.solve()

# 最適な投資の割合を表示
print(weights.value)

plt.bar(range(len(weights.value)), weights.value)
plt.xlabel('Assets')
plt.ylabel('Investment proportion')
plt.show()

import cvxpy as cp
import numpy as np
import matplotlib.pyplot as plt

# 各食品のコスト（円）
costs = np.array([100, 50, 80])

# 各食品のビタミンCとカリウムの含有量
nutrients = np.array([[0.5, 0.2], [0.1, 0.3], [0.2, 0.2]])

# 必要なビタミンCとカリウムの量
required_nutrients = np.array([1.0, 0.5])

# 食品の購入量を表す変数
x = cp.Variable(3, nonneg=True)

# コスト最小化問題を定義
objective = cp.Minimize(costs @ x)
constraints = [nutrients.T @ x >= required_nutrients, x >= 0]
problem = cp.Problem(objective, constraints)

problem.solve()

print(f"Optimal purchase quantity for Apple: {x.value[0]}")
print(f"Optimal purchase quantity for Banana: {x.value[1]}")
print(f"Optimal purchase quantity for Orange: {x.value[2]}")

plt.bar(['Apple', 'Banana', 'Orange'], x.value)
plt.xlabel('Food')
plt.ylabel('Optimal purchase quantity')
plt.show()

from pulp import *
import matplotlib.pyplot as plt
import numpy as np

# 問題の定義
prob = LpProblem("Discount Optimization", LpMaximize)

# 変数の定義
discounts = ['product_1', 'product_2', 'product_3']
discount_vars = LpVariable.dicts("Discount", discounts, 0, 1)

# 利益の定義
profits = {'product_1': 100, 'product_2': 200, 'product_3': 300}

# 目的関数の定義
prob += lpSum([profits[i]*discount_vars[i] for i in discounts])

# 制約条件の定義
prob += lpSum([discount_vars[i] for i in discounts]) <= 1.5, "Total Discount Constraint"

# 問題の解決
prob.solve()

# 結果の表示
for v in prob.variables():
    print(v.name, "=", v.varValue)

# グラフ化
plt.figure(figsize=(8,6))
plt.bar(range(len(discounts)), [value(var) for var in prob.variables()])

plt.xticks(range(len(discounts)), [var.name for var in prob.variables()])
plt.xlabel('Products')
plt.ylabel('Discount Rate')
plt.title('Optimal Discount Rates for Products')
plt.show()

from pulp import *

# 医師と日付のリストを作成
doctors = ['A', 'B', 'C']
days = ['Mon', 'Tue', 'Wed', 'Thu', 'Fri']

# 問題を定義
prob = LpProblem("Doctors_Scheduling", LpMinimize)

# 変数を定義（各医師が各日に働くかどうか）
x = LpVariable.dicts("x", (doctors, days), cat='Binary')

# 目的関数を定義（医師の働く日数を最小化）
prob += lpSum(x[doctor][day] for doctor in doctors for day in days)

# 制約を定義
for day in days:
    prob += lpSum(x[doctor][day] for doctor in doctors) >= 2  # 1日に最低2人の医師が必要

for doctor in doctors:
    prob += lpSum(x[doctor][day] for day in days) <= 4  # 各医師は1週間に最大4日しか働けない
    prob += lpSum(x[doctor][day] for day in days) >= 2  # 各医師は1週間に最低2日は働く必要があります

# 問題を解く
prob.solve()

# 結果を表示
for v in prob.variables():
    print(v.name, "=", v.varValue)

import matplotlib.pyplot as plt
import numpy as np

# 結果を配列に格納
result = np.zeros((len(doctors), len(days)))
for v in prob.variables():
    doctor, day = v.name.split("_")[1:]
    result[doctors.index(doctor)][days.index(day)] = v.varValue

# ヒートマップを表示
plt.imshow(result, cmap='hot', interpolation='nearest')
plt.colorbar(label='Shifts')
plt.xticks(np.arange(len(days)), days)
plt.yticks(np.arange(len(doctors)), doctors)
plt.show()

from pulp import LpProblem, LpMinimize, LpVariable

# 問題の定義
prob = LpProblem("MedicineOptimization", LpMinimize)

# 変数の定義
A = LpVariable("Medicine_A", 0, None)  # 薬剤Aの単位数
B = LpVariable("Medicine_B", 0, None)  # 薬剤Bの単位数

# 目的関数の定義 (薬剤AとBのコストの最小化)
prob += 5 * A + 7 * B

# 制約条件の定義
prob += 3 * A + 4 * B >= 12  # 栄養素Xの最低摂取量
prob += 2 * A + 3 * B >= 15  # 栄養素Yの最低摂取量

# 問題の解法
prob.solve()

print("Optimal Units of Medicine A: ", A.varValue)
print("Optimal Units of Medicine B: ", B.varValue)

import matplotlib.pyplot as plt
import numpy as np

# 栄養素XとYの最低摂取量
min_X = 12
min_Y = 15

# 薬剤AとBの単位数の範囲
range_A = np.linspace(0, 10, 100)
range_B = np.linspace(0, 10, 100)

# 栄養素XとYの摂取量
X = 3 * range_A + 4 * range_B
Y = 2 * range_A + 3 * range_B

plt.figure(figsize=(8, 6))
plt.plot(range_A, (min_X - 3 * range_A) / 4, label="Nutrient X")
plt.plot(range_A, (min_Y - 2 * range_A) / 3, label="Nutrient Y")
plt.xlim(0, 10)
plt.ylim(0, 10)
plt.xlabel("Units of Medicine A")
plt.ylabel("Units of Medicine B")
plt.legend()
plt.grid(True)
plt.show()

from pulp import LpMaximize, LpProblem, LpStatus, lpSum, LpVariable

# 問題の定義
model = LpProblem(name="small-shop-problem", sense=LpMaximize)

# 変数の定義
x = {i: LpVariable(name=f"x{i}", lowBound=0, upBound=10, cat="Integer") for i in range(1, 4)}

# 利益の定義
profit = {1: 100, 2: 200, 3: 300}

# 目的関数
model += lpSum(profit[i] * x[i] for i in range(1, 4))

# 制約条件
model += (lpSum(x[i] for i in range(1, 4)) <= 2000, "budget_constraint")

# 問題の解
status = model.solve()

for var in x.values():
    print(f"{var.name}: {var.value()}")

print(f"Optimal profit: {model.objective.value()}")

import matplotlib.pyplot as plt

# 商品名
items = ['Item A', 'Item B', 'Item C']

# 最適な在庫数
optimal_values = [x[i].value() for i in range(1, 4)]

plt.bar(items, optimal_values)
plt.xlabel('Items')
plt.ylabel('Optimal Stock')
plt.title('Optimal Stock Level for Each Item')
plt.show()

from pulp import *
import numpy as np
import matplotlib.pyplot as plt

# 各投資商品のリターンとリスク
returns = {'A': 0.1, 'B': 0.12, 'C': 0.10, 'D': 0.07}
risks = {'A': 0.05, 'B': 0.07, 'C': 0.06, 'D': 0.03}

# 投資商品のリスト
assets = list(returns.keys())

# リスクの上限
risk_limit = 0.05

# 問題の定義
prob = LpProblem("Efficient Frontier", LpMaximize)

# 変数の定義
x = LpVariable.dicts("x", assets, 0)

# 目的関数
prob += lpSum([returns[i]*x[i] for i in assets])

# 制約条件
prob += lpSum([risks[i]*x[i] for i in assets]) <= risk_limit
prob += lpSum([x[i] for i in assets]) == 1

# 問題の解
prob.solve()

# 結果の表示
for i in assets:
    print(f"Investment in asset {i}: {x[i].varValue}")

# リスクとリターンのリスト
risks_list = np.linspace(0.01, 0.1, 100)
returns_list = []

for risk_limit in risks_list:
    prob = LpProblem("Efficient Frontier", LpMaximize)
    x = LpVariable.dicts("x", assets, 0)
    prob += lpSum([returns[i]*x[i] for i in assets])
    prob += lpSum([risks[i]*x[i] for i in assets]) <= risk_limit
    prob += lpSum([x[i] for i in assets]) == 1
    prob.solve()
    returns_list.append(value(prob.objective))

# グラフの描画
plt.plot(risks_list, returns_list)
plt.xlabel('Risk')
plt.ylabel('Return')
plt.title('Efficient Frontier')
plt.show()

from pulp import LpProblem, LpMinimize, LpVariable, lpSum
import pandas as pd
import matplotlib.pyplot as plt

# データの設定
# 供給地点と供給量
supply_points = {'Factory1': 500, 'Factory2': 600}
# 需要地点と需要量
demand_points = {'Shop1': 300, 'Shop2': 400, 'Shop3': 400}
# 供給地点から需要地点への輸送コスト
costs = {
    'Factory1': {'Shop1': 2, 'Shop2': 4, 'Shop3': 5},
    'Factory2': {'Shop1': 3, 'Shop2': 1, 'Shop3': 3}
}

# 問題の定義
prob = LpProblem('Transportation', LpMinimize)

# 変数の定義
vars = {
    i: {
        j: LpVariable(f"x({i},{j})", lowBound=0)
        for j in demand_points
    }
    for i in supply_points
}

# 目的関数の定義
prob += lpSum(costs[i][j] * vars[i][j] for i in supply_points for j in demand_points)

# 制約条件の定義
for i in supply_points:
    prob += lpSum(vars[i][j] for j in demand_points) <= supply_points[i]
for j in demand_points:
    prob += lpSum(vars[i][j] for i in supply_points) >= demand_points[j]

# 問題の解法
prob.solve()

# 結果の表示
for v in prob.variables():
    print(v.name, "=", v.varValue)

# 結果のグラフ化
df = pd.DataFrame([(v.name, v.varValue) for v in prob.variables()], columns=['Variable', 'Value'])
df.plot(kind='bar', x='Variable', y='Value')
plt.show()