結果

問題 No.1385 Simple Geometry 2
コンテスト
ユーザー mugen_1337
提出日時 2022-12-29 17:08:32
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,823 bytes
コンパイル時間 2,102 ms
コンパイル使用メモリ 195,064 KB
最終ジャッジ日時 2025-02-09 21:50:28
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 TLE * 55
権限があれば一括ダウンロードができます

ソースコード

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;

    double res = 0;
    for(int i = 0; i < N; i++){
        for(int j = i + 1; j < N; j++){
            {
                double theta = 2.0 * pi * (T[j] - T[i]) / L;
                res += 0.5 * sin(theta) * (N - j + i - 1);
                // cout << theta << " : + " << 0.5 * sin(theta) * (N - j + i - 1) << endl;
            }
            {
                double theta = 2.0 * pi * (T[j] - T[i]) / L;
                res -= 0.5 * sin(theta) * (j - i - 1);
                // cout << theta << " : - " << 0.5 * sin(theta) * (N - j + i - 1) << endl;
            }
        }
    }
    double den = double(N) * (N - 1) * (N - 2) / 6.0;
    cout << res / den << endl;
    return 0;
}
0