#881414 (C++17(gcc12)) No.2352 Sharpened Knife in Fall

• 問題: No.2352 Sharpened Knife in Fall
• ユーザー: 沙耶花
• 提出日時: 2023-06-16 21:29:44
• 言語: C++17(gcc12) (gcc 12.3.0 + boost 1.87.0)
• 結果: WA
• 実行時間: -
• コード長: 961 bytes
• コンパイル時間: 4,607 ms
• コンパイル使用メモリ: 266,288 KB
• 実行使用メモリ: 6,944 KB
• 最終ジャッジ日時: 2024-06-24 13:03:45
• 合計ジャッジ時間: 14,032 ms
• ジャッジサーバーID （参考情報）: judge3 / judge4

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 6,812 KB
testcase_01	AC	2 ms; 6,940 KB
testcase_02	AC	2 ms; 6,940 KB
testcase_03	AC	2 ms; 6,940 KB
testcase_04	AC	303 ms; 6,940 KB
testcase_05	AC	2 ms; 6,940 KB
testcase_06	AC	381 ms; 6,944 KB
testcase_07	AC	191 ms; 6,944 KB
testcase_08	AC	381 ms; 6,944 KB
testcase_09	WA	-
testcase_10	AC	382 ms; 6,940 KB
testcase_11	AC	379 ms; 6,944 KB
testcase_12	AC	378 ms; 6,940 KB
testcase_13	AC	376 ms; 6,944 KB
testcase_14	AC	179 ms; 6,940 KB
testcase_15	WA	-
testcase_16	AC	20 ms; 6,940 KB
testcase_17	AC	10 ms; 6,944 KB
testcase_18	AC	137 ms; 6,940 KB
testcase_19	WA	-
testcase_20	WA	-
testcase_21	WA	-

ソースコード

#include <stdio.h>
#include <atcoder/all>
#include <bits/stdc++.h>
using namespace std;
using namespace atcoder;
using mint = modint998244353;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf32 1000000001
#define Inf64 4000000000000000001

long long R,K;

double get(double h){
	double r = R;
	return r * r * acos(1.0 - h/r) - (r-h) * sqrt(2.0 * r * h - h*h);
}

int main(){
	double pi = acos(-1.0);
	
	cin>>R>>K;
	
	double S = R;
	S *= R;
	S *= pi;
	
	S /= K+1;
	
	vector<double> ans;
	double last = 0.0;
	rep(i,1000000){
		double ok = last,ng = 2.0*R;
		rep(_,50){
			double mid = (ok+ng)/2.0;
			if(get(mid)-get(last)<=S)ok = mid;
			else ng = mid;
		}
		if(ok >= R-1e-10)break;
		ans.push_back(-R + ok);
		last = ok;
	}
	int t = ans.size();
	rep(i,t)ans.push_back(-ans[i]);
	if(ans.size()!=K)ans.push_back(0.0);
	sort(ans.begin(),ans.end());
	
	rep(i,K){
		cout<<fixed<<setprecision(10)<<ans[i]<<endl;
	}
	
	//cout<<ans.size()<<endl;
	return 0;
}