2024MathcupA题MR参数的PCI分配基于贪心算法的matlab程序参考

% 假设数据已经导入为表格 'data' data = readtable('data.csv'); % 数据文件名需替换为您的文件路径 % 初始化PCI池和小区编号池 pciPool = 0:1007; cellIdPool = 0:2026; % 计算所有小区的干扰冲突混淆MR之和 totalMR = zeros(height(data), 1); for i = 1:height(data) totalMR(i) = data.ConfusionMR(i) + data.InterferenceConflictMR(i); end data.totalMR = totalMR; % 为每个小区分配PCI，按照MR之和从大到小排序 [~, idx] = sort(totalMR, 'descend'); assignedPCI = zeros(height(data), 1); while ~isempty(cellIdPool) && ~isempty(pciPool) for i = 1:length(idx) currentCell = data.CellID(idx(i)); if ismember(currentCell, cellIdPool) % 从PCI池中分配PCI，并从小区池和PCI池中移除对应编号 assignedPCI(idx(i)) = pciPool(1); pciPool(1) = []; cellIdPool(cellIdPool == currentCell) = []; % 处理邻区 neighbors = data.NeighborCellID{idx(i)}; neighborMRs = zeros(length(neighbors), 1); for j = 1:length(neighbors) neighborIndex = find(data.CellID == neighbors(j)); neighborMRs(j) = data.totalMR(neighborIndex); end [~, neighborOrder] = sort(neighborMRs, 'descend'); for j = 1:length(neighborOrder) neighborCell = neighbors(neighborOrder(j)); if ismember(neighborCell, cellIdPool) && ~isempty(pciPool) assignedPCI(data.CellID == neighborCell) = pciPool(1); pciPool(1) = []; cellIdPool(cellIdPool == neighborCell) = []; end end end % 重新填充PCI池 if isempty(pciPool) pciPool = 0:1007; end end end % 将结果保存回表格 data.AssignedPCI = assignedPCI; writetable(data, 'assigned_pci.csv');