UOJ#179. 线性规划[模板]

传送门 http://uoj.ac/problem/179

震惊,博主竟然还不会线性规划!

单纯形实在学不会啊……背个板子当黑盒用……

学(chao)了NanoApe dalao的板子

 

 1 #include<iostream>
 2 #include<algorithm>
 3 #include<cstring>
 4 #include<cstdio>
 5 #include<cmath>
 6 using namespace std;
 7 const double eps=1e-8;
 8 const int mxn=110;
 9 int read(){
10     int x=0,f=1;char ch=getchar();
11     while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
12     while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
13     return x*f;
14 }
15 
16 int n,m,t;
17 double c[mxn];
18 double a[mxn][mxn];
19 int idx[mxn],idy[mxn];
20 int st[mxn],top=0;
21 void Pivot(int x,int y){//用idy代换idx
22     swap(idy[x],idx[y]);
23     double tmp=a[x][y];a[x][y]=1/a[x][y];
24     int i,j;top=0;
25     for(i=0;i<=n;i++)if(y!=i)a[x][i]/=tmp;
26     for(i=0;i<=n;i++)if((y!=i) && fabs(a[x][i])>eps)st[++top]=i;
27     for(i=0;i<=m;i++){
28         if((i==x)||(fabs(a[i][y])<eps))continue;
29         for(j=1;j<=top;j++)a[i][st[j]]-=a[x][st[j]]*a[i][y];
30         a[i][y]=-a[i][y]/tmp;
31     }
32     return;
33 }
34 void Simplex(){
35     int i,j;
36     while(1){
37         int x=0,y=0;double mn=1e15;
38         for(i=1;i<=n;i++)if(a[0][i]>eps){y=i;break;}
39         if(!y)break;
40         for(i=1;i<=m;i++){
41             if(a[i][y]>eps && (a[i][0]/a[i][y]<mn)){
42                 mn=a[i][0]/a[i][y];
43                 x=i;
44             }
45         }
46         if(!x){printf("Unbounded\n");exit(0);}
47         Pivot(x,y);
48     }
49     return;
50 }
51 int main(){
52     int i,j;
53     n=read();m=read();t=read();
54     for(i=1;i<=n;i++)a[0][i]=read();//目标函数
55     for(i=1;i<=m;i++){
56         for(j=1;j<=n;j++)a[i][j]=read();
57         a[i][0]=read();
58     }
59     for(i=1;i<=n;i++)idx[i]=i;//基变量 
60     for(i=1;i<=m;i++)idy[i]=i+n;//非基变量 
61     while(1){
62         int x=0,y=0;
63         for(i=1;i<=m;i++)if(a[i][0]<-eps && ((!x)||(rand()&1)))x=i;
64         if(!x)break;
65         for(i=1;i<=n;i++)if(a[x][i]<-eps && ((!y)||(rand()&1)))y=i;
66         if(!y){printf("Infeasible\n");return 0;}
67         Pivot(x,y);
68     }
69     Simplex();
70     printf("%.8f\n",-a[0][0]);
71     if(!t)return 0;
72     for(i=1;i<=n;i++)a[0][i]=0;
73     for(i=1;i<=m;i++)if(idy[i]<=n)a[0][idy[i]]=a[i][0];
74     for(i=1;i<=n;i++)printf("%.8f ",a[0][i]);puts("");
75     return 0;
76 }

 

posted @ 2017-05-18 21:00  SilverNebula  阅读(363)  评论(0编辑  收藏  举报
AmazingCounters.com