// Pifagor.cpp : Defines the entry point for the console application.
//

#pragma warning( disable : 4786)
#include <stdio.h>
#include <algorithm>
#include <math.h>
#include <windows.h>
#include <map>


// Greatest Common Divisor
inline unsigned int gcd(unsigned int a, unsigned int b){
  unsigned int c;
  while ( a != 0 ) {
     c = a; a = b%a;  b = c;
  }
  return b;
}

// для хранения ключа
typedef struct {
	unsigned int a;
	unsigned int b;
} sPifagor_t;

struct clsPifagorCompare:
public                                 
std::binary_function<sPifagor_t, sPifagor_t, bool>
{ 
	bool operator() (const sPifagor_t& lhs, const sPifagor_t& rhs) const 
	{
		return lhs.a < rhs.a || (lhs.a == rhs.a && lhs.b < rhs.b);
	}
};

typedef std::map<sPifagor_t, unsigned int, clsPifagorCompare> map_t;
map_t mapRes;
map_t::iterator cur;

unsigned int iStore = 0;

inline void StoreTri(unsigned int a, unsigned int b, unsigned int c)
{
	sPifagor_t sPifagor;
	unsigned int iGCD = gcd(a, b);
	if(iGCD > 1)
	{
		a = a / iGCD;
		b = b / iGCD;
		c = c / iGCD;
	}
	if(a > b)
	{
		sPifagor.a = b;
		sPifagor.b = a;
	}
	else
	{
		sPifagor.a = a;
		sPifagor.b = b;
	}
	
	cur = mapRes.find(sPifagor);
	if(cur != mapRes.end())
	{
		if((*cur).second != c)
		{
			printf("error(%d,%d,%d != %d)\n", a, b, c, (*cur).second);
			MessageBox(NULL, "error", "error", 0);
		}
	}
	else
	{
		mapRes.insert(std::make_pair(sPifagor, c));
		iStore++;
	}

}

inline unsigned int relationPrime(unsigned int m, unsigned int n){
	return gcd(m,n)==1?1:0;
}

// Stuped brute-force pyfagorean triangle generator
int main(int argc,char **argv){

	if (argc>1){
		unsigned int iCount = 0, iNext = 1000000, iStart = GetTickCount();
		unsigned int maxside = atoi(argv[1]);
		for(int i=1;i<maxside;i++){
		    unsigned int ii = i * i;
	            for(int j=1;j<i;j++){
			unsigned int jj = j * j;
	   		if (relationPrime(i,j)) {
				StoreTri(ii - jj,2 * i * j,ii + jj);
				iCount++;
				if(iCount >= iNext) 
				{
					printf("Parse %d Store %d Time %d msec\n", iCount, iStore, GetTickCount() - iStart);
					iNext += 1000000;
				}

			}        
	            }
	            
	        }
		printf("Parse %d Store %d Time %d msec\n", iCount, iStore, GetTickCount() - iStart);
		return 0;
	} else {
		fprintf(stderr,"Pyfagorean Triangle Generator 1.0 of 1 jan 2014, by Mayton\n\nUsage: pyfagor [N]\n\n");
		return 1;	
	}

	
	return 0;
}

// Stuped brute-force pyfagorean triangle generator
int main(int argc,char **argv){

	if (argc>1){
		unsigned int iCount = 0, iNext = 1000000, iStart = GetTickCount();
		unsigned int maxside = atoi(argv[1]);
		for(int i=1;i<maxside;i++){
		    unsigned int ii = i * i;
	            for(int j=1;j<i;j++){
			unsigned int jj = j * j;
	   		if (relationPrime(i,j)) {
				StoreTri(ii - jj,2 * i * j,ii + jj);
				iCount++;
				if(iCount >= iNext) 
				{
					printf("Parse %d Store %d Time %d msec\n", iCount, iStore, GetTickCount() - iStart);
					iNext += 1000000;
				}

			}        
	            }
	            
	        }
		printf("Parse %d Store %d Time %d msec\n", iCount, iStore, GetTickCount() - iStart);
		for(unsigned int a = 1; a <40000; a++)
		{
			for(unsigned int b = 1; b < 40000; b++)
			{
				unsigned int c2 = a*a + b*b;
				unsigned int c = sqrt(c2);
				if(c2 == c*c)
				{
					unsigned int iStorePrev = iStore;
					StoreTri(a, b, c);
					if(iStorePrev != iStore)
					{ 
						printf("skip(%d,%d,%d)\n", a, b, c);
						iCount++;
					}
				}

			}
		}
		printf("Parse %d Store %d Time %d msec\n", iCount, iStore, GetTickCount() - iStart);
		map2file();
		return 0;
	} else {
		fprintf(stderr,"Pyfagorean Triangle Generator 1.0 of 1 jan 2014, by Mayton\n\nUsage: pyfagor [N]\n\n");
		return 1;	
	}

	
	return 0;
}

inline void StoreTri(unsigned int a, unsigned int b, unsigned int c)
{
	sPifagor_t sPifagor;

	unsigned int iGCD = gcd(a, b);
	if(iGCD > 1)
	{
		a = a / iGCD;
		b = b / iGCD;
		c = c / iGCD;
	}

	if(a > b)
	{
		sPifagor.a = b;
		sPifagor.b = a;
	}
	else
	{
		sPifagor.a = a;
		sPifagor.b = b;
	}
	if(sPifagor.b >= 10000000) return;
....
}