結果

問題 No.955 ax^2+bx+c=0
ユーザー ngtkanangtkana
提出日時 2020-03-28 01:45:07
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 4 ms / 2,000 ms
コード長 1,219 bytes
コンパイル時間 2,303 ms
コンパイル使用メモリ 200,732 KB
最終ジャッジ日時 2025-01-09 10:36:45
ジャッジサーバーID
(参考情報)
judge4 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 122
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In lambda function:
main.cpp:17:27: warning: narrowing conversion of ‘c’ from ‘const lint’ {aka ‘const long long int’} to ‘real’ {aka ‘long double’} [-Wnarrowing]
   17 |             return {-real{c}/b};
      |                           ^
main.cpp:17:27: warning: narrowing conversion of ‘(lint)c’ from ‘lint’ {aka ‘long long int’} to ‘real’ {aka ‘long double’} [-Wnarrowing]
main.cpp:20:22: warning: narrowing conversion of ‘b’ from ‘const lint’ {aka ‘const long long int’} to ‘real’ {aka ‘long double’} [-Wnarrowing]
   20 |         real x=-real{b}/(2*a);
      |                      ^
main.cpp:20:22: warning: narrowing conversion of ‘(lint)b’ from ‘lint’ {aka ‘long long int’} to ‘real’ {aka ‘long double’} [-Wnarrowing]

ソースコード

diff #

#include<bits/stdc++.h>
using lint=long long;
using real=long double;
int main(){
    std::cin.tie(nullptr);std::ios_base::sync_with_stdio(false);
    std::cout.setf(std::ios_base::fixed);std::cout.precision(30);
    lint a,b,c;std::cin>>a>>b>>c;
    auto ans=[a,b,c]()->std::vector<real>{
        if(a==0){
            if(b==0){
                if(c==0){
                    std::cout<<-1<<'\n';
                    exit(0);
                }
                return {};
            }
            return {-real{c}/b};
        }
        lint D=b*b-4*a*c;
        real x=-real{b}/(2*a);
        if(D<0)return {};
        if(D==0)return {x};
        real y=real{4}*a*c/(b*b);
        std::vector<real>ans(2);
        if(1-1e-6<y){
            ans.at(0)=(-b-std::sqrt(D))/(2*a);
            ans.at(1)=(-b+std::sqrt(D))/(2*a);
        }else if(std::abs(y)<1e-6){
            ans.at(0)=x*(1+std::sqrt(1-y));
            ans.at(1)=x*(y/2+y*y/8+y*y*y/16);
        }else{
            ans.at(0)=x*(1+std::sqrt(1-y));
            ans.at(1)=x*(1-std::sqrt(1-y));
        }
        return ans;
    }();
    std::sort(ans.begin(),ans.end());
    std::cout<<ans.size()<<'\n';
    for(real x:ans){
        std::cout<<x<<'\n';
    }
}

0