Решение задачи линейного программирования графическим методом - C++ - Ответ 3114492

#include <stdio.h>
#include <values.h>
#include <stdlib.h>
#include <math.h>
#include <conio.h>
#include <dos.h>
#include <graphics.h>
int gd=0,gm;
 
int main()
{clrscr();
 int fmovx=0,fmovy=0,m,n,num=0,i;
 printf("Введите количество ограничивающих функций=>\n");
 scanf("%i",&n);
 double z[2],k[50][3],t[100][2],tud[100][2],buf[50][2];
 printf("Введите коэффициенты, где a*x1+b*x2>=c\n");
 for(i=0;i<n;i++)
    {
     printf("a[%i]=",i+1);scanf("%lf",&k[i][0]);
     printf("b[%i]=",i+1);scanf("%lf",&k[i][1]);
     printf("c[%i]=",i+1);scanf("%lf",&k[i][2]);
    }
 
    printf("Введите коэффициенты целевой ф-ии, где Z=a*x1+b*x2\n");
    printf("a=");scanf("%lf",&z[0]);
    printf("b=");scanf("%lf",&z[1]);
 
    printf("1-задача на min,2-задача на max\n");
    scanf("%i",&m);
 
/*-------------------Поиск точек, состовляющих область---------------------*/
 num=0;
 for(i=0;i<n-1;i++)
     for(int j=i+1;j<n;j++)
         {
 
         if(k[i][0]*k[j][1]!=k[i][1]*k[j][0])
             {
                t[num][0]=(k[i][2]*k[j][1]-k[j][2]*k[i][1])/(k[i][0]*k[j][1]-k[j][0]*k[i][1]);
                t[num][1]=(k[j][2]*k[i][0]-k[j][0]*k[i][2])/(k[i][0]*k[j][1]-k[j][0]*k[i][1]);
             }
         num++;
         }
 
 
    int No=0;
    for(i=0;i<num;i++)
    {
    int z=0;
    for(int j=0;j<n;j++)
        if(k[j][0]*t[i][0]+k[j][1]*t[i][1]-k[j][2]>=0) z++;
 
         if(z==n)
            {
             tud[No][0]=t[i][0];
             tud[No][1]=t[i][1];
             No++;
            }
    }
 
 
 
    for(i=0;i<No;i++)
        {buf[0][0]=tud[i][0];
         buf[0][1]=tud[i][1];
         for(int j=i+1;j<No;j++)
            if(buf[0][0]==tud[j][0]&&buf[0][1]==tud[j][1])
                {
                tud[j][0]=-1;
                tud[j][1]=-1;
                }
        }
 
 
    int j=0;
    for(i=0;i<No;i++)
        if(tud[i][0]!=-1)
        {
            buf[j][0]=tud[i][0];
            buf[j][1]=tud[i][1];
            j++;
        }
 
/*-----------------------------------------------------------------*/
 
for(i=0;i<j;i++)
printf("\nud1=%f ud2=%f",buf[i][0],buf[i][1]);
 
 double minmax;
 int tt=0;
 
 if(m==1)
 {
 minmax=10e308;
 for(i=0;i<j;i++)
     if(z[0]*buf[i][0]+z[1]*buf[i][1]<minmax)
                {minmax=z[0]*buf[i][0]+z[1]*buf[i][1];tt=i;}
 }
 
 if(m==2)
 {
 minmax=-10e308;
 for(i=0;i<j;i++)
     if(z[0]*buf[i][0]+z[1]*buf[i][1]>minmax)
                {minmax=z[0]*buf[i][0]+z[1]*buf[i][1];tt=i;}
 }
 
 printf("\n mm=%lf",minmax);
 getch();
/*-----------------------------------------------------------------*/
 initgraph(&gd,&gm,"c\\bc5\\bgi");
 
 setcolor(15);
 line(320,0,320,480);
 line(0,240,640,240);
 //system koordinat
 int edot=20;
 for(i=1;i<=480/edot;i++)
         {
         line(318,240-edot*i,322,240-edot*i);
         line(318,240+edot*i,322,240+edot*i);
         }
 
 for(i=1;i<=640/edot;i++)
         {
         line(320-edot*i,238,320-edot*i,242);
         line(edot*i+320,238,edot*i+320,242);
         }
 int x1,x2,y1,y2;
 x1=320/edot;
 x2=-320/edot;
 
 //draw line
 setcolor(8);
 for(i=0;i<n;i++)
    {
        if(k[i][0]==0&&k[i][1]!=0)
         {
            moveto(0,240-edot*k[i][2]/k[i][1]);
            lineto(640,240-edot*k[i][2]/k[i][1]);
         } else
        if(k[i][0]!=0&&k[i][1]==0)
         {
            moveto(320-edot*k[i][2]/k[i][0],0);
            lineto(320-edot*k[i][2]/k[i][0],480);
         } else
 
         {
            moveto(0,240-edot*(k[i][2]+320/edot*k[i][0])/k[i][1]);
            lineto(640,240-edot*(k[i][2]-320/edot*k[i][0])/k[i][1]);
         }
    }
 setcolor(BROWN);
 int *pol;pol=new int[2*j+2];
 for(i=0;i<2*j;i+=2)
     {pol[i]=320+buf[i/2][0]*edot;pol[i+1]=240-buf[i/2][1]*edot;}
        pol[2*j]=pol[0];pol[2*j+1]=pol[1];
 char s[10];
 //fill oblast
 setfillstyle(1,BROWN);
 fillpoly(j+1,pol);
 
 setcolor(CYAN);
 line(320,240,z[0]*edot+320,240-z[1]*edot);
 line(z[0]*edot+320,240-z[1]*edot,z[0]*edot+320-3*fabs(z[0])/z[0],240-z[1]*edot);
 line(z[0]*edot+320,240-z[1]*edot,z[0]*edot+320,240+3*fabs(z[1])/z[1]-z[1]*edot);
 
 setcolor(RED);
 line(2*buf[tt][0]*edot+fabs(z[0])/z[0]*50+320,240-2*buf[tt][1]*edot+fabs(z[1])/z[1]*50,320+buf[tt][0]*edot-fabs(z[0])/z[0]*50,240-buf[tt][1]-fabs(z[1])/z[1]*50);
 
 
 getch();
 closegraph;  
}