Вычислительная геометрия

using System;
using System.Collections.Generic;
using System.Drawing;
using System.Drawing.Drawing2D;
using System.Linq;
using System.Windows.Forms;
 
namespace ComputeGeometry
{
    public partial class MainForm : Form
    {
        #region Свойства
 
        private List<PointF> _points;
        private List<Tuple<PointF, PointF>> _connectPoints;
        const float R = 3;
 
        #endregion
 
        public MainForm()
        {
            InitializeComponent();
            rNumericUpDown.Maximum = decimal.MaxValue;
            panel1.Paint += panel1_Paint;
            Load += Form1_Load;
        }
 
        private void Form1_Load(object sender, EventArgs e)
        {
            rNumericUpDown.Value = Math.Min(panel1.ClientSize.Width, panel1.ClientSize.Height) - 15;
        }
 
        private void panel1_Paint(object sender, PaintEventArgs e)
        {
            if (_points == null || _connectPoints == null || _points.Count == 0)
            {
                return;
            }
            e.Graphics.TranslateTransform(
                splitContainer1.Panel2.ClientSize.Width / 2f, splitContainer1.Panel2.ClientSize.Height / 2f);
            e.Graphics.SmoothingMode = SmoothingMode.HighQuality;
 
            //Рисование соединений
            foreach (var tuple in _connectPoints)
            {
                e.Graphics.DrawLine(Pens.Red, tuple.Item1, tuple.Item2);
            }
            //Рисование точек
            foreach (var pt in _points)
            {
                e.Graphics.FillEllipse(Brushes.Blue, pt.X - R, pt.Y - R, 2 * R, 2 * R);
            }
        }
 
        private void buildUpButton_Click(object sender, EventArgs e)
        {
            int n = (int)nNumericUpDown.Value;
            float r = (int)rNumericUpDown.Value;
            //Левая верхняя вершина квадрата
            _points = new List<PointF>
                      {
                          new PointF(
                              -r / 2, -(float) (rNumericUpDown.Value / 2))
                      };
            //Расстояние между точками
            float dr = r / n;
            //векторы сторон квадрата
            var vectors = new List<Point>
                {
                    new Point(1, 0),
                    new Point(0, 1),
                    new Point(-1, 0),
                    new Point(0, -1)
                };
            var pt0 = _points[0];
 
            //Размещение точек на сторонах квадрата
            foreach (var vector in vectors)
            {
                for (int i = 1; i <= n; i++)
                {
                    var pt = pt0;
                    pt.X += i * dr * vector.X;
                    pt.Y += i * dr * vector.Y;
                    _points.Add(pt);
                }
                pt0 = _points.Last();
            }
            _points.RemoveAt(_points.Count - 1);
 
            //Пары точек для соединения
            _connectPoints = new List<Tuple<PointF, PointF>>();
            for (int i = 0; i < _points.Count; i++)
            {
                for (var j = i + 1; j < _points.Count; j++)
                {
                    if (IsFibonacci((j - i) / 4 * n))
                        _connectPoints.Add(new Tuple<PointF, PointF>(_points[i], _points[j]));
                }
            }
            panel1.Invalidate();
        }
 
        //Проверка, является ли число — числом Фибоначчи
        private static bool IsFibonacci(int n)
        {
            int a = 1, b = 1;
            if (n == 1)
            {
                return true;
            }
            while (b < n)
            {
                b += a;
                a = b - a;
            }
            return b == n;
        }
 
        private void rNumericUpDown_ValueChanged(object sender, EventArgs e)
        {
            buildUpButton.PerformClick();
        }
    }
}