結果
問題 | No.453 製薬会社 |
ユーザー | char134217728 |
提出日時 | 2017-09-04 18:48:21 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 2,275 bytes |
コンパイル時間 | 1,312 ms |
コンパイル使用メモリ | 159,604 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-04-24 12:49:53 |
合計ジャッジ時間 | 2,056 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
5,248 KB |
testcase_01 | AC | 1 ms
5,376 KB |
testcase_02 | AC | 1 ms
5,376 KB |
testcase_03 | AC | 1 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 1 ms
5,376 KB |
testcase_06 | AC | 2 ms
5,376 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 1 ms
5,376 KB |
testcase_09 | AC | 2 ms
5,376 KB |
testcase_10 | AC | 1 ms
5,376 KB |
testcase_11 | AC | 2 ms
5,376 KB |
testcase_12 | AC | 2 ms
5,376 KB |
コンパイルメッセージ
main.cpp:77:1: warning: ISO C++ forbids declaration of ‘main’ with no type [-Wreturn-type] 77 | main(){ | ^~~~
ソースコード
#include <bits/stdc++.h> #define FOR(i,a,b) for (int i=(a);i<(b);i++) #define FORR(i,a,b) for (int i=(a);i>=(b);i--) #define pb push_back #define pcnt __builtin_popcount #define show(x) cout<<#x<<" = "<<x<<endl; #define maxs(x,y) x = max(x,y) #define mins(x,y) x = min(x,y) #define fi first #define se second #define rng(a) a.begin(),a.end() #define each(it,c) for(__typeof((c).begin()) it=(c).begin();it!=(c).end();it++) #define sz(x) (int)(x).size() using namespace std; typedef long long ll; typedef pair<int,int> pii; typedef vector<int> vi; typedef vector<vi> vvi; typedef vector<pii> vpii; typedef set<int> si; typedef pair<ll,ll> pll; typedef vector<ll> vl; typedef vector<vl> vvl; typedef vector<pll> vpll; typedef set<ll> sl; template<typename T>string join(vector<T>&v) {stringstream s;FOR(i,0,sz(v))s<<' '<<v[i];return s.str().substr(1);} ll gcd(ll a,ll b){if(a>b)swap(a,b);for(;a>0;b%=a,swap(a,b));return b;} int modpow(ll a,ll n,int m){if(a==0)return a;ll p=1;for(;n>0;n/=2,a=a*a%m)if(n&1)p=p*a%m;return(int)p;} void dout(double d){printf("%.12f\n",d);} const int iinf = 1e9; const ll linf = 1e18; const int mod = 1e9+7; const double pi = acos(-1); const double eps = 1e-10; double r[2], m[2], e[2][2] = { {0.75, 2./7}, {0.25, 5./7} }; const int GN=2; struct Gauss{ int solve(int n, double o[GN][GN], double u[GN], double r[GN]) { int q; double m[GN][GN],v[GN], p; memmove(m,o,sizeof(m)); memmove(v,u,sizeof(v)); FOR(i, 0, n) { p = 0.; q = 0; FOR(j, i, n) if(p < abs(m[i][j])){ p = abs(m[i][j]); q = j; } if(p < eps) return -1; swap(v[i],v[q]); FOR(k, i, n) swap(m[i][k], m[q][k]); v[i]/=m[i][i]; FORR(k, n-1, i) m[i][k]/=m[i][i]; FOR(j, i+1, n){ v[j]-=v[i]*m[j][i]; FORR(k, n-1, i) m[j][k]-=m[i][k]*m[j][i]; } } FORR(i, n-1, 0) { FORR(j, n-1, i+1) v[i]-=m[i][j]*v[j],m[i][j]=0; r[i]=v[i]; } return 0; } }; Gauss g; main(){ cin.tie(0); ios::sync_with_stdio(false); cin >> m[0] >> m[1]; g.solve(2, e, m, r); double ans; if(r[0] >= 0 && r[1] >= 0){ ans = r[0] * 1000 + r[1] * 2000; }else if(r[1] < 0){ ans = m[1] * 4 * 1000; }else{ ans = m[0] * 7 / 2 * 2000; } dout(ans); }