import numpy as np
import math
 
print("=" * 60)
print("🚀 Ideone.com 线性回归实验")
print("=" * 60)
 
# 设置随机种子保证结果可重现
np.random.seed(42)
 
# 生成模拟数据：城市人口(万)与利润(万美元)
population = np.random.normal(8, 2.5, 50)
profit = 2.5 * population + 1.0 + np.random.normal(0, 1.5, 50)
 
print("📊 生成数据统计信息:")
print(f"样本数量: {len(population)}")
print(f"人口 - 均值: {np.mean(population):.2f}万, 标准差: {np.std(population):.2f}")
print(f"利润 - 均值: {np.mean(profit):.2f}万美元, 标准差: {np.std(profit):.2f}")
print(f"人口与利润相关系数: {np.corrcoef(population, profit)[0,1]:.3f}")
 
print("\n📈 数据前10行:")
print("人口(万)\t利润(万美元)")
print("-" * 30)
for i in range(min(10, len(population))):
    print(f"{population[i]:6.2f}\t{profit[i]:8.2f}")

aW1wb3J0IG51bXB5IGFzIG5wCmltcG9ydCBtYXRoCgpwcmludCgiPSIgKiA2MCkKcHJpbnQoIvCfmoAgSWRlb25lLmNvbSDnur/mgKflm57lvZLlrp7pqowiKQpwcmludCgiPSIgKiA2MCkKCiMg6K6+572u6ZqP5py656eN5a2Q5L+d6K+B57uT5p6c5Y+v6YeN546wCm5wLnJhbmRvbS5zZWVkKDQyKQoKIyDnlJ/miJDmqKHmi5/mlbDmja7vvJrln47luILkurrlj6Mo5LiHKeS4juWIqea2pijkuIfnvo7lhYMpCnBvcHVsYXRpb24gPSBucC5yYW5kb20ubm9ybWFsKDgsIDIuNSwgNTApCnByb2ZpdCA9IDIuNSAqIHBvcHVsYXRpb24gKyAxLjAgKyBucC5yYW5kb20ubm9ybWFsKDAsIDEuNSwgNTApCgpwcmludCgi8J+TiiDnlJ/miJDmlbDmja7nu5/orqHkv6Hmga86IikKcHJpbnQoZiLmoLfmnKzmlbDph486IHtsZW4ocG9wdWxhdGlvbil9IikKcHJpbnQoZiLkurrlj6MgLSDlnYflgLw6IHtucC5tZWFuKHBvcHVsYXRpb24pOi4yZn3kuIcsIOagh+WHhuW3rjoge25wLnN0ZChwb3B1bGF0aW9uKTouMmZ9IikKcHJpbnQoZiLliKnmtqYgLSDlnYflgLw6IHtucC5tZWFuKHByb2ZpdCk6LjJmfeS4h+e+juWFgywg5qCH5YeG5beuOiB7bnAuc3RkKHByb2ZpdCk6LjJmfSIpCnByaW50KGYi5Lq65Y+j5LiO5Yip5ram55u45YWz57O75pWwOiB7bnAuY29ycmNvZWYocG9wdWxhdGlvbiwgcHJvZml0KVswLDFdOi4zZn0iKQoKcHJpbnQoIlxu8J+TiCDmlbDmja7liY0xMOihjDoiKQpwcmludCgi5Lq65Y+jKOS4hylcdOWIqea2pijkuIfnvo7lhYMpIikKcHJpbnQoIi0iICogMzApCmZvciBpIGluIHJhbmdlKG1pbigxMCwgbGVuKHBvcHVsYXRpb24pKSk6CiAgICBwcmludChmIntwb3B1bGF0aW9uW2ldOjYuMmZ9XHR7cHJvZml0W2ldOjguMmZ9Iik=