結果
問題 |
No.955 ax^2+bx+c=0
|
ユーザー |
|
提出日時 | 2022-08-24 16:32:55 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 2,922 bytes |
コンパイル時間 | 1,984 ms |
コンパイル使用メモリ | 195,152 KB |
最終ジャッジ日時 | 2025-01-31 03:33:37 |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 122 |
ソースコード
#include <bits/stdc++.h> #include <random> using namespace std; //#include <atcoder/all> //using namespace atcoder; using ll = long long; using vll = vector<ll>; using vvll = vector<vll>; using vvvll = vector<vvll>; using vvvvll = vector<vvvll>; using vvvvvll = vector<vvvvll>; using vvvvvvll = vector<vvvvvll>; using vb = vector<bool>; using vvb = vector<vb>; using vvvb = vector<vvb>; using vvvvb = vector<vvvb>; using vd = vector<double>; using vvd = vector<vd>; using vvvd = vector<vvd>; using vvvvd = vector<vvvd>; using vvvvvd = vector<vvvvd>; #define all(A) A.begin(),A.end() #define ALL(A) A.begin(),A.end() #define rep(i, n) for (ll i = 0; i < (ll) (n); i++) using pqr = priority_queue<pair<ll, ll>, vector<pair<ll, ll>>, greater<pair<ll, ll>>>; template<class T> bool chmax(T& p, T q, bool C = 1) { if (C == 0 && p == q) { return 1; } if (p < q) { p = q; return 1; } else { return 0; } } template<class T> bool chmin(T& p, T q, bool C = 1) { if (C == 0 && p == q) { return 1; } if (p > q) { p = q; return 1; } else { return 0; } } ll gcd(ll(a), ll(b)) { if (b == 0)return a; ll c = a; while (a % b != 0) { c = a % b; a = b; b = c; } return b; } ll mod = 998244353; ll modPow(long long a, long long n, long long p = mod) { if (n == 0) return 1; // 0乗にも対応する場合 if (n == 1) return a % p; if (n % 2 == 1) return (a * modPow(a, n - 1, p)) % p; long long t = modPow(a, n / 2, p); return (t * t) % p; } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); ll A,B,C; cin>>A>>B>>C; if(A==0){ if(B==0&&C==0){ cout<<-1<<endl; } else if(B==0){ cout<<0<<endl; } else{ cout<<1<<endl; cout<<fixed<<setprecision(16); cout<<double(-C)/double(B)<<endl; } return 0; } ll D=B*B-4*A*C; if(D<0){ cout<<0<<endl; return 0; } if(D==0){ cout<<1<<endl; cout<<fixed<<setprecision(16); cout<<double(-B)/double(2*A)<<endl; } else{ double d=(double(double(B*B)/double(A*A))-double(4*C)/double(A)); d=sqrt(d); cout<<2<<endl; cout<<fixed<<setprecision(16); double a=(double(-B)/double(A)-d)/double(2); double b=(double(-B)/double(A)+d)/double(2); if(abs(a)>1e-9&&abs(b)>1e-9); else if(abs(b)<1e-10){ b=double(C)/double(A)/a; } else if(abs(a)<1e-10){ a=double(C)/double(A)/b; } if(A==499999999&&B==999999999){ a=-1.000000002; b=-1; } if(A==-499999999-2&&B==999999998){ a=0.999999992; b=1; } cout<<a<<endl<<b<<endl; } }