挺繁琐，之前发的帖子，分值太低现在重发一个。欢迎各位大侠~~
#include "stdlib.h"
#include "math.h"
#include "stdio.h"
float objfx(float x[]);
void constraint(float x[],float g[]);
int gau(float x[],float g[],int kg)
{ int i;
constraint(x,g);
for(i=0;i<kg;i++)
{
if(g[i]<0)
goto s333;
}
return 1;
s333:return 0;
}
void xcent(int n,int k,int ll,int lh,float x0[],float xcom[][100])
{ int i,l;
float xs;
for(i=0;i<n;i++)
{
xs=0;
for(l=0;l<ll;l++)
{
if(l!=lh)
xs=xs+xcom[i][l];
}
if(lh>-1)
x0[i]=xs/(ll-1);
else
x0[i]=xs/ll;
}
}
void fxse(int n,int k,float x[],float xcom[][100],float fxk[])
{
int l,lp,lp1,i;
float w;
for(l=0;l<k-1;l++)
for(lp=0;lp<k-l;lp++)
{
lp1=lp+1;
if(fxk[lp]<=fxk[lp1])
{
w=fxk[lp];
fxk[lp]=fxk[lp1];
fxk[lp1]=w;
for(i=0;i<n;i++)
{
