[CSP-J 2023] 一元二次方程【民间数据】（模拟题）

[CSP-J 2023] 一元二次方程【民间数据】

题目背景

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

• 计算 Δ=𝑏2−4𝑎𝑐Δ=b2−4ac，则:
  1. 若 Δ<0Δ<0，则该一元二次方程无实数解。
  2. 否则 Δ≥0Δ≥0，此时该一元二次方程有两个实数解 𝑥1,2=−𝑏±Δ2𝑎x1,2​=2a−b±Δ​​。

例如：

• 𝑥2+𝑥+1=0x2+x+1=0 无实数解，因为 Δ=12−4×1×1=−3<0Δ=12−4×1×1=−3<0。
• 𝑥2−2𝑥+1=0x2−2x+1=0 有两相等实数解 𝑥1,2=1x1,2​=1。
• 𝑥2−3𝑥+2=0x2−3x+2=0 有两互异实数解 𝑥1=1,𝑥2=2x1​=1,x2​=2。

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

题目描述

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

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

• 由有理数的定义，存在唯一的两个整数 𝑝p 和 𝑞q，满足 𝑞>0q>0，gcd⁡(𝑝,𝑞)=1gcd(p,q)=1 且 𝑣=𝑝𝑞v=qp​。

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

• 例如：

  • 当 𝑣=−0.5v=−0.5 时，𝑝p 和 𝑞q 的值分别为 −1−1 和 22，则应输出 -1/2；
  • 当 𝑣=0v=0 时，𝑝p 和 𝑞q 的值分别为 00 和 11，则应输出 0。

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

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

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

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

   2. 否则根据上文公式，𝑥x 可以被唯一表示为 𝑥=𝑞1+𝑞2𝑟x=q1​+q2​r​ 的形式，其中：

      • 𝑞1,𝑞2q1​,q2​ 为有理数，且 𝑞2>0q2​>0；
      • 𝑟r 为正整数且 𝑟>1r>1，且不存在正整数 𝑑>1d>1 使 𝑑2∣𝑟d2∣r（即 𝑟r 不应是 𝑑2d2 的倍数）；

   此时：

   1. 若 𝑞1≠0q1​=0，则按有理数的格式输出 𝑞1q1​，并再输出一个加号 +；
   2. 否则跳过这一步输出；

   随后：

   1. 若 𝑞2=1q2​=1，则输出 sqrt({r})；
   2. 否则若 𝑞2q2​ 为整数，则输出 {q2}*sqrt({r})；
   3. 否则若 𝑞3=1𝑞2q3​=q2​1​ 为整数，则输出 sqrt({r})/{q3}；
   4. 否则可以证明存在唯一整数 𝑐,𝑑c,d 满足 𝑐,𝑑>1,gcd⁡(𝑐,𝑑)=1c,d>1,gcd(c,d)=1 且 𝑞2=𝑐𝑑q2​=dc​，此时输出 {c}*sqrt({r})/{d}；

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

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

输入格式

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

接下来 𝑇T 行，每行包含三个整数 𝑎,𝑏,𝑐a,b,c。

输出格式

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

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

样例 #1

样例输入 #1

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

Copy

样例输出 #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

Copy

提示

【样例 #2】

见附件中的 uqe/uqe2.in 与 uqe/uqe2.ans。

【数据范围】

对于所有数据有：1≤𝑇≤50001≤T≤5000，1≤𝑀≤1031≤M≤103，∣𝑎∣,∣𝑏∣,∣𝑐∣≤𝑀∣a∣,∣b∣,∣c∣≤M，𝑎≠0a=0。

测试点编号	𝑀≤M≤	特殊性质 A	特殊性质 B	特殊性质 C
11		是
22	2020	否	否	否
33	103103	是		是
44				否
55		否	是	是
66				否
7,87,8			否	是
9,109,10				否

其中：

• 特殊性质 A：保证 𝑏=0b=0；
• 特殊性质 B：保证 𝑐=0c=0；
• 特殊性质 C：如果方程有解，那么方程的两个解都是整数。

思路

直接模拟！！！

代码

#include<bits/stdc++.h> 
using namespace std;
int n,m,a,b,c,d,k,t;
int main(){
    cin>>n>>m;
    for(int i=1;i<=n;i++){
        cin>>a>>b>>c;
        if(a<0){
            a=-a;
            b=-b;
            c=-c;
        }            
        d=b*b-4*a*c;
        k=1;
        if(d<0){
            cout<<"NO"<<endl;
            continue;     
        }
        for(int i=2;i*i<=d;i++){
            while(d%(i*i)==0){
                k*=i;
                d/=(i*i);        
            }                    
        }
        if(d==0||d==1){
            t=abs(__gcd(2*a,-b+k*d));
            cout<<(-b+k*d)/t;
            if(2*a/t!=1) cout<<"/"<<2*a/t;
            cout<<endl;
            continue;     
        }
        t=abs(__gcd(-b,2*a));
        if(-b/t!=0) cout<<-b/t;            
        if(2*a/t!=1) cout<<"/"<<2*a/t;        
        if(-b/t!=0) cout<<"+";
        t=abs(__gcd(k,2*a));
        if(k/t!=1) cout<<k/t<<"*";    
        cout<<"sqrt("<<d<<")";
        if(2*a/t!=1) cout<<"/"<<2*a/t;
        cout<<endl;
    }
    return 0;
}