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

• 問題: No.955 ax^2+bx+c=0
• ユーザー: HIcoder
• 提出日時: 2023-07-18 16:08:16
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 2 ms / 2,000 ms
• コード長: 2,887 bytes
• コンパイル時間: 1,100 ms
• コンパイル使用メモリ: 102,348 KB
• 実行使用メモリ: 4,348 KB
• 最終ジャッジ日時: 2023-10-18 19:27:07
• 合計ジャッジ時間: 4,428 ms
• ジャッジサーバーID （参考情報）: judge13 / judge12

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

ソースコード

#include<iostream>
#include<set>
#include<algorithm>
#include<vector>
#include<string>
#include<set>
#include<map>
#include<numeric>
#include<queue>
#include<cmath>
using namespace std;
typedef long long ll;
const ll INF=1LL<<60;
typedef pair<ll,ll> P;
typedef pair<int,P> PP;
const ll MOD=1e9+7;

ll gcd(ll x,ll y){
    return y==0?x:gcd(y,x%y);
}

double f(ll n){
    const int N=100;
    double lb=0,ub=n+1;
    for(int t=0;t<N;t++){
        double mid=(lb+ub)/2;

        if(mid*mid<=n){
            lb=mid;
        }else{
            ub=mid;
        }
    }

    return lb;
}

void print(double x){
    printf("%.23f\n",x);
}

int main(){
    ll a,b,c;
    cin>>a>>b>>c;

    if(a==0){

        if(b==0){
            if(c==0){
                //a=b=c=0
                //解xは∞個
                cout<<-1<<endl;
            }else{
                //a=b=0, c!=0
                //解なし
                cout<<0<<endl;
            }
        }else{
            //b!=0
            // bx+c=0
            cout<<1<<endl;

            // ll g=gcd(b,c);
            // b/=g;
            // c/=g;

            {
                ll g=1;
                
                for(ll t=1;t*t<=min({b,c});t++ ){

                    if(b%t==0 && c%t==0){
                        g=t;
                    }

                }

                b/=g;
                c/=g;

            }



            double x=-1.0*c/b;
            //printf("%.12f\n",x);
            print(x);


        }
        
        return 0;
    }

    //ll g=gcd(c,gcd(b,a));
    {
        ll g=1;
        
        for(ll t=1;t*t<=min({a,b,c});t++ ){

            if(a%t==0 && b%t==0 && c%t==0){
                g=t;
            }

        }


        a/=g;
        b/=g;
        c/=g;

    }

    if(a<0){
        a*=(-1);
        b*=(-1);
        c*=(-1);
    }

    //a!=0

    ll D=b*b-4*a*c;

    if(D==0){
        cout<<1<<endl;


        double x=1.0*(-b)/(2*a);
        //printf("%.12f\n",x);
        print(x);


    }else if(D>0){

        double x1,x2;

        if(b>0){
            x1=(1.0*(-b)-sqrt(D))/(2*a);
            //x2=-1.0*b/a-x1;
            x2=1.0*c/(a*x1);

        }else{
            
            x1=(1.0*(-b)+sqrt(D))/(2*a);

            if(b*b==D){
                x1=0;
                x2=-1.0*b/a-x1;
            }else{
                //x1=0になる場合
                x2=1.0*c/(a*x1);
            }
           
            //x2=-1.0*b/a-x1;
            //x2=1.0*c/(a*x1);
        }
        
        //double x2=(1.0*(-b)+sqrt(D))/(2*a);
        //double x2=-1.0*b/a-x1;

        //double x1=(1.0*(-b)-f(D))/(2*a);
        //double x2=(1.0*(-b)+f(D))/(2*a);


       
        cout<<2<<endl;
        //printf("%.12f\n",x1);
        //printf("%.12f\n",x2);
        print(min(x1,x2));
        print(max(x1,x2));

    }else{
        //D<0
        //虚数解
        cout<<0<<endl;
    }
}