public QualityRecoveringOptimizer(Data data, Solution solutionBefore, ProblemFormulation formulation, CCTModel.MIPModelParameters mipModelParameters)
{
this.data = data;
this.formulation = formulation;
this.MIPmodelParameters = mipModelParameters;
this.solutionBefore = solutionBefore;
model = new CCTModel(data, formulation, MIPmodelParameters);
model.SetProximityOrigin(solutionBefore._assignments.ToList());
model.SetMipHints(solutionBefore._assignments.ToList());
}
public void TopSolutionsAnalysis(string filename)
{
var name = Path.GetFileNameWithoutExtension(filename);
var data = Data.ReadXml(dataPath + filename, "120", "4");
var formulation = ProblemFormulation.UD2NoOverbook;
var courestime = new Dictionary <Course, Dictionary <TimeSlot, int> >();
var penaltyforcourse = new Dictionary <Course, int>();
foreach (var course in data.Courses)
{
courestime.Add(course, new Dictionary <TimeSlot, int>());
foreach (var timeSlot in data.TimeSlots)
{
courestime[course].Add(timeSlot, 0);
}
penaltyforcourse[course] = 0;
}
var bestsol = new List <Assignment>();
var penMWD = new List <int>();
var penCC = new List <int>();
for (var i = 0; i < 10; i++)
{
var sol = new Solution(data, formulation);
sol.Read(dataPath + @"ITC2007\sol\" + name + "-UD2-" + i + ".sol");
sol.AnalyzeSolution();
foreach (var assignment in sol._assignments)
{
courestime[assignment.Course][assignment.TimeSlot]++;
}
foreach (var course in data.Courses)
{
penaltyforcourse[course] += sol.PenaltyForCourse[course];
}
penMWD.Add((int)sol.MinimumWorkingDays);
penCC.Add((int)sol.CurriculumCompactness);
if (i == 0)
{
bestsol = sol._assignments.ToList();
}
}
const int minhint = 7;
var maxhint = 7;
var noise = 0.5;
const bool fixhints = true;
var random = new Random(62);
foreach (var course in data.Courses)
{
foreach (var ts in data.TimeSlots)
{
// if (random.Next(10) < courestime[course][ts]) courestime[course][ts] = 1;
// else courestime[course][ts] = 0;
if (courestime[course][ts] < minhint)
{
courestime[course][ts] = 0;
}
if (courestime[course][ts] > maxhint)
{
courestime[course][ts] = 0;
}
//result should be made to double.
// if(courestime[course][ts] > 0) courestime[course][ts] += random.Next(-3,3);
//courestime[course][ts] *= 1 + (int) Math.Round(noise*(2*random.NextDouble()-1));
}
}
var hints = courestime.SelectMany(c => c.Value.Where(cv => cv.Value > 0).Select(cv => Tuple.Create(c.Key, cv.Key, cv.Value))).ToList();
var currtime = new Dictionary <Curriculum, Dictionary <TimeSlot, int> >();
foreach (var curr in data.Curricula)
{
currtime.Add(curr, new Dictionary <TimeSlot, int>());
foreach (var timeSlot in data.TimeSlots)
{
currtime[curr].Add(timeSlot, courestime.Where(c => curr.Courses.Contains(c.Key)).Sum(kv => kv.Value[timeSlot]));
}
}
var chints = currtime.SelectMany(c => c.Value.Select(cv => Tuple.Create(c.Key, cv.Key, cv.Value))).ToList();
Console.WriteLine($"Timeslots: {data.TimeSlots.Count}");
Console.WriteLine($"MWD: avg: {penMWD.Average()} - {string.Join(", ", penMWD)} ");
Console.WriteLine($"CC: avg: {penCC.Average()} - {string.Join(", ", penCC)} ");
foreach (var course in data.Courses.OrderByDescending(c => penaltyforcourse[c]))
{
Console.WriteLine($"\nCourse: {course}\n" +
$"Size: {course.NumberOfStudents}\n" +
$"Lectures: {course.Lectures}\n" +
$"MWD: {course.MinimumWorkingDays}\n" +
$"TSAvaialble:{data.TimeSlots.Count - course.UnavailableTimeSlots.Count}\n" +
$"Unavail: {string.Join(",", course.UnavailableTimeSlots)}\n" +
$"TotalPenalty: {penaltyforcourse[course]}");
foreach (var valuePair in courestime[course].Where(kv => kv.Value > 0).OrderByDescending(kv => kv.Value))
{
Console.WriteLine(valuePair.Key + ": " + valuePair.Value);
}
}
Console.WriteLine("\n\nBy prio: ");
foreach (var hint in hints.Where(h => h.Item3 > 0).OrderByDescending(h => h.Item3))
{
Console.WriteLine(hint.Item1.ID + "_" + hint.Item2 + ": " + hint.Item3);
}
Console.WriteLine("\n\nHistogram: ");
Console.WriteLine("total lectures: " + data.Courses.Sum(c => c.Lectures));
for (var i = 1; i <= 10; i++)
{
Console.WriteLine($"{i}\t{hints.Count(h => h.Item3 == i)}");
}
Console.WriteLine("Unavailabilites\tSame");
foreach (var course in data.Courses)
{
Console.WriteLine($"{course.UnavailableTimeSlots.Count}\t{(double) courestime[course].Sum(t => t.Value)/course.Lectures : 0.00}");
}
Console.WriteLine("Penal\tSame");
foreach (var course in data.Courses)
{
Console.WriteLine($"{penaltyforcourse[course]}\t{(double) courestime[course].Sum(t => t.Value)/course.Lectures : 0.00}");
}
Console.WriteLine("Size\tSame");
foreach (var course in data.Courses)
{
Console.WriteLine($"{course.NumberOfStudents}\t{(double) courestime[course].Sum(t => t.Value)/course.Lectures : 0.00}");
}
// return;
//Console.WriteLine("\n\nCurriculumPatterns");
//foreach (var hint in chints.Where(h => h.Item3 > 0).OrderByDescending(h => h.Item3))
//{
// Console.WriteLine(hint.Item1.ID + "_" + hint.Item2 + ": " + hint.Item3);
//}
//Console.WriteLine("\n\nHistogram: ");
//Console.WriteLine("total curr_lectures: " + data.Curricula.Sum(c => c.Lectures));
//for (var i = 1; i <= 10; i++)
//{
// Console.WriteLine($"{i}\t{chints.Count(h => h.Item3 == i)}");
//}
var ModelParameters = new CCTModel.MIPModelParameters()
{
UseStageIandII = false,
TuneGurobi = true,
UseHallsConditions = true,
UseRoomHallConditions = true,
// Seed = seeds[i],
};
var model = new CCTModel(data, formulation, ModelParameters);
// model.SetMipHints(chints);
Console.WriteLine($"Adding {hints.Count} hints. Total letures: {data.Courses.Sum(c => c.Lectures)} ({(double) hints.Count / data.Courses.Sum(c => c.Lectures) :0.00%})");
model.SetMipHints(hints, fixhints);
// model.MipStart(bestsol);
var r = model.Optimize(300);
// model.Optimize(60*60*3);
var ass = model.GetAssignments();
var finalsol = new Solution(data, formulation);
finalsol.SetAssignments(ass);
finalsol.Save(@"C:\temp\finalsolStageI_" + name + ".sol");
var ModelParameters2 = new CCTModel.MIPModelParameters()
{
UseStageIandII = true,
TuneGurobi = true,
UseHallsConditions = true,
UseRoomHallConditions = true,
// Seed = seeds[i],
};
var model2 = new CCTModel(data, formulation, ModelParameters2);
model2.SetAndFixSolution(ass);
model2.Optimize(300);
ass = model2.GetAssignments();
finalsol.SetAssignments(ass);
finalsol.Save(@"C:\temp\finalsolStageII_" + name + ".sol");
}
