线性规划
int find(){ LDB maxi=-1e15;int po; for (int i=1;i<=ycnt;i++) if (ans[i]>maxi){ maxi=ans[i];po=i; } if (maxi<eps) return(0); int tar=po;maxi=1e15; for (int i=1;i<=xcnt;i++) if ((lim[i][tar])>eps&&lim[i][0]/lim[i][tar]<maxi){ maxi=lim[i][0]/lim[i][tar];po=i; } if (maxi==1e15){ ans[0]=-1e12;return(0); } for (int i=0;i<=ycnt;i++){ LDB t=lim[po][i];lim[po][i]=lim[1][i];lim[1][i]=t; } return(tar); } void work(){ while (int po=find()){ for (int i=2;i<=xcnt;i++) if (fabs(lim[i][po])>eps){ LDB bas=lim[i][po]/lim[1][po]; lim[i][po]=0; for (int j=0;j<=ycnt;j++) lim[i][j]-=lim[1][j]*bas; } LDB bas=ans[po]/lim[1][po]; ans[po]=0; for (int j=0;j<=ycnt;j++) ans[j]-=lim[1][j]*bas; } }
注意po那一列需要消2倍