单纯形模板

最新推荐文章于 2019-09-16 20:42:59 发布

ACgay_caoyue

最新推荐文章于 2019-09-16 20:42:59 发布

阅读量352

点赞数

本文链接：https://blog.csdn.net/KIKO_caoyue/article/details/89929481

版权

本文详细介绍了运筹学中的单纯形法，这是一种解决线性规划问题的有效算法。通过阐述单纯形法的基本原理和步骤，帮助读者理解如何运用单纯形模板来求解实际问题。此外，还探讨了单纯形法的优缺点以及在优化领域的应用。

摘要由CSDN通过智能技术生成

#include <bits/stdc++.h>
#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;
        doub