結果

問題 No.453 製薬会社
ユーザー srup٩(๑`н´๑)۶
提出日時 2016-12-19 22:04:20
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 5 ms / 2,000 ms
コード長 1,473 bytes
コンパイル時間 1,434 ms
コンパイル使用メモリ 157,804 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-14 09:00:58
合計ジャッジ時間 1,994 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 9
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef vector<int> vint;
typedef pair<int,int> pint;
typedef vector<pint> vpint;
#define rep(i,n) for(int i=0;i<(n);i++)
#define reps(i,f,n) for(int i=(f);i<(n);i++)
#define each(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++)
#define all(v) (v).begin(),(v).end()
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))
const int MOD = 1e9 + 7;
const int INF = 1e9;

double c, d;
double ans;
double search(double a){
    double tc = c, td = d;
    double ret = 0.0;
    if(tc >= a * 3.0 / 4.0 && td >= a * 1.0 / 4.0){
        tc -= a * 3.0 / 4.0;
        td -= a * 1.0 / 4.0;
        ret += a * 1000.0;
    }else{
        //aの時に薬品Cが作れないなら、合計も0に
        return 0;
    }

    //残りはすべてBを作る。
    ret += 2000.0 * (min(tc * 7.0 / 2.0, td * 7.0 / 5.0));
    ans = ret;
    return ret;
}

int main(void){
    cin >> c >> d;
    //下凸関数を三分探索

    double l = 0.0;
    double r = 2000.0;
    rep(i, 10000){
        //m は、製品Aを作る量
        double m1 = (l * 2.0 + r) / 3.0;
        double m2 = (l + r * 2.0) / 3.0;
        if(search(m1) < search(m2)){//右が高い
            l = m1;
        }else{//左側が高い
            r = m2;
        }
    }
    printf("%.9f\n", ans);
    return 0;
}
0