[单纯形] UOJ#179. 线性规划

本文介绍了一种解决线性规划问题的方法,通过单纯形法实现求解,包括初始化过程及迭代优化步骤,适用于变量和约束数量较小的情况。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题意

解线性规划,以标准型给出。n个变量, m个约束。
n,m20

题解

单纯形模板题

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=55;
const double eps=1e-8;
int n,m,Type;
double a[maxn][maxn],ans[maxn];
int id[maxn*2];

void Pivot(int l,int e){
    swap(id[l+n],id[e]);
    double t=a[l][e]; a[l][e]=1; 
    for(int i=0;i<=n;i++) a[l][i]/=t;
    for(int i=0;i<=m;i++) if(i!=l&&fabs(a[i][e])>eps){
        t=a[i][e]; a[i][e]=0; 
        for(int j=0;j<=n;j++) a[i][j]-=t*a[l][j]; 
    }
}
bool Init(){
    while(1){
        int l=0,e=0;
        for(int i=1;i<=m;i++) if(a[i][0]<-eps&&(!l||(rand()&1))) l=i;
        if(!l) return true;
        for(int i=1;i<=n;i++) if(a[l][i]<-eps&&(!e||(rand()&1))) e=i;
        if(!e) return printf("Infeasible\n"),false; 
        Pivot(l,e);
    }
    return true;
}
bool Simplex(){
    while(1){
        int l=0,e=0; double _min=1e+50;
        for(int i=1;i<=n;i++) if(a[0][i]>eps){ e=i; break; }
        if(!e) break;
        for(int i=1;i<=m;i++) 
         if(-a[i][e]<-eps&&a[i][0]/a[i][e]<_min) _min=a[i][0]/a[i][e], l=i; 
        if(!l) return printf("Unbounded\n"),false; 
        Pivot(l,e);
    }
    return true;
}
int main(){
    freopen("uoj179.in","r",stdin);
    freopen("uoj179.out","w",stdout);
    scanf("%d%d%d",&n,&m,&Type);
    for(int i=1;i<=n;i++) scanf("%lf",&a[0][i]);
    for(int i=1;i<=m;i++){
        for(int j=1;j<=n;j++) scanf("%lf",&a[i][j]);
        scanf("%lf",&a[i][0]);
    }
    for(int i=1;i<=n;i++) id[i]=i; 
    if(Init()&&Simplex()){
        printf("%.8lf\n",-a[0][0]); 
        if(Type){
            for(int i=1;i<=m;i++) ans[id[i+n]]=a[i][0]; 
            for(int i=1;i<=n;i++) printf("%.8lf ",ans[i]);
        }
    }
    return 0;
} 
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值