#303377 (cLay) No.764 浮動点

• 問題: No.764 浮動点
• ユーザー: tails
• 提出日時: 2018-12-12 01:15:08
• 言語: cLay (20240104-1)
• 結果: AC
• 実行時間: 3 ms / 1,500 ms
• コード長: 462 bytes
• コンパイル時間: 2,663 ms
• コンパイル使用メモリ: 173,984 KB
• 実行使用メモリ: 4,388 KB
• 最終ジャッジ日時: 2023-09-19 00:14:08
• 合計ジャッジ時間: 4,484 ms
• ジャッジサーバーID （参考情報）: judge13 / judge14

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1 ms; 4,380 KB
testcase_01	AC	1 ms; 4,380 KB
testcase_02	AC	1 ms; 4,376 KB
testcase_03	AC	1 ms; 4,380 KB
testcase_04	AC	3 ms; 4,380 KB
testcase_05	AC	2 ms; 4,380 KB
testcase_06	AC	3 ms; 4,388 KB
testcase_07	AC	2 ms; 4,380 KB
testcase_08	AC	2 ms; 4,380 KB
testcase_09	AC	2 ms; 4,380 KB
testcase_10	AC	2 ms; 4,380 KB
testcase_11	AC	2 ms; 4,380 KB
testcase_12	AC	3 ms; 4,380 KB
testcase_13	AC	3 ms; 4,384 KB
testcase_14	AC	3 ms; 4,376 KB
testcase_15	AC	3 ms; 4,380 KB
testcase_16	AC	3 ms; 4,384 KB
testcase_17	AC	3 ms; 4,380 KB
testcase_18	AC	2 ms; 4,380 KB
testcase_19	AC	2 ms; 4,376 KB
testcase_20	AC	3 ms; 4,388 KB
testcase_21	AC	2 ms; 4,384 KB
testcase_22	AC	3 ms; 4,380 KB
testcase_23	AC	3 ms; 4,376 KB
testcase_24	AC	3 ms; 4,384 KB
testcase_25	AC	1 ms; 4,384 KB

ソースコード

int n;
double r,l[2000],a,b,c,d;

double f(double a,double b){
	if(a<0)a=0;
	if(b<0)b=0;
	if(a+b<=r)return 0;
	if(b+r<=a)return b*b*PI;
	if(a+r<=b)return a*a*PI;
	double ta=acos((r*r+a*a-b*b)/(r*a*2));
	double tb=acos((r*r+b*b-a*a)/(r*b*2));
	return a*a*(ta-sin(2*ta)/2)+b*b*(tb-sin(2*tb)/2);
}

{
	rd(n,r,l(n+1));
	b=sum(l(n+1));
	rep(i,n){
		a+=l[i];
		b-=l[i];
		c=max[k=0---i](l[k])*2-a;
		d=max[k=i+1---n](l[k])*2-b;
		wt(f(a,b)-f(a,d)-f(c,b)+f(c,d));
	}
}