結果
問題 | No.764 浮動点 |
ユーザー | paruki |
提出日時 | 2019-01-22 19:46:29 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 8 ms / 1,500 ms |
コード長 | 2,052 bytes |
コンパイル時間 | 1,498 ms |
コンパイル使用メモリ | 172,180 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-16 04:17:30 |
合計ジャッジ時間 | 2,854 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 1 ms
5,376 KB |
testcase_03 | AC | 3 ms
5,376 KB |
testcase_04 | AC | 6 ms
5,376 KB |
testcase_05 | AC | 3 ms
5,376 KB |
testcase_06 | AC | 7 ms
5,376 KB |
testcase_07 | AC | 3 ms
5,376 KB |
testcase_08 | AC | 6 ms
5,376 KB |
testcase_09 | AC | 4 ms
5,376 KB |
testcase_10 | AC | 6 ms
5,376 KB |
testcase_11 | AC | 4 ms
5,376 KB |
testcase_12 | AC | 7 ms
5,376 KB |
testcase_13 | AC | 8 ms
5,376 KB |
testcase_14 | AC | 7 ms
5,376 KB |
testcase_15 | AC | 8 ms
5,376 KB |
testcase_16 | AC | 8 ms
5,376 KB |
testcase_17 | AC | 8 ms
5,376 KB |
testcase_18 | AC | 7 ms
5,376 KB |
testcase_19 | AC | 7 ms
5,376 KB |
testcase_20 | AC | 7 ms
5,376 KB |
testcase_21 | AC | 8 ms
5,376 KB |
testcase_22 | AC | 8 ms
5,376 KB |
testcase_23 | AC | 8 ms
5,376 KB |
testcase_24 | AC | 8 ms
5,376 KB |
testcase_25 | AC | 3 ms
5,376 KB |
ソースコード
#define _USE_MATH_DEFINES #include "bits/stdc++.h" using namespace std; #define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i)) #define rep(i,j) FOR(i,0,j) #define each(x,y) for(auto &(x):(y)) #define mp make_pair #define MT make_tuple #define all(x) (x).begin(),(x).end() #define debug(x) cout<<#x<<": "<<(x)<<endl #define smax(x,y) (x)=max((x),(y)) #define smin(x,y) (x)=min((x),(y)) #define MEM(x,y) memset((x),(y),sizeof (x)) #define sz(x) (int)(x).size() #define RT return using ll = long long; using pii = pair<int, int>; using vi = vector<int>; using vll = vector<ll>; int N, L0; vi L; pair<double, double> radius(int l, int r) { double sm = accumulate(L.begin() + l, L.begin() + r, 0.0); double ma = *max_element(L.begin() + l, L.begin() + r); return mp(sm, max(2 * ma - sm, 0.0)); } double area(double r1, double r2) { if (r1 + r2 <= L0)return 0; if (L0 + r2 <= r1)return r2 * r2 * M_PI; if (L0 + r1 <= r2)return r1 * r1 * M_PI; // 余弦定理から角度を求める double t1 = acos((r1 * r1 + L0 * L0 - r2 * r2) / (2 * r1 * L0)) * 2; double t2 = acos((r2 * r2 + L0 * L0 - r1 * r1) / (2 * r2 * L0)) * 2; // 弓形の面積 // 角度が180度未満の時 // (扇形の面積)ー(三角形の面積) // 角度が180度を越える時 // (扇形の面積)+(三角形の面積) // π<θ<2πのときsin(θ)は負なので // 面積も負になるためうまく計算できる。 double res = 0; res += 0.5*r1*r1*(t1 - sin(t1)); res += 0.5*r2*r2*(t2 - sin(t2)); return res; } void solve() { cin >> N >> L0; L.resize(N + 1); rep(i, N + 1)cin >> L[i]; for (int i = 1; i <= N; ++i) { double R1, r1, R2, r2; tie(R1, r1) = radius(0, i); tie(R2, r2) = radius(i, N + 1); double ans = area(R1, R2) - area(R1, r2) - area(r1, R2) + area(r1, r2); cout << ans << endl; } } int main() { ios::sync_with_stdio(false); cin.tie(0); cout << fixed << setprecision(15); solve(); return 0; }