結果

問題 No.1385 Simple Geometry 2
ユーザー mugen_1337
提出日時 2022-12-29 18:47:54
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 143 ms / 500 ms
コード長 2,086 bytes
コンパイル時間 1,929 ms
コンパイル使用メモリ 198,416 KB
最終ジャッジ日時 2025-02-09 21:56:55
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 65
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define ALL(x) begin(x), end(x)
#define rep(i, n) for (int i = 0; i < (n); i++)
#define mod 1000000007
using ll = long long;
const int INF = 1000000000;
const ll LINF = 1001002003004005006ll;
int dx[] = {1, 0, -1, 0}, dy[] = {0, 1, 0, -1};
// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template <class T>
bool chmax(T &a, const T &b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}
template <class T>
bool chmin(T &a, const T &b) {
    if (b < a) {
        a = b;
        return true;
    }
    return false;
}

struct IOSetup {
    IOSetup() {
        cin.tie(0);
        ios::sync_with_stdio(0);
        cout << fixed << setprecision(12);
    }
} iosetup;

template <typename T>
ostream &operator<<(ostream &os, const vector<T> &v) {
    for (int i = 0; i < (int)v.size(); i++) os << v[i] << (i + 1 == (int)v.size() ? "" : " ");
    return os;
}
template <typename T>
istream &operator>>(istream &is, vector<T> &v) {
    for (T &x : v) is >> x;
    return is;
}
double pi = acos(-1);
signed main() {
    ll N, L; cin >> N >> L;
    vector<ll> T(N); cin >> T;
    for(int i = 0; i < N; i++) T.push_back(T[i] + L);

    vector<double> sins(T.size() + 1), coss(T.size() + 1);
    sins[0] = 0.0, coss[0] = 0.0;
    rep(i, (int)T.size()){
        sins[i + 1] = sins[i] + sin(2.0 * pi * T[i] / L);
        coss[i + 1] = coss[i] + cos(2.0 * pi * T[i] / L);
    }

    double cos_ = 0, sin_ = 0;
    for(int i = 1; i < N; i++){
        sin_ += sin(2.0 * pi * T[i] / L) * (N - 1 - i);
        cos_ += cos(2.0 * pi * T[i] / L) * (N - 1 - i);
    }

    double res = 0;
    for(int i = 0; i < N; i++){
        res += 0.5 * cos(2.0 * pi * T[i] / L) * sin_;
        res -= 0.5 * sin(2.0 * pi * T[i] / L) * cos_;

        sin_ -= sin(2.0 * pi * T[i + 1] / L) * (N - 2);
        sin_ += sins[N + i] - sins[i + 2];

        cos_ -= cos(2.0 * pi * T[i + 1] / L) * (N - 2);
        cos_ += coss[N + i] - coss[i + 2];
    }
    double den = double(N) * (N - 1) * (N - 2) / 6.0;
    cout << res / den << endl;
    return 0;
}
0