Cod sursa(job #1021338)

Utilizator danny794Dan Danaila danny794 Data 3 noiembrie 2013 18:26:15
Problema Algoritmul lui Euclid extins Scor 30
Compilator cpp Status done
Runda Arhiva educationala Marime 1.31 kb
#include <iostream>
#include <fstream>
#include <cstdlib>

using namespace std;

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

typedef struct pairs{ int x; int y; } pairs;

int cmmdc(int a, int b){
	if (a == 0 || b == 0)
		return a + b;
	if (a > b) return cmmdc(a % b, b);
	else return cmmdc(a, b % a);
}

pairs euclid_extended(int a, int b){
	pairs result;
	if (b == 1){
		result.x = 0;
		result.y = -1;
	} else if (a == 1){
		result.x = 1;
		result.y = 0;
	} else {
		int d, c;
		pairs rez;
		if(a > b){
			d = a / b;
			c = a % b;
			result = euclid_extended(b, c);
			rez.x = - result.y;
			rez.y = - (result.x + d * result.y);
		} else {
			d = b / a;
			c = b % a;
			result = euclid_extended(a, c);
			rez.y = result.y;
			rez.x = result.x + d * result.y;
		}
		return rez;
	}

	return result;
}

void solve(int a, int b, int c){
	int cm = cmmdc(abs(a), abs(b));
	if (c % cm != 0)
		g << "0 0";
	else{
		pairs p = euclid_extended(abs(a) / cm, abs(b) / cm);
		p.x = (a / abs(a)) * p.x * (c / cm);
		p.y = - (b / abs(b)) * p.y * (c / cm);
		g << p.x << " " << p.y << '\n';
	}
}

int main() {
	int n, a, b, c;
	f >> n;
	for(int i = 1; i <= n; i++){
		f >> a >> b >> c;
		solve(a, b, c);
//		pairs p = euclid_extended(abs(a), abs(b));
//		g << p.x << "*" << a << " - " << p.y << "*" << b;
	}
	f.close();
	g.close();
	return 0;
}