#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define REP(i,a,n) for(int i=(a); i<(int)(n); i++) #define rep(i,n) REP(i,0,n) #define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it) #define ALLOF(c) (c).begin(), (c).end() typedef long long ll; class Simplex { int N; int M; std::vector< std::vector< double > > tbl; void pivot(int p, int q){ for(int j=0; j<=N; j++){ for(int k=M; k>=0; k--){ if(j != p && k != q){ tbl[j][k] -= tbl[p][k] * tbl[j][q] / tbl[p][q]; } } } for(int j=0; j<=N; j++){ if(j != p){ tbl[j][q] = 0; } } for(int k=0; k<=M; k++){ if(k != q){ tbl[p][k] = tbl[p][k] / tbl[p][q]; } } tbl[p][q] = 1; } public: Simplex(int N, int M): N(N), M(M), tbl(N+1, std::vector(M+1, 0.0)) { if(N > M){ std::cerr << "check: N > M" << std::endl; } } int getN(){ return N; } int getM(){ return M; } double get(int i, int j){ if(0<=i && i=0); q++); for(p=0; (p<=N) && (tbl[p][q]<=0); p++); if(q > M || p > N) break; //std::cerr << "!" << p << " " << q << std::endl; for(int i=p+1; i<=N; i++){ if(tbl[i][q] > 0){ if(tbl[i][M] / tbl[i][q] < tbl[p][M] / tbl[p][q]){ p = i; } } } pivot(p, q); } if(q == M+1) return true; if(p == N+1) return false; return false; } }; int main(){ double C, D; cin >> C >> D; Simplex smp(2, 4); smp.set(0, 0, -1000); smp.set(0, 1, -2000); smp.set(0, 2, 0); smp.set(0, 3, 0); smp.set(0, 4, 0); smp.set(1, 0, 3/4.0); smp.set(1, 1, 2/7.0); smp.set(1, 2, 1); smp.set(1, 3, 0); smp.set(1, 4, C); smp.set(2, 0, 1/4.0); smp.set(2, 1, 5/7.0); smp.set(2, 2, 0); smp.set(2, 3, 1); smp.set(2, 4, D); smp.solve(); printf("%.12lf\n", smp.get(0,4)); return 0; }