using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Windows;
using System.Windows.Controls;
using System.Windows.Data;
using System.Windows.Documents;
using System.Windows.Input;
using System.Windows.Media;
using System.Windows.Media.Imaging;
using System.Windows.Navigation;
using System.Windows.Shapes;


namespace WpfApplication3
{
    /// <summary>
    /// Логика взаимодействия для Window1.xaml
    /// </summary>
    public partial class Window1 : Window
    {
        public Window1()
        {
            InitializeComponent();
        }

        private void btnClose_Click(object sender, RoutedEventArgs e)
        {
            Close();
        }
        private int Lz, P, C, Td, Nr;
        private double Stotn;
        private double[] MTobs = { 0, 0, 0, 0 };

        private void btn_Rashet_Click(object sender, RoutedEventArgs e)
        {
            Lz = int.Parse(txt_Lz.Text); //Среднее число заявок в день
            P = int.Parse(txt_P.Text); //Прибыль от единицы произведенной продукции
            C = int.Parse(txt_C.Text);  //Издержки
            Td = int.Parse(txt_Td.Text); //Время на выполнение всех заявок
            Stotn = double.Parse(txt_SKO.Text); //Среднее квадратичное отклонение
            Nr = int.Parse(txt_Np.Text); //Число реализаций

            MTobs[0] = double.Parse(txt_1.Text); //массивы цехов
            MTobs[1] = double.Parse(txt_2.Text);
            MTobs[2] = double.Parse(txt_3.Text);
            MTobs[3] = double.Parse(txt_4.Text);

           Model7();

        }
        public double[] TH = { 0, 0, 0, 0 };
        public double[] TK = { 0, 0, 0, 0 };
        public double Mprof, Cprof, Sum2, SNM, SigProf, Gprof;
        public double Tz, TzI, TzO, Prof, Disp, N, Tobs, fact, Min, Max, SMT;
        public int Ir, Nz, Nobs, i, j;

        public const int K = 4; //кол-во цехов
        private void Model7()
        {
            Mprof = 0;
            Sum2 = 0;
            for (Ir = 0; Ir < Nr; Ir++)
            {
                Nz = 0; //Количество поступивших заявок на выполнение
                TzO = 0;    //Время поступления предыдущей заявки
                Nobs = 0;   //Количество выполненных заявок
                //Обнуление времени работы цехов
                for (j = 0; j < K; j++)
                {
                    TK[j] = 0;
                    TH[j] = 0;
                    bool t = false;

                    do
                    {
                        Random r = new Random();
                        double z = r.NextDouble();

                        Tz = TzO - Math.Log(z, Math.E) / Lz; //Расчет конца выполнения поступившей заявки
                        TzO = Tz; //Обновляется время поступления будущей заявки
                        if (TzO > Td)
                        {
                            t = true;
                            break;
                        }
                        else
                        {
                            Nz += 1;
                            for (j = 0; j < K; j++)
                            {
                                if (Tz > TK[j])
                                {
                                    TH[j] = Tz;
                                }
                                else
                                {
                                    TH[j] = TK[j];
                                }
                                N = NORM();
                                Tobs = MTobs[j] * (1 + N * Stotn);
                                TK[j] = TK[j] + Tobs;
                                if (TK[j] > Td)
                                {
                                    goto a;
                                }
                                Tz = TK[j];
                            }
                            Nobs += 1;
                        }
                    } while (t == false);
                a: t = true;

                    Prof = P * Nobs - C;
                    Mprof = Mprof + Prof;
                    Sum2 = Sum2 + Prof * Prof;
                }

                Cprof = Mprof / Nr;
                Disp = (Sum2 - Nr * Cprof * Cprof) / (Nr - 1);
                if (Disp > 0)
                {
                    SigProf = Math.Sqrt(Disp);
                }
                else
                {
                    SigProf = 0;
                }
                Gprof = Cprof - 1.28 * SigProf;
                Factor();

                txt_Fact.Text = fact.ToString();
                txt_MinP.Text = Gprof.ToString();
            }
        }

        private double NORM()
        {
            double Sz = 0;
            for (int i = 0; i < 12; i++)
            {
                Random r = new Random();
                double z = r.NextDouble();
               Sz += z;
            }
            return Sz / 12;
        }

        private void Factor()
        {
            Min = MTobs.Min();
            Max = MTobs.Max();
            SMT = MTobs.Sum();
           
            fact = (Max - Min) / SMT;
        }

        private void button1_Click(object sender, RoutedEventArgs e)
        {
            txt_1.Clear();
            txt_2.Clear();
            txt_3.Clear();
            txt_4.Clear();
            txt_Lz.Clear();
            txt_P.Clear();
            txt_C.Clear();
            txt_Td.Clear();
            txt_Fact.Clear();
            txt_MinP.Clear();
            txt_Np.Clear();
            txt_SKO.Clear();

        }
		
	}
}