oporniy_plan/Опорные планы/Plans.cs

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;

namespace Опорные_планы
{
    internal class Plans
    {
        public int[,]? InputArray()
        {
            Console.WriteLine("\nВведите кол-во складов: ");
            int storage = Convert.ToInt32(Console.ReadLine());
            Console.WriteLine("\nВведите кол-во потребителей: ");
            int consumer = Convert.ToInt32(Console.ReadLine());
            int[,] postavki = new int[storage + 1, consumer + 1];

            for (int i = 0; i <= storage; i++)
            {
                for (int j = 0; j <= consumer; j++)
                {
                    if (i == 0 && j == 0)
                    {
                        postavki[i, j] = 0;
                    }
                    else if (j == 0)
                    {
                        Console.WriteLine("\nВведите кол-во груза для складов:");
                        postavki[i, j] = int.Parse(Console.ReadLine());
                    }
                    else if (i == 0)
                    {
                        Console.WriteLine("\nВведите кол-во груза для поставщиков:");
                        postavki[i, j] = int.Parse(Console.ReadLine());
                    }
                }
            }

            for (int i = 1; i <= storage; i++)
            {
                for (int j = 1; j <= consumer; j++)
                {
                    Console.WriteLine("\nВведите тарифный план:");
                    postavki[i, j] = int.Parse(Console.ReadLine());
                }
            }

            int sumStorage = 0;
            int sumConsumer = 0;

            for (int i = 0; i <= storage; i++) //Проверка на открытость
            {
                for (int j = 0; j <= consumer; j++)
                {
                    sumStorage = sumStorage + postavki[i, 0];
                    sumConsumer = sumConsumer + postavki[0, j];
                }
            }

            if (sumStorage == sumConsumer)
            {
                Console.WriteLine("\nЗадача закрытая\n");
                Console.WriteLine("Начальные данные: ");
                for (int i = 0; i <= storage; i++) //Вывод массива

                {
                    for (int j = 0; j <= consumer; j++)
                    {
                        Console.Write("[" + i + "," + j + "]: " + postavki[i, j] + "\t");
                    }
                    Console.WriteLine();
                }
                Console.WriteLine();
                return postavki;
            }
            else
            {
                Console.WriteLine("\nЗадача открытая\n");
                for (int i = 0; i <= storage; i++) //Вывод массива

                {
                    for (int j = 0; j <= consumer; j++)
                    {
                        Console.Write("[" + i + "," + j + "]: " + postavki[i, j] + "\t");
                    }
                    Console.WriteLine();
                }
                Console.WriteLine();
                return postavki;
            }

        }

        public static bool IsOptimal(int?[,] opornyPlan, double?[] potenciali_u, double?[] potenciali_v, int[,] postavki)
        {
            for (int i = 0; i < opornyPlan.GetLength(0); i++)
            {
                for (int j = 0; j < opornyPlan.GetLength(1); j++)
                {
                    if (opornyPlan[i, j] == 0 && potenciali_u[i] + potenciali_v[j] > postavki[i, j])
                    {
                        return false;
                    }
                }
            }

            return true;
        }

        public int?[,] Northwest(int[,] postavki)
        {
            int?[,] opornyPlan = new int?[postavki.GetLength(0), postavki.GetLength(1)];
            int counter = 0;
            for (int i = 1; i < postavki.GetLength(0); i++) //Заполнение опорного плана
            {
                for (int j = 1; j < postavki.GetLength(1); j++)
                {
                    if (postavki[i, 0] >= postavki[0, j]) //Условие, если кол-во груза на складе больше
                    {
                        opornyPlan[i, j] = postavki[0, j];
                        postavki[i, 0] = postavki[i, 0] - postavki[0, j]; //Остаток
                        postavki[0, j] = 0;
                        counter++;
                    }
                    else if (postavki[i, 0] < postavki[0, j]) //Условие, если кол-во груза для потребителя больше
                    {
                        opornyPlan[i, j] = postavki[i, 0];
                        postavki[0, j] = postavki[0, j] - postavki[i, 0]; //Остаток
                        postavki[i, 0] = 0;
                        counter++;
                    }
                }
            }

            for (int i = 1; i < postavki.GetLength(0); i++)
            {
                for (int j = 1; j < postavki.GetLength(1); j++)
                {
                    Console.Write("[" + i + "," + j + "]: " + opornyPlan[i, j] + "\t");
                }
                Console.WriteLine();
            }
            Console.WriteLine();

            int? Lx = 0;
            for (int i = 0; i < postavki.GetLength(0); i++)
            {
                for (int j = 0; j < postavki.GetLength(1); j++)
                {
                    if (i != 0 && j != 0 && opornyPlan[i, j] != 0)
                    {
                        Lx = Lx + postavki[i, j] * opornyPlan[i, j];
                    }
                }
            }
            Console.WriteLine("Целевая функция равна: " + Lx + "");


            int numNonZeroRows = 0;
            int numNonZeroCols = 0;
            for (int i = 1; i < postavki.GetLength(0); i++)
            {
                for (int j = 1; j < postavki.GetLength(1); j++)
                {
                    if (opornyPlan[i, j] != null && opornyPlan[i, j] > 0)
                    {
                        numNonZeroRows++;
                        break; 
                    }
                }
            }

            for (int j = 1; j < postavki.GetLength(1); j++)
            {
                for (int i = 1; i < postavki.GetLength(0); i++)
                {
                    if (opornyPlan[i, j] != null && opornyPlan[i, j] > 0)
                    {
                        numNonZeroCols++;
                        break; 
                    }
                }
            }

            if (counter == (opornyPlan.GetLength(0) - 1 + opornyPlan.GetLength(1) - 1 - 1))
            {
                bool isNonDegenerate = numNonZeroRows == opornyPlan.GetLength(0) - 1 && numNonZeroCols == opornyPlan.GetLength(1) - 1;

                if (!isNonDegenerate)
                {
                    Console.WriteLine("\nОпорный план вырожденный, переделайте его!");
                    Northwest(InputArray());
                }
            }
            else
            {
                Console.WriteLine("\nОпорный план невырожденный\n");
            }

            bool Check = MethodOfPotentials.CalculatePotentials(opornyPlan, postavki);

            if (Check == false)
            {
                Console.WriteLine("\nВведите новый опорный план:");
                Northwest(InputArray());
            }

            return opornyPlan;
        }
        public int[,] Fogel_approximations(int[,] postavki)
        {
            Console.WriteLine("\nВ разработке\n");
            return postavki;
        }


        public void ChoiseMethod()
        {
            while (true)
            {
                Console.WriteLine("Введите метод, с помощью которого хотите вычислить опорный план: \n" +
                    "1 - Метод северо-западного угла\n" +
                    "0 - Выход из программы");
                int Choice = Convert.ToInt32(Console.ReadLine());
                if (Choice == 2)
                {
                    Console.WriteLine("\nВ разработке\n");
                }
                else if (Choice == 1)
                {
                    Northwest(InputArray());
                }
                else if (Choice == 3)
                {
                    Fogel_approximations(InputArray());
                }
                else if (Choice == 0)
                {
                    break;
                }
                else
                {
                    Console.WriteLine("Метод не существует!");
                }
            }
        }
    }
}