Cod sursa(job #566371)

Utilizator maritimCristian Lambru maritim Data 28 martie 2011 21:51:30
Problema Trapez Scor 40
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.37 kb
#include<stdio.h>
#include<malloc.h>
#include<stdlib.h>

typedef struct
{
	double x;
	double y;
} xy;

xy A[1001];
xy B[1000001];
int N;
int nr = 0;
int nr1 = 0;
int nr2 = 0;
int acc = 0;

void citire(void)
{
	FILE *f = fopen("trapez.in","r");
	
	fscanf(f,"%d ",&N);
	for(int i=1;i<=N;i++)
		fscanf(f,"%d %d",&A[i].x,&A[i].y);
	
	fclose(f);
}

void add(int a,int b)
{
	B[++acc].x = A[a].y-A[b].y;
	B[acc].y = A[a].x-A[b].x;
}

void creare(void)
{
	for(int i=1;i<=N;i++)
		for(int j=i+1;j<=N;j++)
			add(i,j);
}

void calc(void)
{
	int j;
	int i = 1;
	int d;
	while(i<acc)
	{
		j = i+1;
		d = 0;
		while(j<=acc && B[i].x/B[i].y == B[j].x/B[j].y)
		{
			d ++;
			j ++;
		}
		nr += (d*(d+1))/2;
		i = j;
	}
	nr += (nr1*(nr1-1))/2;
	nr += (nr2*(nr2-1))/2;
}

long int poz(int li,int ls)
{
	xy a;
	int c;
	int i1 = 0;
	int j1 = -1;
	while(li<ls)
	{
		if(double(B[li].x/B[li].y)>double(B[ls].x/B[ls].y))
		{
			a = B[li];
			B[li] = B[ls];
			B[ls] = a;
			c = -i1;
			i1 = -j1;
			j1 = c;
		}
		li += i1;
		ls += j1;
	}
	return li;
}

void quicksort(int li,int ls)
{
	if(li<ls)
	{
	long int k = poz(li,ls);
	quicksort(li,k-1);
	quicksort(k+1,ls);
	}
}

int main()
{
	FILE *f = fopen("trapez.out","w");
	
	citire();
	creare();
	quicksort(1,acc);
	calc();
	fprintf(f,"%d",nr);
	
	fclose(f);
	return 0;
}