結果
問題 |
No.948 Bomb vs Dush
|
ユーザー |
|
提出日時 | 2019-12-08 16:50:01 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 6 ms / 2,000 ms |
コード長 | 746 bytes |
コンパイル時間 | 2,074 ms |
コンパイル使用メモリ | 165,708 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-12-30 18:56:49 |
合計ジャッジ時間 | 3,609 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 38 |
ソースコード
#include <bits/stdc++.h> using namespace std; constexpr long double PI = 3.1415926535897932384626433; //Euclidean distance template<class T> T EucDist(T x1, T y1, T x2, T y2) { return sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2)); } int main() { double r,R,ok = 180., ng = 360., md,x,y,d; cin >> r >> R; //10000回程度二分探索 for(int i = 0; i < 10000; ++i){ md = (ok+ng)/2.; x = R*cos(md/180.*PI); y = R*sin(md/180.*PI); d = EucDist(x,y,R,0.); //2円が共有範囲を持たない⇔(半径+半径)<(中心間の距離) を確かめる (2*r<d?ok:ng)=md; } //面積を求める double area = PI*((R+r)*(R+r)-(R-r)*(R-r))*ok/360. + PI*r*r; double t = PI*ok/180.; printf("%.10f\n",t); printf("%.10f\n",area); }