結果

問題 No.2352 Sharpened Knife in Fall
ユーザー akiaa11akiaa11
提出日時 2023-06-27 15:41:57
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,967 ms / 3,000 ms
コード長 2,546 bytes
コンパイル時間 2,058 ms
コンパイル使用メモリ 202,340 KB
実行使用メモリ 4,692 KB
最終ジャッジ日時 2023-09-17 13:05:11
合計ジャッジ時間 32,980 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,384 KB
testcase_04 AC 1,904 ms
4,492 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 1,962 ms
4,376 KB
testcase_07 AC 983 ms
4,380 KB
testcase_08 AC 1,963 ms
4,500 KB
testcase_09 AC 1,962 ms
4,692 KB
testcase_10 AC 1,967 ms
4,384 KB
testcase_11 AC 1,962 ms
4,388 KB
testcase_12 AC 1,962 ms
4,444 KB
testcase_13 AC 1,959 ms
4,376 KB
testcase_14 AC 924 ms
4,380 KB
testcase_15 AC 1,849 ms
4,532 KB
testcase_16 AC 95 ms
4,380 KB
testcase_17 AC 44 ms
4,380 KB
testcase_18 AC 705 ms
4,376 KB
testcase_19 AC 1,540 ms
4,436 KB
testcase_20 AC 604 ms
4,384 KB
testcase_21 AC 577 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define SELECTER(_1,_2,_3,SELECT,...) SELECT
#define rep1(i,n) for(int i=0;i<(int)n;++i)
#define rep2(i,a,n) for(int i=(int)a;i<(int)n;++i)
#define rep(...) SELECTER(__VA_ARGS__,rep2,rep1)(__VA_ARGS__)
#define RSELECTER(_1, _2, _3, RSELECT, ...) RSELECT
#define rrep1(i,n) for(int i=(int)n-1;i>=0;--i)
#define rrep2(i,a,n) for(int i=(int)n-1;i>=(int)a;--i)
#define rrep(...) RSELECTER(__VA_ARGS__, rrep2, rrep1)(__VA_ARGS__)
#define all(a) a.begin(),a.end()
#define rall(a) a.rbegin(),a.rend()
#define fi first
#define se second
#define PrintR LogOutput
#ifdef _DEBUG
#define Log(...) LogOutput(__VA_ARGS__)
#else
#define Log(...)
#endif
using namespace std;
using ll=long long;
using pii=pair<int,int>;
using pll=pair<long long,long long>;
using pdd=pair<long double,long double>;
using veci=vector<int>;
using vecpii=vector<pair<int,int>>;
using vecll=vector<long long>;
using vecpll=vector<pair<long long,long long>>;
using vecpdd=vector<pair<long double,long double>>;
using vecs=vector<string>;
using vecb=vector<bool>;
using vecd=vector<long double>;
template<class T> inline bool chmax(T& a,T b){if(a<b) {a=b;return true;} return false;}
template<class T> inline bool chmin(T& a,T b){if(a>b) {a=b;return true;} return false;}

constexpr ll INF=numeric_limits<ll>::max() / 4;
constexpr ll MOD=998244353;
using tp=tuple<int,int,int>;
using tpll=tuple<ll,ll,ll>;

template<typename T>
T LogOutput(T a, bool b = true){
	if(b) cout << a << endl;
	return a;
}
template<typename... Args>
void LogOutput(Args&&... args){
	stringstream ss;
	((ss << LogOutput(args, false) << ' '), ...);
	cout << ss.str().substr(0, ss.str().length() - 1) << endl;
}
template<typename T>
void LogOutput(vector<T>& data){
	for(int i = 0 ; i < data.size() ; ++i){
		cout << LogOutput(data[i], false) << " \n"[i == data.size() - 1];
	}
}
template<typename T>
void LogOutput(vector<vector<T>>& data){
	for(auto d : data) LogOutput(d);
}

const double EPS = 1e-9;
int R, k;
double s;

bool f(double rad){
	if(double(R) * R * (rad - sin(rad)) / 2 <= s + EPS) return true;
	else return false;
}

int main(){
	ios::sync_with_stdio(false);
	std::cin.tie(nullptr);
	cin>>R>>k;
	s = (M_PI * R * R) / double(k + 1);
	double t = s;
	vecd ans;
	rep(i,  k / 2){
		double l = 0, r = M_PI;
		rep(j, 1000){
			double mid = (l + r) / 2;
			if(f(mid)) l = mid;
			else r = mid;
		}
		ans.push_back(R * cos(l / 2));
		s += t;
	}
	rep(i, ans.size()) cout<<fixed<<setprecision(6)<<-ans[i]<<endl;
	if(k % 2) cout<<0<<endl;
	rrep(i, ans.size()) cout<<ans[i]<<endl;
}
0