#include <iostream>
#include <pari/pari.h>

int main (int argc, char **argv){
    // pari_init нужна обязательно!!!
    // первый аргумент - это размер стека, второй - число, определяющее таблицу
    // простых чисел (в ней будут простые числа, меньшие указанного (вроде бы))
    // подробнее на стр.47
    pari_init(4096*4096*2, 10000);

    // mkmat2 создает матрицу из двух столбцов
    // mkcol2 создает столбец(вектор) из двух элементов
    // подробнее на стр.159
    // gen_0 это определенное в библиотеке название для 0
    // gen_1 аналогично для 1
    // подробнее на стр.12
    GEN a = mkmat2(mkcol2(gen_0, gen_1), mkcol2(gen_1, gen_0));

    // zeromatcopy(m, n) создает изменяемую нулевую матрицу m на n
    GEN b = zeromatcopy(2, 2);

    b = qfgaussred(a);

    // output(x) печатает объект x (число, вектор, вообще любой объект PARI)
    // подробнее на стр.36
    output(b);

    // pari_close как и pari_init нужна обязательно
    pari_close();
}