# 工厂的工人数数据
factory_sizes <- c(58, 908, 418, 442, 615, 1972, 613, 734, 514, 846)
cumulative_sizes <- cumsum(factory_sizes)  # 累计工人数

# 随机数
set.seed(123)  # 设置随机种子，以确保结果可重复
random_numbers <- sample(1:cumulative_sizes[length(cumulative_sizes)], 2)

# 累积和法：选取的样本
selected_factories_cumulative <- which(random_numbers <= cumulative_sizes)
cat("累积和法选取的工厂编号为：", selected_factories_cumulative, "\n")

# 最大规模法：
# 找到最大的规模，选择对应的工厂
sorted_factories <- order(factory_sizes, decreasing = TRUE)
selected_factories_max_size <- sorted_factories[1:2]
cat("最大规模法选取的工厂编号为：", selected_factories_max_size, "\n")

# 输出随机数，查看选择的工厂编号
cat("生成的随机数为：", random_numbers, "\n")