#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define rep(i,n) for(int i = 0;i=0;i--) #define ALL(a) a.begin(),a.end() template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } typedef long long ll; typedef long double ld; struct Simplex { using F = double; static constexpr F EPS = 1e-9; int N, M; vector idx;vector> a; int i_ch, j_ch; void initialize(const vector>& A, const vector& b, const vector& c) { N = c.size(), M = A.size(); a.assign(M + 2, vector(N + 2)); i_ch = M; rep(i,M) { rep(j,N) a[i][j] = -A[i][j]; a[i][N] = 1, a[i][N + 1] = b[i]; if (a[i_ch][N + 1] > a[i][N + 1]) i_ch = i; }rep(j,N) a[M][j] = c[j]; a[M + 1][N] = -1; idx.resize(N + M + 1);iota(ALL(idx), 0); }void solve() { vector jupd; for (j_ch = N;;) { if (i_ch < M) { swap(idx[j_ch], idx[i_ch + N + 1]); a[i_ch][j_ch] = 1 / a[i_ch][j_ch]; jupd.clear(); rep(j, N+2) { if (j != j_ch) { a[i_ch][j] *= -a[i_ch][j_ch]; if (abs(a[i_ch][j]) > EPS) jupd.push_back(j); } }rep(i,M+2) { if (abs(a[i][j_ch]) < EPS or i == i_ch) continue; for (auto j : jupd) a[i][j] += a[i][j_ch] * a[i_ch][j]; a[i][j_ch] *= a[i_ch][j_ch]; } }j_ch = -1; rep(j,N+1) { if (j_ch < 0 or idx[j_ch] > idx[j]) { if (a[M + 1][j] > EPS or (abs(a[M + 1][j]) < EPS and a[M][j] > EPS)) j_ch = j; } }if (j_ch < 0) break; i_ch = -1; rep(i,M) { if (a[i][j_ch] < -EPS) { if (i_ch < 0) { i_ch = i; }F com = a[i_ch][N + 1] / a[i_ch][j_ch] - a[i][N + 1] / a[i][j_ch]; if (com < -EPS) i_ch = i; else if (com < EPS && idx[i_ch] > idx[i]) i_ch = i; } }if (i_ch < 0) { is_infty = 1; break; } }if (a[M + 1][N + 1] < -EPS) { infeasible = 1; return; }x.assign(N, 0); rep(i,M) { if (idx[N + 1 + i] < N) x[idx[N + 1 + i]] = a[i][N + 1]; }ans = a[M][N + 1]; } Simplex(const vector>& A, const vector& b, const vector& c) { is_infty = infeasible = false; initialize(A, b, c); solve(); }bool is_infty,infeasible; vector x; F ans; }; int main() { double c, d; cin >> c >> d; vector> a = { {0.75,2.0/7.0},{0.25,5.0 / 7.0} }; vector b = { c,d }, C = { 1000,2000 }; Simplex sim(a, b, C); cout <