интерактивный курс: введение в программирование на PARI/GP

[Soul Friend]

Как отсортировать (от меньшего значения к большему) в коде вывод последовательностей функции от двух переменных?

for(i=1, 10, for(j=1, 5, print(i^2+j!)))

[dmd]

Soul Friend

S=Set(); for(i=1, 10, for(j=1, 5, S=setunion(S, [i^2+j!]))); print(S)

[Soul Friend]

Спасибо, итого у вас получилось 777 сообщений.
Только что сам решил как сделать, и хотел написать ответ, но вы опередили. Но у меня немного по другому:

а=vector(50);
for(n=1, 50, for(i=1, 10, for(j=1, 5, a[n]=i^2+j! ; n++)));
vecsort(a)

Но в вашем случае можно не заботиться о размере вектора, и повторений нет.

[Dmitriy40]

Я-то думал надо отсортировать матрицу, а не вектор ... С вектором разумеется vecsort решает задачу.
У vecsort есть третий параметр, им можно например исключить дубли.
Чтобы не задавать размер вектора в начале можно использовать concat:

v=[];for(i=1, 10, for(j=1, 5, v=concat(v,[i^2+j!])));v=vecsort(v)

Soul Friend в сообщении #1329527 писал(а):

и повторений нет.

Таки есть, числа 10 и 124 повторяются дважды.

[maxal]

Soul Friend в сообщении #1329410 писал(а):

Как отсортировать (от меньшего значения к большему) в коде вывод последовательностей функции от двух переменных?

for(i=1, 10, for(j=1, 5, print(i^2+j!)))

как вариант:

vecsort(concat(vector(10,i,vector(5,j,i^2+j!))))

[Soul Friend]

Спасибо, пригодятся все варианты.

Dmitriy40 в сообщении #1329592 писал(а):

Таки есть, числа 10 и 124 повторяются дважды.

У меня все code-ы записаны, у dmd нет повторений.
Для удобства вывода в В-файл испорчу красивый код от maxal:

b=vecsort(concat(vector(10,i,vector(5,j,i^2+j!))));
for(i=1, 50, print(i, "  ", b[i]))

Сократить бы это всё в одну строчку.

[kotenok gav]

Soul Friend в сообщении #1329617 писал(а):

Сократить бы это всё в одну строчку.

for(i=1, 50, print(i, "  ", vecsort(concat(vector(10,i,vector(5,j,i^2+j!))))[i]))

[maxal]

kotenok gav, так можно, но неэфективно. Вектор в цикле будет каждый раз пересчитываться.

Вот ещё вариант кода:

vecsort(concat(Vec(matrix(10,5,i,j,i^2+j!))))

[Dmitriy40]

Был бы forall/foreach, выполняющийся над вектором для каждого элемента, получилось бы в строчку. Но я ничего подобного в хелпе не нашёл. Если только свою функцию задать ...
PS. Вместо числа 50 в цикле for удобно использовать конструкцию #b, выдающую размер вектора.

[maxal]

Dmitriy40, см. apply()

[Dmitriy40]

Спасибо!!
Вот тогда решение в одну строчку (воспользуюсь Вашим трюком с матрицей как более наглядным), даже два варианта:

[print1(x," ")|x<-vecsort(concat(Vec(matrix(10,5,i,j,i^2+j!))))];
apply(x->print1(x," "),vecsort(concat(Vec(matrix(10,5,i,j,i^2+j!)))));

Уф, таки получилось и с apply(), но с этими стрелочками туда-сюда фиг разберёшься.

[maxal]

Dmitriy40, если построчно и с номерами, то

k=0; [print(k++," ",x)|x<-vecsort(concat(Vec(matrix(10,5,i,j,i^2+j!))))];

[Soul Friend]

Онлайн калькулятор PARI/GP :
https://pari.math.u-bordeaux.fr/gp.html
Там и справочник есть, только на английском.

[wrest]

Как проверить что натуральное $n=a^b; a,b>1$ т.е. $n$ является степенью (pure power, perfect power)?

[Dmitriy40]

wrest

? ?ispower
ispower(x,{k},{&n}): if k > 0 is given, return true (1) if x is a k-th power, false (0) if not. If k is omitted, return the maximal k >= 2 such that x = n^k is a perfect power, or 0 if no such k exist. If n is present, and the function returns a non-zero result, set n to the k-th root of x.