結果
| 問題 |
No.1385 Simple Geometry 2
|
| コンテスト | |
| ユーザー |
蜜蜂
|
| 提出日時 | 2021-02-02 22:09:41 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 278 ms / 500 ms |
| コード長 | 1,243 bytes |
| コンパイル時間 | 2,087 ms |
| コンパイル使用メモリ | 167,972 KB |
| 実行使用メモリ | 23,544 KB |
| 最終ジャッジ日時 | 2024-07-04 12:58:33 |
| 合計ジャッジ時間 | 21,584 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 65 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ld = long double;
#define fi first
#define se second
#define pb push_back
const double pi=3.141592653589793;
ll inputinf=1e12;
int main(){
int n;
ll l;
cin>>n>>l;
if(n>500000){
cout<<"WARNING : N"<<endl;
}
if(l>inputinf){
cout<<"WARNING : L"<<endl;
}
ll t[n];
double tsin[n],tcos[n],tsinsum[n+1],tcossum[n+1],p;
tsinsum[0]=0,tcossum[0]=0;
for(int i=0;i<n;i++){
cin>>t[i];
if(t[i]>=l||t[i]<0){
cout<<"WARNING : t["<<i<<"]"<<endl;
}
if(i>1){
if(t[i-1]>=t[i]){
cout<<"WARNING : t["<<i<<"]"<<endl;
}
}
p=2*t[i];
p*=pi;
p/=l;
tsin[i]=sin(p),tcos[i]=cos(p);
tsinsum[i+1]=tsinsum[i]+tsin[i];
tcossum[i+1]=tcossum[i]+tcos[i];
}
double sinsum=0,cossum=0;
for(int i=1;i<n;i++){
sinsum+=tsin[i]*(i-1);
cossum+=tcos[i]*(i-1);
}
double ans=0;
for(int i=0;i<n;i++){
ans+=tsin[i]*cossum-tcos[i]*sinsum;
sinsum-=tsinsum[n],cossum-=tcossum[n];
sinsum+=tsin[i]*(n-1)+tsin[(i+1)%n],cossum+=tcos[i]*(n-1)+tcos[(i+1)%n];
}
ans*=3;
ans/=n;
ans/=n-1;
ans/=n-2;
cout<<fixed<<setprecision(15);
cout<<ans<<endl;
}
蜜蜂