private void ApplyAlgo(TspAlgo algo)
{
var n = m_Points.Count;
if (n < 3)
{
algo.Path = null;
return;
}
/*var m = new double[n,n];
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
{
m[i, j] = m_Points[i].Dist(m_Points[j]);
}
}
double avgDist = 0;
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
{
avgDist += m[i, j];
}
}
avgDist /= (n * n);
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
{
m[i, j] /= avgDist;
}
}*/
//var measure = new MatrixMeasureFactory().CreateMatrixMeasure(m_Points);
var measure = new PointsMeasure(m_Points);
algo.Path = algo.Solver.GetPath(n, measure);
algo.Length = GetLength(n, m_Points, algo.Path);
}
public static void Main(string[] args)
{
Thread.CurrentThread.CurrentCulture = CultureInfo.InvariantCulture;
//var lines = args[0].Split('\n');
string text;
using (var reader = File.OpenText(args[0]))
{
text = reader.ReadToEnd();
}
var lines = text.Split('\n');
var n = int.Parse(lines[0]);
var points = new Point2DReal[n];
for (int i = 0; i < n; i++)
{
var line = lines[i + 1].Split();
var x = double.Parse(line[0]);
var y = double.Parse(line[1]);
points[i] = new Point2DReal(x, y);
}
/*var matrix = new double[n,n];
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
{
matrix[i, j] = points[i].Dist(points[j]);
}
}
double avgDist = 0;
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
{
avgDist += matrix[i, j];
}
}
avgDist /= (n * n);
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
{
matrix[i, j] /= avgDist;
}
}*/
IMeasure measure;
if (n <= 6000)
{
measure = new MatrixMeasureFactory().CreateMatrixMeasure(points);
}
else
{
measure = new PointsMeasure(points);
}
var path = new Opt3().GetPath(n, measure);
var ans = points[path[0]].Dist(points[path[n - 1]]);
for (int i = 0; i < n - 1; i++)
{
ans += points[path[i]].Dist(points[path[i + 1]]);
}
Console.WriteLine("{0} 0", ans);
Console.WriteLine(string.Join(" ", path));
}
