Реализации алгоритмов/Решето Аткина: различия между версиями
Содержимое удалено Содержимое добавлено
м →Rust |
|||
Строка 252:
// Просеивание:
// n = 4x² + y²
n = (x2 <<
if n <= limit as u64 && (n %
item = sieve[n as usize];
sieve.set(n as usize, !item);
Строка 259:
// n = 3x² + y²
n -= x2;
if n <= limit as u64 && n %
item = sieve[n as usize];
sieve.set(n as usize, !item);
Строка 265:
// n = 3x² - y² (при x > y)
if x2 > y2 {
n -= y2 <<
if n <= limit as u64 && n %
item = sieve[n as usize];
sieve.set(n as usize, !item);
Строка 272:
}
y2 += dy2;
dy2 +=
}
x2 += dx2;
dx2 +=
}
}
Строка 282:
let mut r = 5u32;
let mut r2 = (r * r) as u64;
let mut dr2 = (r <<
while r2 < limit as u64 {
if sieve[r as usize] {
Строка 291:
}
}
r +=
r2 += dr2;
dr2 +=
}
}
// 2 и 3 — заведомо простые
if limit >
sieve.set(2, true);
}
if limit >
sieve.set(3, true);
}
| |||