結果
問題 | No.453 製薬会社 |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
#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;}