結果

問題 No.764 浮動点
ユーザー omochanal
提出日時 2018-12-20 21:17:43
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 56 ms / 1,500 ms
コード長 3,140 bytes
コンパイル時間 1,678 ms
コンパイル使用メモリ 165,928 KB
実行使用メモリ 6,940 KB
最終ジャッジ日時 2024-09-25 09:11:35
合計ジャッジ時間 3,333 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define ADD(a, b) a = (a + ll(b)) % mod
#define MUL(a, b) a = (a * ll(b)) % mod
#define MAX(a, b) a = max(a, b)
#define MIN(a, b) a = min(a, b)
#define rep(i, a, b) for(int i = int(a); i < int(b); i++)
#define rer(i, a, b) for(int i = int(a) - 1; i >= int(b); i--)
#define all(a) (a).begin(), (a).end()
#define sz(v) (int)(v).size()
#define pb push_back
#define sec second
#define fst first
#define debug(fmt, ...) Debug(__LINE__, ":", fmt, ##__VA_ARGS__)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pi;
typedef pair<ll, ll> pl;
typedef pair<int, pi> ppi;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef vector<vl> mat;
typedef complex<double> comp;
void Debug() {cout << '\n'; }
template<class FIRST, class... REST>void Debug(FIRST arg, REST... rest){
	cout<<arg<<" ";Debug(rest...);}
template<class T>ostream& operator<<(ostream& out,const vector<T>& v) {
	out<<"[";if(!v.empty()){rep(i,0,sz(v)-1)out<<v[i]<<", ";out<<v.back();}out<<"]";return out;}
template<class S, class T>ostream& operator<<(ostream& out,const pair<S, T>& v){
	out<<"("<<v.first<<", "<<v.second<<")";return out;}
const int MAX_N = 300010;
const int MAX_V = 100010;
const double eps = 1e-6;
const ll mod = 1000000007;
const int inf = 1 << 30;
const ll linf = 1LL << 60;
const double PI = 3.14159265358979323846;
///////////////////////////////////////////////////////////////////////////////////////////////////

int N;
double R;
double L[MAX_N];

double Spre(double r1, double r2) {
	if(r1 + r2 <= R) return 0.0;
	else if(r1 + R <= r2) {
		return 1.0 / 2.0 * r1 * r1 * (2 * PI);
	}
	else if(r2 + R <= r1) {
		return 1.0 / 2.0 * r2 * r2 * (2 * PI);
	}
	else {
		double the1 = acos((-r2 * r2 + r1 * r1 + R * R) / (2 * r1 * R));
		double the2 = acos((-r1 * r1 + r2 * r2 + R * R) / (2 * r2 * R));
		// debug(r1, r2, the1, the2);
		return 1.0 / 2.0 * r1 * r1 * (2 * the1)
			+ 1.0 / 2.0 * r2 * r2 * (2 * the2)
			- 1.0 / 2.0 * r1 * r1 * sin(2 * the1)
			- 1.0 / 2.0 * r2 * r2 * sin(2 * the2);
	}
}

double S(double o1, double o2, double p1, double p2) {
	return Spre(o1, p1) + Spre(o2, p2) - Spre(o2, p1) - Spre(o1, p2);
}

void solve() {
	cin >> N;
	cin >> R;
	rep(i, 0, N + 1) cin >> L[i];
	rep(i, 1, N + 1) {
		vector<double> vec1;
		vector<double> vec2;
		rep(j, 0, i) vec1.pb(L[j]);
		rep(j, i, N + 1) vec2.pb(L[j]);
		sort(all(vec1));
		sort(all(vec2));
		int n1 = sz(vec1), n2 = sz(vec2);
		double o1 = accumulate(all(vec1), 0.0);
		double o2 = max(0.0, vec1.back() - accumulate(vec1.begin(), vec1.begin() + n1 - 1, 0.0));
		double p1 = accumulate(all(vec2), 0.0);
		double p2 = max(0.0, vec2.back() - accumulate(vec2.begin(), vec2.begin() + n2 - 1, 0.0));
		cout << S(o1, o2, p1, p2) << "\n";
	}
}

int main() {
#ifndef LOCAL
	ios::sync_with_stdio(false);
    cin.tie(0);
#endif
    cout << fixed;
	cout.precision(20);
	srand((unsigned int)time(NULL));
#ifdef LOCAL
	//freopen("in.txt", "wt", stdout); //for tester
    freopen("in.txt", "rt", stdin);
#endif	
	solve();
#ifdef LOCAL
    cerr << "Time elapsed: " << 1.0 * clock() / CLOCKS_PER_SEC << " s.\n";
#endif
	return 0;
}
0