2100：【23CSPJ普及组】一元二次方程(uqe)

【题目描述】

众所周知，对一元二次方程 ax2+bx+c=0,(a≠0)��2+��+�=0,(�≠0)，可以用以下方式求实数解：

∙∙ 计算 Δ=b2−4acΔ=�2−4��，则:

1. 若 Δ<0Δ<0，则该一元二次方程无实数解。

2. 否则 Δ≥0Δ≥0，此时该一元二次方程有两个实数解 x1,2=−b±Δ√2a�1,2=−�±Δ2�。

例如：

∙x2+x+1=0∙�2+�+1=0 无实数解，因为 Δ=12−4×1×1=−3<0Δ=12−4×1×1=−3<0。

∙x2−2x+1=0∙�2−2�+1=0 有两相等实数解 x1,2=1�1,2=1。

∙x2−3x+2=0∙�2−3�+2=0 有两互异实数解 x1=1,x2=2�1=1,�2=2。

在题面描述中 a� 和 b� 的最大公因数使用 gcd(a,b)gcd(�,�) 表示。例如 1212 和 1818 的最大公因数是 66，即 gcd(12,18)=6gcd(12,18)=6。

现在给定一个一元二次方程的系数 a,b,c�,�,�，其中 a,b,c�,�,� 均为整数且 a≠0�≠0。你需要判断一元二次方程 ax2+bx+c=0��2+��+�=0 是否有实数解，并按要求的格式输出。

在本题中输出有理数 v� 时须遵循以下规则：

∙∙ 由有理数的定义，存在唯一的两个整数 p� 和 q�，满足 q>0�>0，gcd(p,q)=1gcd(�,�)=1 且 v=pq�=��。

∙∙ 若 q=1�=1，则输出 {p}，否则输出 {p}/{q}，其中 {n} 代表整数 n� 的值；

∙∙ 例如：

∙∙ 当 v=−0.5�=−0.5 时，p� 和 q� 的值分别为 −1−1 和 22，则应输出 -1/2；

∙∙ 当 v=0�=0 时，p� 和 q� 的值分别为 00 和 11，则应输出 0。

对于方程的求解，分两种情况讨论：

1. 若 Δ=b2−4ac<0Δ=�2−4��<0，则表明方程无实数解，此时你应当输出 NO；

2. 否则 Δ≥0Δ≥0，此时方程有两解（可能相等），记其中较大者为 x�，则：

(1). 若 x� 为有理数，则按有理数的格式输出 x�。

(2). 否则根据上文公式，x� 可以被唯一表示为 x=q1+q2r√�=�1+�2� 的形式，其中：

     ∙q1,q2∙�1,�2 为有理数，且 q2>0�2>0；

     ∙r∙� 为正整数且 r>1�>1，且不存在正整数 d>1�>1 使 d2∣r�2∣�（即 r� 不应是 d2�2 的倍数）；

  此时：

  1. 若 q1≠0�1≠0，则按有理数的格式输出 q1�1，并再输出一个加号 +；

  2. 否则跳过这一步输出；

  随后：

  1. 若 q2=1�2=1，则输出 sqrt({r})；

  2. 否则若 q2�2 为整数，则输出 {q2}*sqrt({r})；

  3. 否则若 q3=1q2�3=1�2 为整数，则输出 sqrt({r})/{q3}；

  4. 否则可以证明存在唯一整数 c,d�,� 满足 c,d>1,gcd(c,d)=1�,�>1,gcd(�,�)=1 且 q2=cd�2=��，此时输出{c}*sqrt({r})/{d}；

  上述表示中 {n} 代表整数 {n} 的值，详见样例。

  如果方程有实数解，则按要求的格式输出两个实数解中的较大者。否则若方程没有实数解，则输出 NO。

【输入】

输入的第一行包含两个正整数 T,M�,�，分别表示方程数和系数的绝对值上限。

接下来 T� 行，每行包含三个整数 a,b,c�,�,�。

【输出】

输出 T� 行，每行包含一个字符串，表示对应询问的答案，格式如题面所述。

每行输出的字符串中间不应包含任何空格。

【输入样例】

9 1000
1 -1 0
-1 -1 -1
1 -2 1
1 5 4
4 4 1
1 0 -432
1 -3 1
2 -4 1
1 7 1

【输出样例】

1
NO
1
-1
-1/2
12*sqrt(3)
3/2+sqrt(5)/2
1+sqrt(2)/2
-7/2+3*sqrt(5)/2

【提示】

【数据范围】

对于所有数据有：1≤T≤50001≤�≤5000，1≤M≤1031≤�≤103，∣a∣,∣b∣,∣c∣≤M∣�∣,∣�∣,∣�∣≤�，a≠0�≠0。

测试点编号	M≤	特殊性质 A	特殊性质 B	特殊性质 C
1	1	是	是	是
2	20	否	否	否
3	103	是	否	是
4	103	是	否	否
5	103	否	是	是
6	103	否	是	否
7,8	103	否	否	是
9,10	103	否	否	否

其中：

特殊性质 A：保证 b=0；

特殊性质 B：保证 c=0；

特殊性质 C：如果方程有解，那么方程的两个解都是整数。

废话不多说，AC代码启动！

#include<bits/stdc++.h>
using namespace std;
int t,m,a,b,c;
int gcd(int p,int q){
    int u;
    while(q!=0){
        u=p%q;
        p=q;
        q=u;
    }
    return p;
}
void work1(int p,int q){
    if(p*q<0) cout<<'-';
    if(p<0) p=0-p;
    if(q<0) q=0-q;
    int g=gcd(p,q);
    if(g==q) cout<<p/g;	
    else cout<<p/g<<'/'<<q/g;
    return;
}
void work2(int p,int q){
    if(p<0) p=-p;
    if(q<0) q=-q;
    int u=1;
    for(int i=sqrt(p); i>=2; i--) {
        if(p%(i*i)==0){
            p=p/(i*i);
            u=u*i;
            break;
        }
    }
    int g=gcd(u,q);
    u=u/g;
    q=q/g;
    if(u!=1) cout<<u<<"*";
    if(p!=1) cout<<"sqrt("<<p<<")";
    if(q!=1) cout<<"/"<<q;
}
int main(){
//    freopen("uqe2.in","r",stdin);
//    freopen("uqe2.ans","w",stdout);
    cin>>t>>m;
    while(t--) {
        cin>>a>>b>>c;
        int d=b*b-4*a*c;
        if(d<0){
            cout<<"NO";
        }else{ 
            int q=sqrt(d);
            if(q*q==d) { 
                double ans1=1.0*(-b+q)/(2*a);
                double ans2=1.0*(-b-q)/(2*a);
                if(ans1>ans2) {
                    if(-b+q==0) cout<<0;
                    else work1(-b+q,2*a);
                } else {
                    if(-b-q==0) cout<<0;
                    else work1(-b-q,2*a);
                }
            } else { 
                if(b!=0) work1(-b,2*a);
                if(b!=0) cout<<"+";
                work2(d,2*a);
            }

        }
        cout<<"\n";
    }
    return 0;
}

给个关注，www