单纯形法的c语言实现

单纯形法的c语言实现

#include<stdio.h>
#define m 3//行
#define n 5//列
int min(float a[n],int l){
float min=+9999;
int i,t;
for(i=0;i<l;i++)
{
if(a[i]<min)
{
min=a[i];
t=i;
}
}
return t;
}

int pd(float a[n]){
int i,j=0;
for(i=0;i<n;i++)
if(a[i]>=0)
j++;
if(j==n)
return 0;
else
return 1;
}

void main(){
float s[m][n+3]={{0.25,0.5,1,0,0,120,3,0},{0.5,0.5,0,1,0,150,4,0},{0.25,0,0,0,1,50,5,0}};
float x[5]={-12,-15,0,0,0};
float y[n+1];
float q[m];
int i,j,my,mq;
float tem;
float chu,ru,cs;
printf("------------------------------------------------------------\n");
for(i=0;i<n;i++)
y[i]=-1;
printf(“cb xb b x1 x2 x3 x4 x5 “);
while(pd(y))
{
printf(”\n------------------------------------------------------------\n”);
for(i=0;i<n;i++)
{
tem=0;
for(j=0;j<m;j++)
tem=tem+s[j][n+2]*s[j][i];
y[i]=x[i]-tem;
//printf("%f",y[i]);
}
for(i=0;i<3;i++)
{
printf("%-8.2f",s[i][n+2]);
printf("%-8.2f",s[i][n+1]);
printf("%-8.2f",s[i][n]);
for(j=0;j<5;j++)
printf("%-8.2f",s[i][j]);
printf("\n");
}
printf(" -z: “);
y[n]=0;
for(i=0;i<m;i++)
y[n]=y[n]+s[i][n]*s[i][n+2];
printf(”%-8.2f",-y[n]);

	for(i=0;i<5;i++)
	{
		printf("%-8.2f",y[i]);
	}				
	my=min(y,5);
	ru=my+1;//确定入基变量
	//printf("%d",my);
	for(i=0;i<m;i++)
	{
		q[i]=s[i][5]/s[i][my];
	}
	mq=min(q,3);
	chu=s[mq][6];//确定出基变量
	s[mq][n+1]=ru;
	s[mq][n+2]=x[(int)ru-1];
	cs=1/s[mq][my];
	for(i=0;i<n+1;i++)
	{
		s[mq][i]=s[mq][i]*cs;
	}
	for(i=0;i<m;i++)
	{
		if(i!=mq&&s[i][my]!=0)
		{
			tem=s[i][my];
			for(j=0;j<n+1;j++)
			{
				s[i][j]=s[i][j]-(s[mq][j]*tem);
				//printf("%f\n",tem);
				
			}
		}
	}		
}
printf("\n------------------------------------------------------------\n");
printf("最优解为:");
for(i=0;i<m;i++)
	printf("x%.0f=%.2f ",s[i][n+1],s[i][n]);
printf("\n最优值为:%.2f\n",-y[n]);

}