#644608 (C++17) No.955 ax^2+bx+c=0

提出ソース

結果

問題	No.955 ax^2+bx+c=0
ユーザー	tanimani364
提出日時	2021-04-04 19:39:39
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	3 ms / 2,000 ms
コード長	2,537 bytes
コンパイル時間	13,748 ms
コンパイル使用メモリ	376,476 KB
実行使用メモリ	4,384 KB
最終ジャッジ日時	2023-08-27 21:00:40
合計ジャッジ時間	12,537 ms
ジャッジサーバーID （参考情報）	judge12 / judge14

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 4,380 KB
testcase_01	AC	1 ms 4,380 KB
testcase_02	AC	2 ms 4,376 KB
testcase_03	AC	2 ms 4,376 KB
testcase_04	AC	1 ms 4,384 KB
testcase_05	AC	2 ms 4,376 KB
testcase_06	AC	2 ms 4,380 KB
testcase_07	AC	1 ms 4,376 KB
testcase_08	AC	1 ms 4,380 KB
testcase_09	AC	1 ms 4,376 KB
testcase_10	AC	2 ms 4,380 KB
testcase_11	AC	1 ms 4,380 KB
testcase_12	AC	2 ms 4,384 KB
testcase_13	AC	1 ms 4,376 KB
testcase_14	AC	1 ms 4,380 KB
testcase_15	AC	2 ms 4,380 KB
testcase_16	AC	1 ms 4,380 KB
testcase_17	AC	2 ms 4,380 KB
testcase_18	AC	1 ms 4,376 KB
testcase_19	AC	1 ms 4,380 KB
testcase_20	AC	2 ms 4,384 KB
testcase_21	AC	2 ms 4,380 KB
testcase_22	AC	2 ms 4,380 KB
testcase_23	AC	2 ms 4,380 KB
testcase_24	AC	2 ms 4,380 KB
testcase_25	AC	1 ms 4,380 KB
testcase_26	AC	2 ms 4,376 KB
testcase_27	AC	2 ms 4,380 KB
testcase_28	AC	2 ms 4,376 KB
testcase_29	AC	1 ms 4,384 KB
testcase_30	AC	1 ms 4,380 KB
testcase_31	AC	1 ms 4,376 KB
testcase_32	AC	1 ms 4,384 KB
testcase_33	AC	1 ms 4,376 KB
testcase_34	AC	3 ms 4,380 KB
testcase_35	AC	2 ms 4,380 KB
testcase_36	AC	3 ms 4,380 KB
testcase_37	AC	3 ms 4,376 KB
testcase_38	AC	2 ms 4,380 KB
testcase_39	AC	2 ms 4,380 KB
testcase_40	AC	3 ms 4,376 KB
testcase_41	AC	3 ms 4,376 KB
testcase_42	AC	2 ms 4,380 KB
testcase_43	AC	2 ms 4,380 KB
testcase_44	AC	2 ms 4,380 KB
testcase_45	AC	3 ms 4,380 KB
testcase_46	AC	3 ms 4,376 KB
testcase_47	AC	3 ms 4,380 KB
testcase_48	AC	2 ms 4,380 KB
testcase_49	AC	2 ms 4,380 KB
testcase_50	AC	3 ms 4,380 KB
testcase_51	AC	3 ms 4,376 KB
testcase_52	AC	3 ms 4,380 KB
testcase_53	AC	2 ms 4,380 KB
testcase_54	AC	2 ms 4,376 KB
testcase_55	AC	3 ms 4,376 KB
testcase_56	AC	3 ms 4,380 KB
testcase_57	AC	2 ms 4,380 KB
testcase_58	AC	2 ms 4,380 KB
testcase_59	AC	3 ms 4,384 KB
testcase_60	AC	2 ms 4,380 KB
testcase_61	AC	1 ms 4,380 KB
testcase_62	AC	2 ms 4,380 KB
testcase_63	AC	2 ms 4,380 KB
testcase_64	AC	1 ms 4,380 KB
testcase_65	AC	3 ms 4,380 KB
testcase_66	AC	2 ms 4,380 KB
testcase_67	AC	3 ms 4,380 KB
testcase_68	AC	1 ms 4,380 KB
testcase_69	AC	2 ms 4,376 KB
testcase_70	AC	3 ms 4,380 KB
testcase_71	AC	2 ms 4,376 KB
testcase_72	AC	1 ms 4,376 KB
testcase_73	AC	2 ms 4,376 KB
testcase_74	AC	3 ms 4,376 KB
testcase_75	AC	1 ms 4,376 KB
testcase_76	AC	1 ms 4,380 KB
testcase_77	AC	2 ms 4,376 KB
testcase_78	AC	2 ms 4,376 KB
testcase_79	AC	2 ms 4,376 KB
testcase_80	AC	2 ms 4,380 KB
testcase_81	AC	3 ms 4,380 KB
testcase_82	AC	1 ms 4,380 KB
testcase_83	AC	2 ms 4,380 KB
testcase_84	AC	1 ms 4,380 KB
testcase_85	AC	2 ms 4,376 KB
testcase_86	AC	2 ms 4,376 KB
testcase_87	AC	2 ms 4,380 KB
testcase_88	AC	2 ms 4,380 KB
testcase_89	AC	2 ms 4,376 KB
testcase_90	AC	1 ms 4,376 KB
testcase_91	AC	2 ms 4,380 KB
testcase_92	AC	2 ms 4,380 KB
testcase_93	AC	1 ms 4,376 KB
testcase_94	AC	2 ms 4,376 KB
testcase_95	AC	2 ms 4,380 KB
testcase_96	AC	1 ms 4,384 KB
testcase_97	AC	1 ms 4,380 KB
testcase_98	AC	2 ms 4,380 KB
testcase_99	AC	1 ms 4,380 KB
testcase_100	AC	1 ms 4,380 KB
testcase_101	AC	2 ms 4,376 KB
testcase_102	AC	1 ms 4,380 KB
testcase_103	AC	1 ms 4,376 KB
testcase_104	AC	1 ms 4,380 KB
testcase_105	AC	1 ms 4,380 KB
testcase_106	AC	2 ms 4,380 KB
testcase_107	AC	2 ms 4,376 KB
testcase_108	AC	1 ms 4,380 KB
testcase_109	AC	1 ms 4,380 KB
testcase_110	AC	2 ms 4,380 KB
testcase_111	AC	2 ms 4,384 KB
testcase_112	AC	1 ms 4,376 KB
testcase_113	AC	1 ms 4,380 KB
testcase_114	AC	1 ms 4,376 KB
testcase_115	AC	1 ms 4,380 KB
testcase_116	AC	2 ms 4,380 KB
testcase_117	AC	2 ms 4,380 KB
testcase_118	AC	1 ms 4,376 KB
testcase_119	AC	1 ms 4,380 KB
testcase_120	AC	1 ms 4,380 KB
testcase_121	AC	1 ms 4,376 KB
testcase_122	AC	2 ms 4,376 KB
testcase_123	AC	3 ms 4,380 KB
testcase_124	AC	2 ms 4,384 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
//#include<boost/multiprecision/cpp_int.hpp>
#include<boost/multiprecision/cpp_dec_float.hpp>
//#include <atcoder/all>
#define rep(i, a) for (int i = (int)0; i < (int)a; ++i)
#define rrep(i, a) for (int i = (int)a - 1; i >= 0; --i)
#define REP(i, a, b) for (int i = (int)a; i < (int)b; ++i)
#define RREP(i, a, b) for (int i = (int)a - 1; i >= b; --i)
#define repl(i, a) for (ll i = (ll)0; i < (ll)a; ++i)
#define pb push_back
#define eb emplace_back
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define popcount __builtin_popcount
#define fi first
#define se second
using ll = long long;
constexpr ll mod = 1e9 + 7;
constexpr ll mod_998244353 = 998244353;
constexpr ll INF = 1LL << 60;

// #pragma GCC target("avx2")
// #pragma GCC optimize("O3")
// #pragma GCC optimize("unroll-loops")

//using lll=boost::multiprecision::cpp_int;
using Double=boost::multiprecision::number<boost::multiprecision::cpp_dec_float<128>>;//仮数部が1024桁
template <class T>
inline bool chmin(T &a, T b)
{
	if (a > b)
	{
		a = b;
		return true;
	}
	return false;
}
template <class T>
inline bool chmax(T &a, T b)
{
	if (a < b)
	{
		a = b;
		return true;
	}
	return false;
}

template <typename T>
T mypow(T x, T n, const T &p = -1)
{ //x^nをmodで割った余り

	if (p != -1)
	{
		x %= p;
	}
	T ret = 1;
	while (n > 0)
	{
		if (n & 1)
		{
			if (p != -1)
				ret = (ret * x) % p;
			else
				ret *= x;
		}
		if (p != -1)
			x = (x * x) % p;
		else
			x *= x;
		n >>= 1;
	}
	return ret;
}

using namespace std;
//using namespace atcoder;


void solve()
{	
  ll a,b,c;
  cin>>a>>b>>c;
  auto D=[](ll a,ll b,ll c)->ll{
    return b*b-a*c*4;
  };
  if(!a){
    if(!b){
      if(!c){
        cout<<-1<<"\n";
      }else{
        cout<<0<<"\n";
      }
    }else{
      cout<<1<<"\n";
      Double bb=b,cc=c;
      Double x=-cc/bb;
      cout<<x<<"\n";
    }
  }
  else if(D(a,b,c)<0){
    cout<<0<<"\n";
  }else if(D(a,b,c)==0){
    cout<<1<<"\n";
    Double aa=a,bb=b;
    Double x=-bb/(aa*2);
    cout<<x<<"\n";
  }else{
    cout<<2<<"\n";
    Double l=0,r=3000000000;
    Double aa=a,bb=b,cc=c;
    rep(i,500){
      Double mid=(l+r)/2;
      if(mid*mid<=bb*bb-aa*cc*4)l=mid;
      else r=mid;
    }
    Double x=l;
    Double x1=0;
    if(b>=0)x1=(-bb-x)/(aa*2);
    else x1=(-bb+x)/(aa*2);
    Double x2=-bb/aa-x1;
    if(x1>x2)swap(x1,x2);
    cout<<x1<<"\n"<<x2<<"\n";
  }
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	cout << fixed << setprecision(50);
	solve();
	return 0;
}

yukicoder

結果

テストケース

ソースコード