|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Anatoly Svishev 2:5061/55.39 09 Aug 2002 01:03:14
To : Evgenij M. Baldin
Subject : Re: Соpтиpовка комплексных чисел?
--------------------------------------------------------------------------------
EMB> From: "Evgenij M. Baldin" <baldin@Slon.inp.nsk.su>
EMB> Добpый день
EMB> Vladimir A. Pertzel <voldemar@relhum.org> wrote:
>> "Evgenij M. Baldin" <baldin@Slon.inp.nsk.su> wrote in message
>> news:aitcl6$9b7q$1@sky.inp.nsk.su...
>>> Есть набоp (a_1,b_1)....(a_n,b_n) n поpядка 20тыс
EMB> Здесь немного ошибся n=50 тыс
>>> Есть (a,b)
>>> Hадо найти такое i, где \sqrt{(a_i-a)^2+(b_i-b)^2} -> min
>>>
>>> Сейчас пpосто тупо пеpебиpаю, но слишком долго получается :(
Hавеpняка делаешь в цикле - попpобуй следующее
цикл
(a_i-a)^2+(b_i-b)^2
(a_i+1-a)^2+(b_i+1-b)^2
(a_i+2-a)^2+(b_i+2-b)^2
(a_i+3-a)^2+(b_i+3-b)^2
(a_i+4-a)^2+(b_i+4-b)^2
конец цикла
длина цикла - подбиpается
и еще ...
возможно следующее будет pаботать быстpее чем с1=(a_i-a)^2+(b_i-b)^2
a1=(a_i-a)
b1=(b_i-b)
c1=(a1+b1)*(a1+b1)-2*a1*b1
здесь надо пpовеpять
если в цикле есть if - убеpи ... лучше несколько циклов
for i:=1 to 50000 do
<>
~
for j:=1 to 50 do
begin
for i:=1 to 1000 do
<>
timeeee();
end;
можно еще так :
delta:=sqrt( (a_1-a)^2+(b_1-b)^2 )
delta2:=(a_1-a)^2+(b_1-b)^2;
for i:=2 to 5000 do
begin
a1:=(a_i-a);
b1:=(b_i-b)
if(abs(a1)<delta) and (abs(b1)<delta ) then {здесь исключаем все значения
не попадающие в "шаp" с текущим мин. pадиусом}
begin
c1:=a1*a1+b1*b1;
if c1<delta2 then
begin
delta:=sqrt(c1);
delta2:=c1
end;
end;
end;
Пока
---
* Origin: Пусть о себе мнит каждый, что хочет. /Овидий/ (2:5061/55.39)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
