結果

問題 No.1042 愚直大学
コンテスト
ユーザー tunny
提出日時 2020-05-01 23:13:55
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 636 ms / 2,000 ms
コード長 921 bytes
コンパイル時間 1,808 ms
コンパイル使用メモリ 165,672 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-25 14:18:47
合計ジャッジ時間 8,522 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=int(a);i<=int(b);++i)
using namespace std;
template<class A>void pr(A a){cout << (a) << endl;}
template<class A,class B>void pr(A a,B b){cout << a << " "  ;pr(b);}
template<class A,class B,class C>void pr(A a,B b,C c){cout << a << " " ;pr(b,c);}
template<class A,class B,class C,class D>void pr(A a,B b,C c,D d){cout << a << " " ;pr(b,c,d);}
typedef long long ll;
typedef pair<ll, ll> l_l;
typedef pair<int, int> i_i;
int INF=numeric_limits<int>::max();
// for(int j=0;j<N;++j) {
double func(double x, ll p, ll q) {
  return x*x-q*x*log2(x)-p;
}

double dfunc(double x, ll q) {
  return 2*x-q*(log2(x)+1.0/log(2))-1;
}

int main(void) {
  ll p,q;cin>> p >> q;
  double x = 100*q + p;
  int cnt=10000000;
  while(abs(func(x,p,q))>0.000001) {
    x = x - func(x,p,q) / dfunc(x,q);
    cnt--;
    if(cnt<0) {
      break;
    }
  }
  printf("%8.8f\n", x);
  return 0;
}
0