Описание тега pari-gp

Questions about the computer algebra system PARI/GP, including questions about the GP scripting language and the gp interactive shell. Please use the tag [pari] for the PARI C library.
3 ответа

Вызов PARI/GP из Python

Я хотел бы вызвать PARI/GP из Python только для вычисления функции nextprime(n) для разных nс, что я определяю. К сожалению, я не могу установить pari-python для установки, поэтому я подумал, что просто вызову его с помощью командной строки через os…
18 фев '13 в 19:03
0 ответов

Увеличьте точность в PARI/GP с помощью bnfinit

Есть ли способ повысить точность еще больше в PARI/GP, если флаг =2 должен увеличивать точность до тех пор, пока не будет вычислен альфа (второй вектор)? Сколько памяти это потребует? Я настаиваю на вычислении главных генераторов. Конкретная ошибка,…
07 янв '18 в 06:23
1 ответ

Как увеличить размер стека в Pari/Gp, чтобы команда работала

Я работаю с GP и минимальными полиномами следующим образом, работающими на ASUS x75: (19:25) gp > elt=Mod(a*x^3+b*x^2+c*x+d,('x^5-1)/('x-1)) %122 = Mod(a*x^3 + b*x^2 + c*x + d, x^4 + x^3 + x^2 + x + 1) (19:25) gp > (poly=minpoly(elt,x='x)) %12…
14 июн '17 в 02:29
0 ответов

Модульные формы в пари / гп

Я опробую новый пакет модульных форм Pari, и столкнулся с небольшой проблемой, которую не смог решить. Я хочу использовать модульную параметризацию эллиптической кривой E, заданную функцией elltaniyama(E), а затем использовать некоторые новые функци…
29 авг '18 в 15:52
1 ответ

Pari-gp znorder выдает ошибку компиляции

Я запускаю этот простой запрос: for(x=1,10, p=prime(x);a=Mod(100, p);print(a);print(znorder(a))) Это дает мне: Mod(0, 2) *** at top-level: ...od(100,p);print(a);print(znorder(a));) *** ^------------- *** znorder: not an element of (Z/nZ)* in order. …
09 фев '14 в 10:00
1 ответ

Sympy поддержка модульных групп

Sympy поддерживает модульные группы? Я ищу что-то похожее на pari/gp"s Mod функция.
09 окт '18 в 16:42
1 ответ

Как мне настроить gp, чтобы разрешить возобновление кода после того, как я прервал работу с Ctrl-C?

У меня есть ПК с Linux (Ubuntu), где я установил PARI/GP из командной строки (при вводе gp в оболочке он говорит, что программы нет, но дает вам команду, которую вы можете использовать, чтобы установить все это), а затем у меня есть машина Windows, …
19 окт '14 в 19:24
2 ответа

Как мне объявить и вызвать массив в PARI/GP?

Как мне объявить и вызвать массив в PARI/GP? Например, у меня есть следующее в Java: int[] myArray = new int[5]; for(int i = 0; i < 5; i++){ myArray[i] = i + 5; } Как мне сделать то же самое при использовании PARI/GP?
12 ноя '14 в 16:05
1 ответ

PARI: похожая функция, такая как Integer.parseInt()

Я хочу преобразовать текст привет в десятичный ascii в PARI/GP. После этого я объединю значения. Я инициализирую Vecsmall(привет), после этого я запускаю цикл для конкатенации десятичных значений ascii, я хочу использовать это конкатенированное знач…
10 ноя '14 в 14:13
1 ответ

Новое в PARI/GP. Как использовать PARI / GP для с h(p-1/q) mod p

В настоящее время я делаю сценарий для алгоритма цифровой подписи, пока не наткнулся на камень преткновения, который с большим количеством проб и ошибок, я до сих пор не могу пройти через него. Есть 3 целых числа, которые хранятся в p,q,h. Используя…
12 фев '14 в 14:43
1 ответ

Pari/GP: пользовательские функции

Я определил пару функций Arity 1, скажем, func1(-) а также func2(-), Я проверил их и увидел, что они действительно делают то, что должны. Я хочу определить третью функцию, скажем func3(-), что выводит разницу func1(-) а также func2(-), Вот что я дел…
09 фев '13 в 15:15
1 ответ

Неправильный результат в PARI-реализации

Я попытался реализовать алгоритм для расчета мощности башни по модулю m. Ниже процедуры башня должна вычислить 2^3^...^14^15 (мод м), а башня2 должна вычислить 15^14^...^3^2 (мод м). Но для m = 163 tower2 дает неправильный ответ. Я обнаружил, что не…
28 янв '14 в 18:42
1 ответ

Как сделать многострочную вложенность для циклов в PARI/GP?

Как я могу сделать вложенные циклы в PARI/GP, которые занимают несколько строк на каждом уровне? Мне часто приходится делать несколько вещей внутри для циклов, и для удобства чтения я не люблю писать свои циклы в одной строке. Для цикла над одной пе…
03 апр '13 в 17:18
1 ответ

Построение нескольких списков в пари

У меня есть два (пары) списков, которые я хочу построить. Я знаю, что я могу построить каждый отдельно, используя функцию plothraw. Но как я могу нарисовать их на одном и том же изображении, чтобы я получил две кривые разных цветов?
12 янв '18 в 17:09
3 ответа

Как проверить, является ли число целым числом в пари /GP?

Я пытаюсь написать заявление if, как это if(denominator([(i-1)! + 1] / i)-1,print(hi),print(ho)) i может быть любое целое число, например, 10, когда я установил i до 10 это дает эту ошибку. ? [(x-1)! + 1] / x *** this should be an integer: [(x-1)!+1…
22 апр '16 в 01:26
1 ответ

Округление числа в PARI/GP

Я знаю, что команда "\p x" устанавливает точность для всех вычислений, но я ищу что-то немного другое. Учитывая число, вычисленное с высокой точностью, я хочу округлить его до более низкой точности только для одной части моего кода. Причиной этого я…
13 июн '12 в 23:23
2 ответа

Итерация по различным перестановкам вектора в Pari/GP

Я хочу перебрать все различные перестановки вектора. Я пытался сделать это с помощью vecextract() в комбинации с numtoperm() создать вектор перестановок и vecsort(,,,8) удалить эквивалентные перестановки. К сожалению, это плохо масштабируется: макси…
07 янв '19 в 04:06
1 ответ

В pari-gp есть какие-нибудь вещи, чтобы отобразить конечное поле на его некоторое расширение?

Скажем, у меня есть многочлен p(x) над GF(2) и многочлен g(x) над GF(4). Например, gf2 = gf(2, 1); gf4 = gf(2, 2); p2 = Polrev(vector(5, i, random(gf2))); p4 = Polrev(vector(7, i, random(gf4))); p2 * p4 >> *** at top-level: p2*p4 >> *** …
17 фев '16 в 22:04
3 ответа

Есть ли способ рационализировать десятичное число в пари /GP?

Я ищу автоматический способ преобразования десятичного значения в дробь в калькуляторе PARI/GP (gp). Например, я хочу превратить 0,759765625 в 389/512. Я знаю, что могу сделать это вручную, набрав "759765625/10^9", но я бы хотел более простой способ…
03 июн '09 в 21:45
0 ответов

Pari/GP Исключения / разрывные циклы

ellpow(E, P, m) всегда будет выдавать исключение: *** ellpow: impossible inverse modulo: Mod(x, y). *** Break loop: type 'break' to go back to GP где х и у являются целыми числами. Я хочу перехватить значение x, не заканчивая программу, чтобы исполь…