Cod sursa(job #631700)

Utilizator andunhillMacarescu Sebastian andunhill Data 9 noiembrie 2011 17:01:57
Problema Ecuatie Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.18 kb
#include<fstream>
#include<cmath>
using namespace std;

ifstream f("ecuatie.in");
ofstream g("ecuatie.out");

int A, B, C, K;

bool verf(int k, int p1, int p2, int q1, int q2)
{	if(k < K) return 0;

	g<<"(";
	if(p1 == 1 || p1 == -1) (p1 == 1) ? g<<"x" : g<<"-x";
	else g<<p1<<"x";
	(q1 < 0) ? g<<q1 : g<<"+"<<q1;
	g<<")";
	
	g<<"(";
	if(p2 == 1 || p2 == -1) (p2 == 1) ? g<<"x" : g<<"-x";
	else g<<p2<<"x";
	(q2 < 0) ? g<<q2 : g<<"+"<<q2;
	g<<")";
	
	return 1;
}

int main()
{	int k, p1, p2, q1, q2;
	int delta, sqrdel, sqrA;
	
	f>>A>>B>>C>>K;
	
	delta = B * B - 4 * A * C;
	if(delta < 0) { g<<-1; return 0; }
	
	sqrdel = sqrt((double)delta);
	sqrA = sqrt((double)A);
	
	for(p1 = -sqrdel , k = 0; p1 <= sqrdel; p1++)
	{	if(p1 == 0 || A % p1) continue;
		
		p2 = A / p1;
		
		if(delta == 0)
		{	q1 = B / (2 * p2);
			q2 = C / q1;
			k++;
			if(verf(k, p1, p2, q1, q2)) return 0;
		}
		else
		{	q1 = (B - sqrdel) / (2 * p2);
			q2 =C / q1;
			k++;
			if(verf(k, p1, p2, q1, q2)) return 0;
			
			q1 = (B + sqrdel) / (2 * p2);
			q2 =C / q1;
			k++;
			if(verf(k, p1, p2, q1, q2)) return 0;
		}
	}
	
	g<<-1;
	
	f.close();
	g.close();
	return 0;
}