#include <stdio.h>

#include <math.h>

struct POINT{

double x;

double y; };

typedef struct POINT Point;

#define PI 3.14

double dist( Point p1, Point p2 );

void GetBorder( int, Point* AllTree, double *, double *, double *, double * );

int InnerPoint( Point , double, double, double, double, double );

int InnerPoint( Point p1, double Circle, double Left, double Right, double Top, double Bottom )

{

double r;

r = Circle / 2;

if( fabs(p1.x - Left) <= r || fabs(p1.x - Right) <= r ||

fabs(p1.y - Top) <= r || fabs(p1.y - Bottom) <= r )

{

return 0;//boundary point

}

return 1;//interior point

}

double dist( Point p1, Point p2 )

{

double distance;

float f;

f = (p1.x-p2.x)*(p1.x-p2.x) + (p1.y-p2.y)*(p1.y-p2.y);

distance = sqrt( f );

return distance;

}

void GetBorder( int Num, Point *AllTree, double *Left, double *Right, double *Top, double *Bottom )

{

int i;

double minLeft, maxRight, maxTop, minBottom;

minLeft = 1000;

maxRight = -1000;

maxTop = -1000;

minBottom = 1000;

for( i = 0; i < Num; i++ )

{

if( AllTree[i].x <= minLeft )

{

minLeft = AllTree[i].x;

}

if( AllTree[i].x >= maxRight )

{

maxRight = AllTree[i].x;

}

if( AllTree[i].y <= minBottom )

{

minBottom = AllTree[i].y;

}

if( AllTree[i].y >= maxTop )

{

maxTop = AllTree[i].y;

}

}

*Left = minLeft;

*Right = maxRight;

*Top = maxTop;

*Bottom = minBottom;

return;

}

int main( )

{

FILE *file;

FILE * file1;

int num;

double *pCircle;

double *pSize;

Point *pTree;

int *pGrade;

int i;

int j;

double Left, Right, Top, Bottom;

double *pPressure, r, r2, distance,temp;

file = fopen( "//root//Project//ganranyali//yali//sample5.txt", "r" );

file1 = fopen( "//root//Project//ganranyali//yali//pressure5.txt", "w" );

if( NULL == file )

{

printf( "fail to open file!\n" );

return 0;

}

fscanf( file, "%d\n", &num );

pCircle = (double*)malloc( num * sizeof(double) );

pSize = (double*)malloc( num * sizeof(double) );

pTree = (Point*)malloc( num*sizeof(Point) );

pGrade = (int*)malloc( num*sizeof(int) );

pPressure = (double*)malloc( num*sizeof(double) );

for( i = 0; i < num; i++ )

{

fscanf( file, "%lf ", &temp );

pCircle[i] = temp;

fscanf( file, "%lf ", &temp );

pSize[i] = temp;

fscanf( file, "%lf ", &temp );

pTree[i].x = temp ;

fscanf( file, "%lf ", &temp);

pTree[i].y = temp;

fscanf( file, "%d\n", &pGrade[i] );

// pGrade[i] = temp;

}

fclose( file );

GetBorder( num, pTree, &Left, &Right, &Top, &Bottom );

printf( "Left = %f\n Right = %f\n Top = %f\n Bottom = %f\n", Left, Right, Top, Bottom );

for( i = 0; i < num; i++ )

{

pPressure[i] = -1;

}

for( i = 0; i < num; i++ )

{

if( i == 13 )

{

j = 0;

}

r = pCircle[i] / 2;

if( InnerPoint( pTree[i], pCircle[i], Left, Right, Top, Bottom ) )

{

pPressure[i] = 0;

for( j = 0; j < num; j++ )

{

if( j != i )

{

distance = dist(pTree[i], pTree[j]);

r2 = pCircle[j] / 2;

if( distance <= 0.000001 )

{

distance = 0.1;

}

if( distance <= (r + r2) )

{

if( i == 23 & pGrade[j] != 0 )

printf( "%d\n", j );

pPressure[i] += (pSize[i]*pSize[j]*pGrade[j]) / ( 10000 * PI * PI *distance );

}

}

}

}

}

//print result of pressure

if( NULL == file1 )

{

printf( "fail to open output file\n" );

return 0;

}

for( i = 0; i < num; i++ )

{

fprintf( file1, "%d\t", i+1 );

fprintf( file1, "%f\n", pPressure[i] );

}

fclose( file1 );

}