Вычислительная геометрия на плоскости (си)

[Stilet]

Ну и где файл?

[rubakKa]

#include <stdio.h>
main(){
FILE *in,*out;
in=fopen("input.txt","r");
out=fopen("output.txt","w");
 
float s;
int n,min;
char x[50];
char y[50];
fscanf(in,"%d",&n);
for(int i=0;i<n;i++)
{
fscanf(in,"%d",&x[i]);
fscanf(in,"%d",&y[i]);
}
x[n]=x[0];
y[n]=y[0];
 
min=y[0];
for(int i=1;i<n;i++)
if(min>y[i])min=y[i];
 
for(int i=0;i<n+1;i++)
y[i]=y[i]-min;
 
s=0;
for(int i=0;i<n;i++){
s=s+((y[i+1]+y[i])*(x[i+1]-x[i])/2);
}
if(s<0)
s=-s;
fprintf(out,"%1.1f",s);
fclose(in);
fclose(out);
}

Отредактировал.. Почти все нормально: заместо 3.5 3.0.. Где я напутал?

[V0id]

ИМХО, должно быть как-то так...

#include <stdio.h>
#include <math.h>

float round(float x)
{
	float n;
	float fract = modf(x*10, &n);
	if (fract >= 0.5)
		n += 1;
	return n / 10;
}

int main()
{
	int n;
	float s, min, x[500], y[500];

	FILE *in,*out;
	in=fopen("input.txt","rt");
	out=fopen("output.txt","wt");

	if (in && out)
	{
		fscanf(in,"%d",&n);
		for(int i=0;i<n;i++)
			fscanf(in,"%f %f",&x[i], &y[i]);
		x[n]=x[0];
		y[n]=y[0];
		min=y[0];
		for(int i=1;i<n;i++)
			if(min>y[i]) 
				min=y[i];
		for(int i=0;i<n+1;i++)
			y[i]=y[i]-min;
		s=0.0;
		for(int i=0;i<n;i++)
			s=s+((y[i+1]+y[i])*(x[i+1]-x[i])/2);
		s = fabs(s);
		fprintf(out,"%.1f",round(s));
	}
	else
		printf("Can't open files!");

	fclose(in);
	fclose(out);
	return 0;
}

Но в условии задачи написано, что вершин может быть до 50000 штук, а массивы создаются лишь на 500 элементов. Исправляйте.

Но тестировалось лишь для одного простейшего случая, так что насчет работоспособности не уверен.

[rubakKa]

Сообщение от V0id

ИМХО, должно быть как-то так...

cannot convert `float*' to `double*' for argument `2' to `double modf(double, double*)'

Какую функцию выполняет код:

float round(float x)
{
	float n;
	float fract = modf(x*10, &n);
	if (fract >= 0.5)
		n += 1;
	return n / 10;
}

[V0id]

Сообщение от rubakKa

cannot convert `float*' to `double*' for argument `2' to `double modf(double, double*)'

Какую функцию выполняет код:

float round(float x)
{
	float n;
	float fract = modf(x*10, &n);
	if (fract >= 0.5)
		n += 1;
	return n / 10;
}

Округление до десятых.

[rubakKa]

Координаты записываются в целом виде.. Тогда массивы не float, а long. А зачем тут округление?
Можно просто

fprintf(out,"%1.1f",s);

[V0id]

А зачем тут округление?

Цитирую задание: "Его следует округлить до ближайшего числа с одной цифрой после запятой."

Сообщение от rubakKa

Координаты записываются в целом виде.. Тогда массивы не float, а long.

Конечно, но тогда, действительно, зачем тут округление?

Можно просто

fprintf(out,"%1.1f",s);

Это не округление.

[rubakKa]

#include <stdio.h>
main(){
FILE *in,*out;
in=fopen("input.txt","r");
out=fopen("output.txt","w");
 
float s;
int n,min;
long x[50000];
long y[50000];
fscanf(in,"%d",&n);
for(int i=0;i<n;i++)
{
fscanf(in,"%d",&x[i]);
fscanf(in,"%d",&y[i]);
}
x[n]=x[0];
y[n]=y[0];
 
min=y[0];
for(int i=1;i<n;i++)
if(min>y[i])min=y[i];
 
for(int i=0;i<n+1;i++)
y[i]=y[i]-min;
s=0.0;
for(int i=0;i<n;i++){
s=s+(((float)y[i+1]+(float)y[i])*((float)x[i+1]-(float)x[i])/2);
}
if(s<0)
s=-s;
fprintf(out,"%1.1f",s);
fclose(in);
fclose(out);
}

Вот теперь работает!!!
Спасибо V0id и Гром

[rubakKa]

Задача из той же басни:


Группа экспертов обнаружила на территории нежилого массива множество опасных участков, соприкосновение с которыми небезопасно для жизни человека.
В целях безопасности требуется создать защитный периметр в форме выпуклого многоугольника, который бы смог обезопасить проникновение человека в эту зону.
По заданным координатам опасных участков требуется вычислить минимально возможную площадь опасной зоны, которая попадет в защитный периметр.

Входные данные
В первой строке входного файла INPUT.TXT записано натуральное число N – количество опасных участков. В каждой из N последующих строк находятся два числа Xi и Yi – координаты участков, размерами которых можно пренебречь. При этом участки могут повторяться.
Все числа целые, не превосходящие 1000 по абсолютной величине.
Выходные данные
В выходной файл OUTPUT.TXT необходимо вывести одно число – площадь опасной зоны, округленной до целого значения.

[Lyubitel']

Привет всем!
Мне досталась очень хорошая работа, но не могу разобраться в одном куске кода. Помогите, прошу, объясните как именно рассчитывается координаты нормали в if.

Vertex bufNorm;
for (i=0; i<kol_t/2; i++)
{
GetNormal(t[i], bufNorm); // я полагаю, считается при каждой i-ой итерации - нормаль
t[i].N.Set(bufNorm.x, bufNorm.y, bufNorm.z, 1); //могу предположить, что устанавливается нормаль для i-го треугольника

if(t[i].N.x*InsPoint1.x + t[i].N.y*InsPoint1.y + t[i].N.z*InsPoint1.z + (-t[i].N.x*t[i].A.x - t[i].N.y*t[i].A.y - t[i].N.z*t[i].A.z) > 0)
t[i].N.Set(-t[i].N.x, -t[i].N.y, -t[i].N.z, 1);
} //здесь не могу ничего разобрать

InsPoint1.Set(b[n+2].x, b[n+2].y, b[n+2].z, 1);//для основания конуса (это комментарий уже существовал в программе)
Triangle *t;
Vertex N; //думаю, что координаты текущей нормали, но могу ошибаться