結果

問題 No.1179 Quadratic Equation
ユーザー leafirbyleafirby
提出日時 2020-08-21 21:57:12
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,064 bytes
コンパイル時間 2,206 ms
コンパイル使用メモリ 206,716 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-15 05:37:01
合計ジャッジ時間 2,504 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 WA -
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 1 ms
5,248 KB
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 AC 2 ms
5,248 KB
testcase_12 WA -
testcase_13 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#pragma GCC target("avx")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")    
using namespace std;
#define lli long long int
#define uli unsigned long long int
#define INF 9999999999
#define rep(i,m,n) for(lli i = m;i < n;i++)
#define rrep(i,m,n) for(lli i=m-1;i>=n;i--)
#define pb(n) push_back(n)
#define UE(N) N.erase(unique(N.begin(),N.end()),N.end());
#define Sort(n) sort(n.begin(), n.end())
#define Rev(n) reverse(n.begin(),n.end())
#define Out(S) cout << S << endl
#define NeOut(S) cout << S
#define HpOut(S) cout << setprecision(50) << S << endl
#define Vec(K,L,N,S) vector<L> K(N,S)
#define DV(K,L,N,M,S) vector<vector<L>> K(N,vector<L>(M,S))
#define TV(K,L,N,M,R,S) vector<vector<vector<L>>> K(N,vector<vector<L>>(M,vector<L>(R,S)))
#define pint pair<lli,lli>
#define paf(L,R) pair<L,R>
#define mod 998244353
#define MAX 10000000
#define ALL(a)  a.begin(),a.end()
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))
#define __float128 _Float128
int main(){
  ios::sync_with_stdio(false);
  cin.tie(0);
  lli A,B,C,D,E,F,K,N,M,L,X,Y,Z,H,W,sum=0,num=0,flag=0;string S,T;
  cin >> A >> B >> C;
  Vec(ans,__float128,0,0.0);
    __float128 limit=(__float128)-B/(2*A);
    __float128 eps=0.000000000001;
    __float128 ok=limit;
    __float128 ng=-pow(10,9)*3;
    Z=103;
    while(Z--){
      __float128 mid=(ok+ng)/2.0;
      if(A>0&&A*mid*mid+B*mid+C>0)ng=mid;
      else if(A<0&&A*mid*mid+B*mid+C<0)ng=mid;
      else ok=mid;
    }
    if(A*ok*ok+B*ok+C<0.000000000001&&-0.000000000001<A*ok*ok+B*ok+C)ans.pb(ok);
    ok=limit;
    ng=pow(10,9)*3;
    Z=103;
    while(Z--){
      __float128 mid=(ok+ng)/2.0;
      if(A>0&&A*mid*mid+B*mid+C>0)ng=mid;
      else if(A<0&&A*mid*mid+B*mid+C<0)ng=mid;
      else ok=mid;
    }
    if(A*ok*ok+B*ok+C<0.000000000001&&-0.000000000001<A*ok*ok+B*ok+C)ans.pb(ok);
  if(ans.size()==2){
    if(ans[1]-ans[0]<0.000000000001)ans.pop_back();
  }
  if(ans.size())for(auto v:ans)HpOut((long double)v);
  else cout << "imaginary" << endl;
}
0