結果

問題 No.1042 愚直大学
コンテスト
ユーザー s0j1san
提出日時 2020-05-01 21:41:13
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 3,126 bytes
コンパイル時間 1,799 ms
コンパイル使用メモリ 167,460 KB
実行使用メモリ 6,816 KB
最終ジャッジ日時 2024-12-25 01:38:57
合計ジャッジ時間 2,709 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 21 WA * 2
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'void bisection(double, double, double, double*, int*)':
main.cpp:86:15: warning: 's' may be used uninitialized [-Wmaybe-uninitialized]
   86 |     *solution = s;
      |     ~~~~~~~~~~^~~
main.cpp:75:12: note: 's' was declared here
   75 |     double s;
      |            ^

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

//using Graph = vector<vector<edge>>;
using dou =long double;
string yes="yes";
string Yes="Yes";
string YES="YES";
string no="no";
string No="No";
string NO="NO";

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; }
#define rep(i, n)         for(int i = 0; i < (int)(n); i++)
#define brep(n)           for(int bit=0;bit<(1<<n);bit++)
#define erep(i,container) for (auto i : container)
#define irep(i, n)        for(int i = n-1; i >= (int)0; i--)
#define rrep(i,m,n) for(ll i = m; i < (ll)(n); i++)
#define reprep(i,j,h,w) rep(i,h)rep(j,w)
#define all(x) (x).begin(),(x).end()
#define VEC(type,name,n) std::vector<type> name(n);rep(i,n)std::cin >> name[i];
#define pb push_back
#define pf push_front
#define query int qq;std::cin >> qq;rep(qqq,qq)
#define lb lower_bound
#define ub upper_bound
#define fi first
#define se second
#define itn int
#define mp make_pair
#define sum(a) accumulate(all(a),0ll)
#define keta fixed<<setprecision
#define vvector(name,typ,m,n,a)vector<vector<typ> > name(m,vector<typ> (n,a))
#define vvvector(name,t,l,m,n,a) vector<vector<vector<t> > > name(l, vector<vector<t> >(m, vector<t>(n,a)));
#define vvvvector(name,t,k,l,m,n,a) vector<vector<vector<vector<t> > > > name(k,vector<vector<vector<t> > >(l, vector<vector<t> >(m, vector<t>(n,a)) ));
#define case std::cout <<"Case #" <<qqq+1<<": "
#define RES(a,i,j) a.resize(i);rep(ii,i)a[ii].resize(j); 
#define as assign
#define ffor for(;;)
#define ppri(a,b) std::cout << a<<" "<<b << std::endl
#define pppri(a,b,c) std::cout<<a<<" "<<b<<" "<<c<<std::endl
#define aall(x,n) (x).begin(),(x).begin()+(n)
#define ssum(a) accumulate(a,0ll) 
#define gre(n) greater<n>()

typedef long long ll;
typedef pair<int,int> P;
typedef pair<ll,ll> PL;
//typedef long long T;
const int INF = 1000000000;
const double dINF = 1000000000;
const ll INF64 = 922337203685477580ll;
const ll MOD = 1000000007ll;
const dou pi=3.141592653589793;
//メモ
//二次元配列のresizeのしかた
//vv.resize(n);
//for(size_t i=0; i<n; i++){
// vv[i].resize(n);
//}

//メモ
//負の辺を持つ最短経路問題では到達不可経路でも初期化値より小さくなりうるので注意しよう！
double p,q;
double f(double x)
{
    return x*x-p-q*x*log2(x);
}

// 2分法
void bisection(double a, double b, double eps, double *solution, int *N)
{
    int i = 0;
    double s;

    // 解が収束条件を満たせば終了
    while (!(abs(a-b)<eps)){
        i++;
        s = (a+b)/2.0;
        if(f(s) * f(a)<0) b=s;
        else a = s;
        if(i==1000) break; // 1000回繰り返したら強制終了
    };
    *N = i; 
    *solution = s; 
}



int main()
{   
    std::cin >> p>>q;
    double solution;
    int N;
    // 2分法
    bisection(0.0000001, 1.e50, 1.0e-10, &solution, &N);
    std::cout <<keta(10)<< solution << std::endl;
   // printf("解:%f (繰り返し回数：%d )",solution,N); // 解:1.165901 (繰り返し回数：18 )

    return 0;
}
0