結果

問題 No.764 浮動点
コンテスト
ユーザー rickytheta
提出日時 2018-12-12 00:51:13
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 6 ms / 1,500 ms
コード長 1,350 bytes
コンパイル時間 1,256 ms
コンパイル使用メモリ 159,704 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-24 17:37:47
合計ジャッジ時間 2,230 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:25:8: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   25 |   scanf("%d",&n);
      |   ~~~~~^~~~~~~~~
main.cpp:26:18: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   26 |   REP(i,n+2)scanf("%lf",l+i);
      |             ~~~~~^~~~~~~~~~~

ソースコード

diff # 

#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