結果

問題 No.1385 Simple Geometry 2
ユーザー cardano1016cardano1016
提出日時 2021-02-07 15:33:20
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 346 ms / 500 ms
コード長 2,767 bytes
コンパイル時間 966 ms
コンパイル使用メモリ 108,016 KB
実行使用メモリ 11,728 KB
最終ジャッジ日時 2024-07-04 13:02:46
合計ジャッジ時間 24,204 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 5 ms
6,940 KB
testcase_09 AC 5 ms
6,940 KB
testcase_10 AC 5 ms
6,940 KB
testcase_11 AC 5 ms
6,944 KB
testcase_12 AC 5 ms
6,944 KB
testcase_13 AC 342 ms
11,516 KB
testcase_14 AC 336 ms
11,656 KB
testcase_15 AC 335 ms
11,548 KB
testcase_16 AC 340 ms
11,452 KB
testcase_17 AC 335 ms
11,440 KB
testcase_18 AC 335 ms
11,516 KB
testcase_19 AC 325 ms
11,504 KB
testcase_20 AC 337 ms
11,572 KB
testcase_21 AC 339 ms
11,556 KB
testcase_22 AC 316 ms
11,448 KB
testcase_23 AC 340 ms
11,608 KB
testcase_24 AC 321 ms
11,556 KB
testcase_25 AC 318 ms
11,524 KB
testcase_26 AC 338 ms
11,660 KB
testcase_27 AC 344 ms
11,652 KB
testcase_28 AC 321 ms
11,520 KB
testcase_29 AC 338 ms
11,612 KB
testcase_30 AC 337 ms
11,528 KB
testcase_31 AC 337 ms
11,608 KB
testcase_32 AC 332 ms
11,476 KB
testcase_33 AC 336 ms
11,520 KB
testcase_34 AC 331 ms
11,488 KB
testcase_35 AC 338 ms
11,516 KB
testcase_36 AC 346 ms
11,488 KB
testcase_37 AC 344 ms
11,444 KB
testcase_38 AC 337 ms
11,632 KB
testcase_39 AC 338 ms
11,512 KB
testcase_40 AC 339 ms
11,548 KB
testcase_41 AC 338 ms
11,664 KB
testcase_42 AC 331 ms
11,484 KB
testcase_43 AC 340 ms
11,692 KB
testcase_44 AC 339 ms
11,728 KB
testcase_45 AC 341 ms
11,504 KB
testcase_46 AC 337 ms
11,572 KB
testcase_47 AC 340 ms
11,516 KB
testcase_48 AC 342 ms
11,668 KB
testcase_49 AC 340 ms
11,452 KB
testcase_50 AC 339 ms
11,444 KB
testcase_51 AC 340 ms
11,624 KB
testcase_52 AC 339 ms
11,660 KB
testcase_53 AC 341 ms
11,568 KB
testcase_54 AC 342 ms
11,488 KB
testcase_55 AC 342 ms
11,512 KB
testcase_56 AC 343 ms
11,444 KB
testcase_57 AC 341 ms
11,432 KB
testcase_58 AC 341 ms
11,568 KB
testcase_59 AC 292 ms
11,696 KB
testcase_60 AC 322 ms
11,500 KB
testcase_61 AC 323 ms
11,604 KB
testcase_62 AC 323 ms
11,488 KB
testcase_63 AC 343 ms
11,524 KB
testcase_64 AC 342 ms
11,488 KB
testcase_65 AC 343 ms
11,632 KB
testcase_66 AC 340 ms
11,444 KB
testcase_67 AC 334 ms
11,724 KB
権限があれば一括ダウンロードができます

ソースコード

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<double> coss(n),sins(n);
    rep(i,n){
        double t;
        cin >> t;
        coss[i] = cos(2*PI*t/(double)l);
        sins[i] = sin(2*PI*t/(double)l);
    }
    double sumsin = 0.0,sumcos = 0.0;
    rep(i,n){
        sumsin += sins[i];
        sumcos += coss[i];
    }
    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("%.10lf\n",ans/(double)k/2.0);
}
0