#include <stdlib.h>
#define CORE_LEVEL 3   

#include "CORE/CORE.h" 

#define NUMOFLINE 50

class Line
{
public:
    double a,b,c;
};



Line generatedLines[NUMOFLINE+2];


void createLines()
{
    int i;
    for(i=0;i<NUMOFLINE;i++)
    {
       generatedLines[i].a=rand()%100+1;
       generatedLines[i].b=rand()%100+1;
       generatedLines[i].c=rand()%100+1;
    }
}



void intersectLines()
{
    int i,j;
	double PX, PY;
    double Delta,result;
    double correctResult,totalResult,r;

    correctResult=0; totalResult=0;
    for(i=0;i<NUMOFLINE;i++)
    {
       for(j=0;j<NUMOFLINE;j++)
       {
          if(i==j)continue;
          Delta=generatedLines[i].a*generatedLines[j].b
		  	-generatedLines[j].a*generatedLines[i].b;
		 if(Delta==0.0)continue;
          PX=(-generatedLines[j].b*generatedLines[i].c
			 +generatedLines[i].b*generatedLines[j].c)/Delta;
          PY=(generatedLines[j].a*generatedLines[i].c
			 -generatedLines[i].a*generatedLines[j].c)/Delta;
		 totalResult++;
          result=generatedLines[i].a*PX
		  	+generatedLines[i].b*PY+generatedLines[i].c;
          if(result==0.0)correctResult++;

		  //std::cout << totalResult << correctResult << std::endl;
       }
    }
	 r=correctResult/totalResult;
	 std::cout << "Correct percent " << r << std::endl;
}


int main()
{

     createLines();
     intersectLines();
     return(0);     
}