結果

問題 No.2352 Sharpened Knife in Fall
ユーザー lgswdnlgswdn
提出日時 2023-06-16 22:31:14
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 359 ms / 3,000 ms
コード長 1,788 bytes
コンパイル時間 2,107 ms
コンパイル使用メモリ 206,632 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-09-06 20:47:28
合計ジャッジ時間 12,026 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 195 ms
4,380 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 331 ms
4,380 KB
testcase_07 AC 163 ms
4,380 KB
testcase_08 AC 359 ms
4,384 KB
testcase_09 AC 337 ms
4,384 KB
testcase_10 AC 330 ms
4,384 KB
testcase_11 AC 331 ms
4,380 KB
testcase_12 AC 345 ms
4,380 KB
testcase_13 AC 338 ms
4,384 KB
testcase_14 AC 145 ms
4,384 KB
testcase_15 AC 296 ms
4,380 KB
testcase_16 AC 17 ms
4,384 KB
testcase_17 AC 10 ms
4,384 KB
testcase_18 AC 133 ms
4,384 KB
testcase_19 AC 266 ms
4,384 KB
testcase_20 AC 102 ms
4,384 KB
testcase_21 AC 99 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0