結果

問題 No.2352 Sharpened Knife in Fall
コンテスト
ユーザー ueta
提出日時 2023-06-17 00:41:56
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 702 ms / 3,000 ms
コード長 4,319 bytes
コンパイル時間 1,931 ms
コンパイル使用メモリ 200,908 KB
最終ジャッジ日時 2025-02-14 21:46:20
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <cstdint>
#define rp(n) for(int i=0;i<n;i++)
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define Rep(i,sta,n) for(int i=sta;i<n;i++)
#define Per(i,sta,n) for(int i=n-1;i>=sta;i--)
#define fore(i, a) for(auto &&i : a)
#define foreP(a, b, v) for(auto &&[a, b] : v)
#define foreP3(a, b, c, v) for(auto &&[a, b, c] : v)
#define fitr(itr, m) for(auto &&itr=m.begin();itr!=m.end();itr++)
#define all(v) (v).begin(),(v).end()
#define Sum(a) accumulate(all(a),0LL)
#define VC vector
#define TP tuple

using namespace std;
mt19937 engine(time(0));

using ll = long long;
using ld = long double;
constexpr long long mod = 998244353;
// constexpr long long mod = 1000000007;
// constexpr long long mod = 100000000000000003;
constexpr long double PI=3.1415926535897932384626433832795028841971;
constexpr long long inf = INT_MAX;
constexpr long long infll = LONG_MAX;

template<typename T>using MaxHeap = priority_queue<T>;
template<typename T>using MinHeap = priority_queue<T, vector<T>, greater<T>>;
template<typename T>void chmin(T& a, T b) {a = min(a, b);}
template<typename T>void chmax(T& a, T b) {a = max(a, b);}

template<typename T>void debug(T a){cerr << a <<  endl;}
template<typename T,typename S>void debug(T a, S b){cerr << a << " " << b << endl;}
template<typename T,typename S,typename U>void debug(T a,S b,U c){cerr <<"( "<< a << " " << b << " " << c <<" )" << endl;}
template<typename T,typename S,typename U,typename W>void debug(T a,S b,U c,W d){cerr <<"( "<< a << " " << b << " " << c<< " " << d <<" )" << endl;}
template<typename T>void debug(vector<T> v){for(int i=0;i<v.size();i++)cerr << v[i] << " ";cerr << endl;}
template<typename T>void debug(set<T> m){for(auto itr: m){cerr << itr << " ";}cerr << endl;}
template<typename T,typename S>void debug(map<T,S> m){for(auto itr=m.begin();itr!=m.end();itr++)cerr <<"("<< itr->first << " => " << itr->second<<")" << endl;}
template<typename T,typename S>void debug(pair<T,S> p){cerr <<"( "<< p.first << ", " << p.second <<" )" << endl;}
template<typename T,typename S>void debug(tuple<T,S> t){cerr <<"( "<< get<0>(t) << ", " << get<1>(t) <<" )" << endl;}
template<typename T,typename S,typename U>void debug(tuple<T,S,U> t){cerr <<"( "<< get<0>(t) << ", " << get<1>(t)<< ", " << get<2>(t) <<" )" << endl;}
template<typename T,typename S,typename U,typename W>void debug(tuple<T,S,U,W> t){cerr <<"( "<< get<0>(t) << ", " << get<1>(t)<< ", " << get<2>(t)<< ", " << get<3>(t) <<" )" << endl;}

template<typename T>void out(T a){ cout << a << endl;}
template<typename T,typename S>void out(T a,S b){ cout << a <<" " << b << endl;}
template<typename T,typename S,typename U>void out(T a,S b,U c){ cout << a <<" " << b <<" " << c << endl;}
template<typename T>void out(vector<T> v){ for(int i=0;i<v.size();i++){cout << v[i];if(i==v.size()-1) cout << endl;else cout << " ";}}
template<typename T>void out(set<T> m){for(auto itr: m){cout << itr; if((itr+1)==m.endl()) cout << endl;else cout << " ";}}
template<typename T,typename S>void out(pair<T,S> p){cout << p.first << " " << p.second << endl;}
template<typename T,typename S>void out(tuple<T,S> t){cout << get<0>(t) << " " << get<1>(t) << endl;}
template<typename T,typename S,typename U>void out(tuple<T,S,U> t){cout << get<0>(t) << " " << get<1>(t)<< " " << get<2>(t) << endl;}
void yes(bool f=1,string yes="Yes",string no="No"){if(f){cout << yes << endl;}else{cout << no << endl;}}


struct  S
{
	int x;
	int y;
	S(int x, int y) : x(x),y(y){}
	bool operator<(const struct S& other) const
	{
		return x < other.x;
	}
	void debug(){
		cerr <<"debug " << x <<"-" << y << endl;
	}
};

void solve(){
	ld R;
	int k;
	cin >> R >> k;
	VC<ld> p(k,0.0);
	rep(i,k){
		ld s = PI/2.0*(i+1)/(k+1);
		ld l = -1,r = 1;
		while (r-l>0.0000000001)
		{
			ld c = (r-l)/2+l;
			// debug(c);
			ld sita=acos(c);
			ld ss;
				if(c<0)
					ss = PI-sita/2.0-abs(c)*sqrt(1.0-c*c)/2.0;
				else
					ss = sita/2.0-abs(c)*sqrt(1.0-c*c)/2.0;
			if(s<=ss){
				l = c;
			}else{
				r = c;
			}
		}
		p[i] = l*R;
	}

	rep(i,k/2){
		p[i+k/2] = -(ld)1*p[i];
	}
	sort(all(p));
	rep(i,k)out(p[i]);
	return;
}

int main()
{
	cin.tie(nullptr);
    ios::sync_with_stdio(false);
    cout<<fixed << std::setprecision(30);
	int t=1;
	// cin >> t;
	while(t--)
		solve();
	return 0;
}
0