結果

問題 No.764 浮動点
ユーザー chocoruskchocorusk
提出日時 2018-12-12 00:30:13
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 3 ms / 1,500 ms
コード長 1,458 bytes
コンパイル時間 871 ms
コンパイル使用メモリ 96,868 KB
実行使用メモリ 4,372 KB
最終ジャッジ日時 2023-10-24 22:29:26
合計ジャッジ時間 1,789 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,372 KB
testcase_01 AC 1 ms
4,372 KB
testcase_02 AC 1 ms
4,372 KB
testcase_03 AC 2 ms
4,372 KB
testcase_04 AC 3 ms
4,372 KB
testcase_05 AC 2 ms
4,372 KB
testcase_06 AC 3 ms
4,372 KB
testcase_07 AC 2 ms
4,372 KB
testcase_08 AC 3 ms
4,372 KB
testcase_09 AC 2 ms
4,372 KB
testcase_10 AC 3 ms
4,372 KB
testcase_11 AC 3 ms
4,372 KB
testcase_12 AC 3 ms
4,372 KB
testcase_13 AC 3 ms
4,372 KB
testcase_14 AC 3 ms
4,372 KB
testcase_15 AC 3 ms
4,372 KB
testcase_16 AC 3 ms
4,372 KB
testcase_17 AC 3 ms
4,372 KB
testcase_18 AC 2 ms
4,372 KB
testcase_19 AC 2 ms
4,372 KB
testcase_20 AC 3 ms
4,372 KB
testcase_21 AC 2 ms
4,372 KB
testcase_22 AC 3 ms
4,372 KB
testcase_23 AC 3 ms
4,372 KB
testcase_24 AC 3 ms
4,372 KB
testcase_25 AC 2 ms
4,372 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
using namespace std;
typedef long long int ll;
//typedef pair<int, int> P;
const double PI=acos(-1.0);
int n;
double l[1002];
double arg(double a, double b, double c){
	return acos((a*a+b*b-c*c)/(2*a*b));
}
double triarea(double a, double b, double c){
	double t=arg(a, b, c);
	return a*b*sin(t)/2;
}
double intersection(double r1, double r2){
	if(r1+r2<=l[0]) return 0;
	if(r2-r1>=l[0]) return r1*r1*PI;
	if(r1-r2>=l[0]) return r2*r2*PI;
	double s1=r1*r1*PI, s2=r2*r2*PI;
	double t1=arg(r1, l[0], r2), t2=arg(r2, l[0], r1);
	double s=s1*(PI-t1)/PI+s2*(PI-t2)/PI+2*triarea(r1, r2, l[0]);
	return s1+s2-s;
}
int main()
{
	cin>>n;
	for(int i=0; i<=n+1; i++) cin>>l[i];
	double mx1[1002], s[1002], mx2[1002]; mx1[1]=l[1], s[0]=0, s[1]=l[1];
	for(int i=2; i<=n+1; i++){
		mx1[i]=max(mx1[i-1], l[i]);
		s[i]=s[i-1]+l[i];
	}
	mx2[n]=l[n+1];
	for(int i=n-1; i>=1; i--){
		mx2[i]=max(mx2[i+1], l[i+1]);
	}
	for(int i=1; i<=n; i++){
		double s1=s[i], t1=max(0.0, 2*mx1[i]-s1);
		double s2=s[n+1]-s[i], t2=max(0.0, 2*mx2[i]-s2);
		double ans=intersection(s1, s2)+intersection(t1, t2)-intersection(s1, t2)-intersection(s2, t1);
		printf("%.7lf\n", ans);
	}
	return 0;
}
0