Как я могу определить порядок простоты в pari/gp?

Я глубоко интересуюсь теорией чисел и хочу проверить некоторые из своих идей в pari/gp, но я не знаком с этим программным обеспечением. В частности, я хочу определить функцию «простого порядка», которая отображает целое число n в то, что я называю его порядком простоты, который равен 0 тогда и только тогда, когда n составное и равно наименьшему целому числу k, такому что k-я итерация Функция подсчета простых чисел, оцениваемая при n, в противном случае является составной.

Как определить такую ​​функцию в pari/gp?


person Sylvain JULIEN    schedule 16.09.2017    source источник


Ответы (1)


Пожалуйста, просмотрите это:

primeorder(x) = {
  if(!isprime(x), return(0));

  my(k=1, p=primepi(x));
  while(isprime(p), p=primepi(p); k++);
  return(k);
}

Я настоятельно рекомендую прочитать руководство по PARI/GP.

person Piotr Semenov    schedule 17.09.2017
comment
Кажется, это уже известно в OEIS: oeis.org/ - person Sylvain JULIEN; 17.09.2017