单纯形模板

#include 
#define INF 0x3f3f3f3f
const int maxn=1e3+7,maxm=1e4+7;
int n,m,x,y,p,id,k;
double ans,c[maxn],b[maxm],a[maxm][maxn];  //a数组是方程组系数矩阵, b数组是目标函数系数矩阵, c数组是增广矩阵最后一列 
void pvt(int id,int p) 
{
    a[id][p]=1/a[id][p];
	b[id]*=a[id][p];
    for(int i=1;i<=n;i++) 
	{
		if(i^p) a[id][i]*=a[id][p];
	}
    for(int i=1;i<=m;i++)
    {
    	if((i^id) && a[i][p]) 
		{
	        for(int j=1;j<=n;j++) 
			{
				if(j^p) a[i][j]-=a[i][p]*a[id][j];
			}
	        b[i]-=a[i][p]*b[id];
			a[i][p]*=-a[id][p];
	    }
	}
    for(int i=1;i<=n;i++) 
	{
		if(i^p) c[i]-=c[p]*a[id][i];
	}
    ans+=c[p]*b[id];
	c[p]*=-a[id][p];
}
double solve() 
{
    while(true)
	{
        for(p=1;p<=n;p++) 
		{
			if(c[p]>0) break;
		}
        if(p==n+1) return ans;
        double fz=INF;
        for(int i=1;i<=m;i++) 
		{
			if(a[i][p]>0 && b[i]/a[i][p]