結果
問題 | No.1381 Simple Geometry 1 |
ユーザー |
![]() |
提出日時 | 2021-02-07 21:17:30 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,651 bytes |
コンパイル時間 | 2,357 ms |
コンパイル使用メモリ | 192,732 KB |
最終ジャッジ日時 | 2025-01-18 14:22:38 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 29 |
ソースコード
#include <bits/stdc++.h>using namespace std;using ll = long long;using ull = unsigned long long;#define v(t) vector<t>#define p(t) pair<t, t>#define p2(t, s) pair<t, s>#define vp(t) v(p(t))#define rep(i, n) for (int i=0,i##_len=((int)(n)); i<i##_len; ++i)#define rep2(i, a, n) for (int i=((int)(a)),i##_len=((int)(n)); i<=i##_len; ++i)#define repr(i, n) for (int i=((int)(n)-1); i>=0; --i)#define rep2r(i, a, n) for (int i=((int)(n)),i##_len=((int)(a)); i>=i##_len; --i)#define repi(itr, c) for (__typeof((c).begin()) itr=(c).begin(); itr!=(c).end(); ++itr)#define repir(itr, c) for (__typeof((c).rbegin()) itr=(c).rbegin(); itr!=(c).rend(); ++itr)#define sz(x) ((int)(x).size())#define all(x) (x).begin(), (x).end()#define rall(x) (x).rbegin(), (x).rend()#define SORT(v, n) sort(v, v+n);#define VSORT(v) sort(v.begin(), v.end());#define RSORT(x) sort(rall(x));#define pb push_back#define eb emplace_back#define INF (1e9)#define LINF (1e18)#define PI (acos(-1))#define EPS (1e-7)#define DEPS (1e-11)int main(){ll x, y, z, w;cin >> x >> y >> z >> w;double low = y, high = x;while (fabs(low-high)>DEPS) {double cent = (low+high) / 2;double cval = ((cent+w)*(cent+w) - (double)z*z) * (cent*cent - (double)y*y);if (cval >= (double)x*x) high = cent;else low = cent;}double weight = sqrt(low*low-(double)y*y), height = sqrt((low+w)*(low+w)-(double)z*z);double area = (double)x - ((double)z*height/2.0) - ((double)y*weight/2.0) - ((weight-z)*(height-y)/2.0);printf("%.10f\n", area);return 0;}