Карта связей

class Node
{
    public Point Position;
    public Point Disp;
}

public class Edge
{
    public Node From;
    public Node To;
}

public class App
{
    double fr(double z, double k) { return k * k / z; }
    double fa(double z, double k) { return z * z / k; }

    void Spread()
    {
        int maxwidth = m_Nodes.Count * 30;
        int maxheight = maxwidth;

        double temperature = maxwidth / 10.0;
        double t = temperature;

        double area = (double)maxwidth * maxheight;
        double k = Math.Sqrt(area / m_Nodes.Count);

        int itercount = 30;

        Random rnd = new Random(0);
        foreach (Node v in m_Nodes)
        {
            v.Position.X = rnd.Next(maxwidth);
            v.Position.Y = rnd.Next(maxheight);
        }

        for (int i = 0; i < itercount; i++)
        {
            Invalidate();
            Application.DoEvents();
            System.Threading.Thread.Sleep(1000);

            foreach (Node v in m_Nodes)
            {
                v.Disp = Point.Empty;
                foreach (Node u in m_Nodes)
                {
                    if (v != u)
                    {
                        Point delta = new Point();
                        delta.X = v.Position.X - u.Position.X;
                        delta.Y = v.Position.Y - u.Position.Y;

                        double l = Math.Sqrt(
                            Math.Pow(delta.X, 2.0) + Math.Pow(delta.Y, 2.0)
                        );

                        double d = fr(l, k) / l;

                        v.Disp.X += (int)(delta.X * d);
                        v.Disp.Y += (int)(delta.Y * d);
                    }
                }
            }

            foreach (Edge r in m_Edges)
            {
                Point delta = new Point();
                delta.X = r.From.Position.X - r.To.Position.X;
                delta.Y = r.From.Position.Y - r.To.Position.Y;

                double l = Math.Sqrt(
                    Math.Pow(delta.X, 2) + Math.Pow(delta.Y, 2)
                );

                double d = fa(l, k) / l;

                r.From.Disp.X -= (int)(delta.X * d);
                r.From.Disp.Y -= (int)(delta.Y * d);

                r.To.Disp.X += (int)(delta.X * d);
                r.To.Disp.Y += (int)(delta.Y * d);
            }

            foreach (Node v in m_Nodes)
            {
                double l = Math.Sqrt(
                    Math.Pow(v.Disp.X, 2.0) + Math.Pow(v.Disp.Y, 2.0)
                );

                double d = Math.Min(l, temperature) / l;
                v.Position.X += (int)(v.Disp.X * d);
                v.Position.Y += (int)(v.Disp.Y * d);

                v.Position.X = Math.Min(
                    maxwidth,
                    (Math.Max(1, v.Position.X))
                );
                v.Position.Y = Math.Min(
                    maxheight,
                    (Math.Max(1, v.Position.Y))
                );
            }
            t -= temperature / itercount;
        }
    }
}