글 읽기 - 반례좀 알 수 있을까요??

nlkey2022 6년 전 0

반례가 뭐가 있을까요?..

#include <iostream>
#include <cmath>
using namespace std;


int main() {
	

	int testcase = 0;
	cin >> testcase;
	
	double result[100][3];
	
	for(int i=0; i<100; i++)
		for(int j=0; j<3; j++)
			result[i][j] = -99999999999999999999999;
	

	
	for(int ttt=0; ttt<testcase; ttt++){
		double a,b,c,d;
		cin >>a>>b>>c>>d;
		
		
		if(d == 0){
			double x1 = 0;
			double x2 = (-b+sqrt(b*b-4*a*c))/(2*a);
			double x3 = (-b-sqrt(b*b-4*a*c))/(2*a);
			//	cout <<"x2 : "<<x2 <<" x3 : "<<x3<<endl;
				if(x1 > x2){
					if(x2 > x3){
						// x1 > x2 > x3
						result[ttt][0] = x3;
						result[ttt][1] = x2;
						result[ttt][2] = x1;
					}else if(x2 < x3){
						if(x1 > x3){
							// x1 > x3 > x2
							result[ttt][0] = x2;
							result[ttt][1] = x3;
							result[ttt][2] = x1;
						}else if(x1 < x3){
							// x3 > x1 > x2
							result[ttt][0] = x2;
							result[ttt][1] = x1;
							result[ttt][2] = x3;	
						}
					}
					else{
						// x2 == x3 < 0
						result[ttt][0] = x2;
						result[ttt][1] = 0;
					}
				}else if(x1 < x2){
					if(x2 < x3){
						// x3 > x2 > x1
						result[ttt][0] = x1;
						result[ttt][1] = x2;
						result[ttt][2] = x3;
					}else if(x2 > x3){
						if(x1 > x3){
							// x2 > x1 > x3
							result[ttt][0] = x3;
							result[ttt][1] = x1;
							result[ttt][2] = x2;
						}else if(x1 < x3){
							// x2 > x3 > x1
							result[ttt][0] = x1;
							result[ttt][1] = x3;
							result[ttt][2] = x2;
						}
					}
					else{
						// x2 == x3 > 0
						result[ttt][0] = 0;
						result[ttt][1] = x2;
					}
				}
				else{
					//x1==x2==0
					if(x3 > 0){
						result[ttt][0] = 0;
						result[ttt][1] = x3;	
					}else if(x3 < 0){
						result[ttt][0] = x3;
						result[ttt][1] = 0;
					}else{
						//x1==x2==x3
						result[ttt][0] = 0;
					}
				}
		}
		else if(d != 0){
		
			b /= a;
			c /= a;
			d /= a;
			double disc, q, r, dum1, s, t, term1, r13; 
			q = (3.0*c - (b*b))/9.0;
			r = -(27.0*d) + b*(9.0*c - 2.0*(b*b));
			r /= 54.0;
			
			
			disc = q*q*q + r*r;
			term1 = (b/3.0);
			if (disc > 0) { // one root real, two are complex
			
				s = r + sqrt(disc);
				s = ((s < 0) ? -pow(-s, (1.0/3.0)) : pow(s, (1.0/3.0)));
				t = r - sqrt(disc);
				t = ((t < 0) ? -pow(-t, (1.0/3.0)) : pow(t, (1.0/3.0)));
				
				double x = -term1 + s + t;
				
				result[ttt][0] = x;
			}
			else if(disc == 0){
				r13 = ((r < 0) ? -pow(-r,(1.0/3.0)) : pow(r,(1.0/3.0)));
				double x1 = -term1 + 2.0*r13;
				double x2 = -(r13 + term1);
				
	
				
	
				if(x1>x2) {
					result[ttt][0] = x2;
					result[ttt][1] = x1;
				}
				else if(x1 < x2){
					result[ttt][0] = x1;
					result[ttt][1] = x2;
				}
				else result[ttt][0] = 0;
	
			}
			else{
				q = -q;
				dum1 = q*q*q;
				dum1 = acos(r/sqrt(dum1));
				r13 = 2.0*sqrt(q);
				
				double x1,x2,x3;
				x1 = -term1 + r13*cos(dum1/3.0);
				x2 = -term1 + r13*cos((dum1 + 2.0*M_PI)/3.0);
				x3 = -term1 + r13*cos((dum1 + 4.0*M_PI)/3.0);
				
	
				
				if(x1 > x2){
					if(x2 > x3){
						// x1 > x2 > x3
						result[ttt][0] = x3;
						result[ttt][1] = x2;
						result[ttt][2] = x1;
					}else if(x2 < x3){
						if(x1 > x3){
							// x1 > x3 > x2
							result[ttt][0] = x2;
							result[ttt][1] = x3;
							result[ttt][2] = x1;
						}else if(x1 < x3){
							// x3 > x1 > x2
							result[ttt][0] = x2;
							result[ttt][1] = x1;
							result[ttt][2] = x3;	
						}
					}else{
						cout <<"is true?"<<endl;
					}
				}else if(x1 < x2){
					if(x2 < x3){
						// x3 > x2 > x1
						result[ttt][0] = x1;
						result[ttt][1] = x2;
						result[ttt][2] = x3;
					}else if(x2 > x3){
						if(x1 > x3){
							// x2 > x1 > x3
							result[ttt][0] = x3;
							result[ttt][1] = x1;
							result[ttt][2] = x2;
						}else if(x1 < x3){
							// x2 > x3 > x1
							result[ttt][0] = x1;
							result[ttt][1] = x3;
							result[ttt][2] = x2;
						}
					}
					else{
						cout <<"is true?"<<endl;
					}
				}
			}
		}
	}
	
	for(int ttt=0; ttt<testcase; ttt++){
		
		cout << fixed;
		cout.precision(4);
	
		for(int i=0; i<3; i++){
			if(result[ttt][i] != -99999999999999999999999){
			  	// minus zero
			  	if(abs(result[ttt][i] - 0) <= 0.000000000001) result[ttt][i] = 0;
			  	
				cout << result[ttt][i];
			
			}
			if(i!=2) {
				if(result[ttt][i+1] != -99999999999999999999999) cout <<" ";
			}
		}
		cout << endl;
	}
	
	return 0;
}

chogahui05 6년 전 0

이거 양아치 문제 같네요..

사실 푸는 방법은 정해져 있긴 합니다. Det로 근 갯수 판별하고

범위를 나눠서 이분탐색을 하면 되긴 됩니다만..

극대, 극솟값이 딱 떨어지지 않는다는 게 문제긴 합니다..

댓글을 작성하려면 로그인해야 합니다.

문제

문제

출처

ICPC

반례좀 알 수 있을까요??

nlkey2022 6년 전 0

chogahui05 6년 전 0

Baekjoon Online Judge

채점 현황

문제

유저 대회

고등학교 대회

출처

대학교 대회

도움말