結果

問題 No.764 浮動点
コンテスト
ユーザー rickytheta
提出日時 2018-12-12 00:51:13
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 5 ms / 1,500 ms
コード長 1,350 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 886 ms
コンパイル使用メモリ 175,932 KB
実行使用メモリ 6,400 KB
最終ジャッジ日時 2026-04-11 20:47:46
合計ジャッジ時間 1,949 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;

#define REP(i,n) for(int i=0;i<(int)(n);i++)
#define FOR(i,a,b) for(int i=(int)(a);i<(int)(b);i++)
#define FORR(i,a,b) for(int i=(int)(b)-1;i>=(int)(a);i--)

int n;
double l[1252];

const double PI = 4.0*atan(1.0);

double calc(double r1, double r2, double d){
  if(r1+r2<d)return 0.0;
  if(d+r1<=r2)return r1*r1*PI;
  if(d+r2<=r1)return r2*r2*PI;
  double cos1 = (d*d + r1*r1 - r2*r2) / (2*d*r1);
  double cos2 = (d*d + r2*r2 - r1*r1) / (2*d*r2);
  double theta1 = acos(cos1);
  double theta2 = acos(cos2);
  return r1*r1*theta1 + r2*r2*theta2 - d*r1*sin(theta1);
}

int main(){
  scanf("%d",&n);
  REP(i,n+2)scanf("%lf",l+i);
  FOR(i,1,n+1){
    double maxlen1 = 0.0, minlen1 = 0.0;
    {
      double sm = 0.0, mx = -1e252;
      REP(j,i)sm += l[j+1];
      REP(j,i)mx = max(mx, l[j+1]);
      maxlen1 = sm;
      minlen1 = max(0.0, mx-(sm-mx));
    }
    double maxlen2 = 0.0, minlen2 = 0.0;
    {
      double sm = 0.0, mx = -1e252;
      FOR(j,i,n+2)sm += l[j+1];
      FOR(j,i,n+2)mx = max(mx, l[j+1]);
      maxlen2 = sm;
      minlen2 = max(0.0, mx-(sm-mx));
    }
    double ans = 0.0;
    ans += calc(maxlen1, maxlen2, l[0]);
    ans -= calc(minlen1, maxlen2, l[0]);
    ans -= calc(maxlen1, minlen2, l[0]);
    ans += calc(minlen1, minlen2, l[0]);
    printf("%.9f\n",ans);
  }
  return 0;
}
0