用 WorkBuddy 做高考志愿数据分析：从 560 分到院校匹配方案

原始信息（选科/分数/城市/专业偏好）
        ↓
① 结构化学情画像  → 约束条件表
        ↓
② 院校-专业数据整理 → 候选池
        ↓
③ 录取概率量化  → 冲/稳/保 三档
        ↓
④ 生成对比报告 + 可视化看板

帮我整理一份高考志愿分析的学情画像，已知：
- 省份：湖北，科类：理科（选科物化政）
- 预估分数：560（2026 届）
- 目标专业：计算机类（含软件工程、AI、网络安全方向）
- 优先城市：武汉、北京
- 诉求：按"冲稳保"输出院校梯度，并标注选科是否符合
请以表格输出约束条件清单。

school,city,major,score_2023,score_2024,score_2025,require_chem,require_geo
武汉理工大学,武汉,计算机类,554,558,561,是,否
武汉科技大学,武汉,软件工程,538,542,547,是,否
北京信息科技大学,北京,计算机科学与技术,549,553,556,是,否
北京石油化工学院,北京,计算机类,527,531,535,是,否

读取 schools.csv，做三件事：
1. 计算 2023-2025 三年均分与波动；
2. 用 560 分对照，给每所学校标注 冲/稳/保；
3. 检查"物化政"选科是否满足专业要求，不满足的标红。
输出带结论的表格。

import random
def admit_prob(my_score, hist_scores, n=10000):
    cnt = 0
    for _ in range(n):
        sim = my_score + random.uniform(-8, 8)      # 分数波动
        line = random.choice(hist_scores) + random.uniform(-5, 5)
        if sim >= line:
            cnt += 1
    return cnt / n

print("武汉理工 录取概率:", round(admit_prob(560, [554,558,561]), 3))

把上面的分析汇总成一份《2026 志愿填报梯度方案》，
包含：冲（2 所）/ 稳（3 所）/ 保（2 所），
每所注明城市、专业、预估概率、选科符合度，
并用图表展示分数分布。