Главная страница
    Top.Mail.Ru    Яндекс.Метрика
Форум: "Основная";
Текущий архив: 2002.02.07;
Скачать: [xml.tar.bz2];

Вниз

Люди помогите !!   Найти похожие ветки 

 
Masa   (2002-01-22 12:06) [0]

Нужно найти простые числа их множества.
Пытаюсь найти вот таким способом. Не выходит.
Что здесь не так, помогите !!

Read(N);
MN := [2..N];
PR := [];

FOR k:=2 to N do
if K in MN Then
begin
PR := PR+[K];
for M:=1 TO N DIV K do
MN := MN -[k*n];
end;


 
MBo   (2002-01-22 12:51) [1]

исправь
MN := MN -[k*M];


 
Masa   (2002-01-22 13:04) [2]

>MBo
Так не выводит вообще ничего.


 
MBo   (2002-01-22 13:10) [3]

procedure TForm1.Button1Click(Sender: TObject);
var n,k,m:integer;
mn,pr:set of byte;
begin
N:=20;
MN := [2..N];
PR := [];

FOR k:=2 to N do
if K in MN Then
begin
PR := PR+[K];
for M:=1 TO N DIV K do
MN := MN -[k*m];
end;
for m:=2 to n do
if m in pr then memo1.lines.add(inttostr(m));
end;


 
Digitman   (2002-01-22 13:11) [4]

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

если параметр IsRange = True, то ф-ция возвращает все простые числа в диапазоне (1..значение_параметра_Count) и параметр Count возвратит общее число простых чисел в результирующем массиве, принадлежащих заданному диапазону (макс.значение Count = 536870910)

если параметр IsRange = False, ф-ция возвращает первые Count-простых чисел, начиная с 1 (макс.значение Count = 536870910)

результат ф-ции = nil, если - ошибка (резуль.массив не сформирован). Если <> nil, то после использования его необходимо освободить



type

TLongWordArray = array[0..$1FFFFFFE] of LongWord;
PLongWordArray = ^TLongWordArray;

TByteBoolArray = array[0..$1FFFFFFE] of ByteBool;
PByteBoolArray = ^TByteBoolArray;

implementation

function Sieving(var Count: LongWord; IsRange: Boolean): PLongWordArray;
var
i,j,k,n: LongWord;
Sieve: PByteBoolArray;
SieveLimit: LongWord;
RangeFactor: Double;
begin
if Count < 1 then
Result:= nil
else begin
if not IsRange then begin
SieveLimit:= Round(Count * RangeFactor);
GetMem(Sieve, SieveLimit);
try
GetMem(Result, Count * SizeOf(LongWord));
try
FillChar(Sieve^, SieveLimit, True);
j:= 0;
for i:= 0 to Pred(Count) do begin
while not Sieve[j] do Inc(j);
k:= j;
Inc(j);
n:= j;
Result[i]:= n;
if n > 1 then while SieveLimit - k > n do begin
Inc(k, n);
Sieve[k]:= False;
end;
end;
except
FreeMem(Result);
raise;
end;
finally
FreeMem(Sieve);
end;
end else begin
if Count = 1 then
SieveLimit:= Count
else
//


 
Digitman   (2002-01-22 13:14) [5]

P.S.
Галиматья в комментарии в коде так выглядит в оригинале :
//ф-ция распределения простых чисел близка к Y(X) = X / Ln(X)



Страницы: 1 вся ветка

Форум: "Основная";
Текущий архив: 2002.02.07;
Скачать: [xml.tar.bz2];

Наверх





Память: 0.45 MB
Время: 0.007 c
1-18796
Tadeush
2002-01-24 14:57
2002.02.07
Формы


3-18783
Val
2002-01-14 11:17
2002.02.07
Коннект к базе


14-18973
V-o-y-a-g-e-r
2001-12-19 15:03
2002.02.07
Компонент Вольтметр !!!


1-18910
Dtr
2002-01-24 10:01
2002.02.07
PopUpMenu


1-18831
vopros
2002-01-21 16:03
2002.02.07
Календарь компонента





Afrikaans Albanian Arabic Armenian Azerbaijani Basque Belarusian Bulgarian Catalan Chinese (Simplified) Chinese (Traditional) Croatian Czech Danish Dutch English Estonian Filipino Finnish French
Galician Georgian German Greek Haitian Creole Hebrew Hindi Hungarian Icelandic Indonesian Irish Italian Japanese Korean Latvian Lithuanian Macedonian Malay Maltese Norwegian
Persian Polish Portuguese Romanian Russian Serbian Slovak Slovenian Spanish Swahili Swedish Thai Turkish Ukrainian Urdu Vietnamese Welsh Yiddish Bengali Bosnian
Cebuano Esperanto Gujarati Hausa Hmong Igbo Javanese Kannada Khmer Lao Latin Maori Marathi Mongolian Nepali Punjabi Somali Tamil Telugu Yoruba
Zulu
Английский Французский Немецкий Итальянский Португальский Русский Испанский