結果
問題 | No.2352 Sharpened Knife in Fall |
ユーザー |
|
提出日時 | 2023-06-16 22:31:14 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 377 ms / 3,000 ms |
コード長 | 1,788 bytes |
コンパイル時間 | 2,162 ms |
コンパイル使用メモリ | 202,548 KB |
最終ジャッジ日時 | 2025-02-14 06:29:54 |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 19 |
ソースコード
#include<bits/stdc++.h> #define fi first #define se second #define eb emplace_back #define mp make_pair using namespace std; typedef long double ld; typedef long long ll; typedef unsigned long long ull; typedef __int128 i128; template<typename T,typename U> T ceil(T x, U y) {return (x>0?(x+y-1)/y:x/y);} template<typename T,typename U> T floor(T x, U y) {return (x>0?x/y:(x-y+1)/y);} template<class T,class S> bool chmax(T &a,const S b) {return (a<b?a=b,1:0);} template<class T,class S> bool chmin(T &a,const S b) {return (a>b?a=b,1:0);} int popcnt(int x) {return __builtin_popcount(x);} int popcnt(ll x) {return __builtin_popcountll(x);} int topbit(int x) {return (x==0?-1:31-__builtin_clz(x));} int topbit(ll x) {return (x==0?-1:63-__builtin_clzll(x));} int lowbit(int x) {return (x==0?-1:__builtin_ctz(x));} int lowbit(ll x) {return (x==0?-1:__builtin_ctzll(x));} #define rep(i,a,b) for(int i=(a);i<=(b);i++) #define per(i,a,b) for(int i=(a);i>=(b);i--) typedef pair<int,int> pii; typedef vector<int> vi; typedef vector<pii> vp; int read() { int x=0,w=1; char c=getchar(); while(!isdigit(c)) {if(c=='-') w=-1; c=getchar();} while(isdigit(c)) {x=x*10+(c-'0'); c=getchar();} return x*w; } const ld eps=1e-10; int k; ld r,g,s,pi; vector<ld> ans; ld calc(ld h) { ld ang=acos(h/r); ld t=sqrt(r*r-h*h)*h; return s/2-s*ang/pi+t; } signed main() { r=read(), k=read(), pi=acos(-1), s=pi*r*r, g=s/(k+1); rep(i,1,k/2) { ld x=0, y=r; ld p=i; if(k%2==0) p-=0.5; while(y-x>eps) { ld mid=(x+y)/2; if(calc(mid)>g*p) y=mid; else x=mid; } ans.eb(x); } sort(ans.begin(),ans.end(),greater<ld>()); for(ld x:ans) printf("-%.10Lf\n",x); if(k&1) puts("0"); sort(ans.begin(),ans.end()); for(ld x:ans) printf("%.10Lf\n",x); return 0; }