結果

問題 No.1385 Simple Geometry 2
ユーザー cardano1016cardano1016
提出日時 2021-02-07 15:32:22
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 2,802 bytes
コンパイル時間 906 ms
コンパイル使用メモリ 108,096 KB
実行使用メモリ 19,552 KB
最終ジャッジ日時 2024-07-04 13:03:50
合計ジャッジ時間 33,781 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18 TLE * 47
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <iomanip>
#include <algorithm>
#include <assert.h>
#include <complex>
#include <utility>
#include <vector>
#include <string>
#include <stack>
#include <queue>
#include <tuple>
#include <cmath>
#include <bitset>
#include <cctype>
#include <set>
#include <map>
#include <unordered_map>
#include <numeric>
#include <functional>
#include <chrono>
#define _overload3(_1,_2,_3,name,...) name
#define _rep(i,n) repi(i,0,n)
#define repi(i,a,b) for(ll i=a;i<b;++i)
#define rep(...) _overload3(__VA_ARGS__,repi,_rep,)(__VA_ARGS__)
#define _rrep(i,a) rrepi(i,a,0)
#define rrepi(i,a,b) for(ll i=a-1;i>=b;--i)
#define rrep(...) _overload3(__VA_ARGS__,rrepi,_rrep,)(__VA_ARGS__)
#define all(x) (x).begin(),(x).end()
#define PRINT(V) cout << V << "\n"
#define SORT(V) sort((V).begin(),(V).end())
#define RSORT(V) sort((V).rbegin(), (V).rend())
using namespace std;
using ll = long long;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }
inline void Yes(bool condition){ if(condition) PRINT("Yes"); else PRINT("No"); }
template<class itr> void cins(itr first,itr last){
    for (auto i = first;i != last;i++){
        cin >> (*i);
    }
}
template<class itr> void array_output(itr start,itr goal){
    string ans = "",k = " ";
    for (auto i = start;i != goal;i++) ans += to_string(*i)+k;
    if (!ans.empty()) ans.pop_back();
    PRINT(ans);
}
ll gcd(ll a, ll b) {
    return a ? gcd(b%a,a) : b;
}

const ll INF = 1e18;
const ll MOD = 1000000007;
const ll MOD2 = 998244353;
const ll MOD3 = 1e6;
const ll EPS = 1e-10;
int sgn(const double a){
    return (a < -EPS ? -1 : (a > EPS ? +1 : 0));
}
typedef pair<int,int> pi;
typedef pair<ll,ll> P;
typedef tuple<ll,ll,ll> tri;
typedef pair<double,double> point;
typedef complex<double> Point;
typedef string::const_iterator State;
const ll MAX = 105;
constexpr ll nx[4] = {-1,0,1,0};
constexpr ll ny[4] = {0,1,0,-1};

const double PI = 3.1415926535;

int main(){
    ll n,l;
    cin >> n >> l;
    vector<long double> coss(n),sins(n);
    rep(i,n){
        long double t;
        cin >> t;
        coss[i] = cos(2*PI*t/(long double)l);
        sins[i] = sin(2*PI*t/(long double)l);
    }
    long double sumsin = 0.0,sumcos = 0.0;
    rep(i,n){
        sumsin += sins[i];
        sumcos += coss[i];
    }
    long double ans = 0.0,nowsin = 0.0,nowcos = 0.0;
    rep(i,n-1){
        nowsin += (n-i-2)*sins[i+1];
        nowcos += (n-i-2)*coss[i+1];
    }
    rep(i,n){
        ans += nowsin*coss[i]-nowcos*sins[i];
        nowsin = nowsin+sumsin-(n-1)*sins[(i+1)%n]-sins[i];
        nowcos = nowcos+sumcos-(n-1)*coss[(i+1)%n]-coss[i];
    }
    ll k = n*(n-1)*(n-2)/6;
    printf("%.12Lf\n",ans/(long double)k/2.0);
}
0