disp("1")
disp("a")

A = -19:34
B = []

for p = -40:70
    for q = 1:3
        B($+1) = p/q;
    end
end
disp(A)
disp(B)
disp("b")

uniaoAB = union(A,B);
interAB = intersect(A,B);

C = setdiff(uniaoAB,interAB);

disp(C)
disp("2")

n=1;
while abs(((1+(5/n))^n)-exp(5))>= 10^(-1) then
    n=n+1;
end
disp(n);


disp("3")
disp("a")
disp("i.")

sheet = readxls("Material/Data.xls")

MS = sheet(1).value;


disp("ii.")

function output = issimetrica(A)
    [a,b] = size(A);
    output = %t;

    for i = 1:a
        for j = 1:b
            if A(i,j)<>A(j,i) then
                output = %f;
            end
        end
    end
endfunction

disp(issimetrica(MS))
disp("b")

function output = grau_out(Adj,vert)
    [a,b] = size(Adj);
    output = 0;
    for i = 1:b

        if Adj(vert,i)<>0 then
            output = output+1;
        end
    end
endfunction


[a,b] = size(MS);
for i = 1:b
    lista($+1) = grau_out(MS,i);
end
disp(lista);

disp("c")

T = zeros(a, b)
for i = 3:6
    T = T + MS^i
end
disp(T);
caminhos = T(7,4);
disp("O número de caminhos de comprimento 3 do vertice 7 para o 4 é "+string(caminhos));

disp("d");

exec("Material/Algoritmo_1.sce");

P = Warshall(MS)
disp("Não é fortemente conexo pois P9 contém zeros.")
Uni = bool2s(P+P');
disp(Uni)
disp("Não é unilateralmente conexo pois P9+P9´ contém zeros fora da diagonal principal.")

A1 = MS+MS';
P2 = Warshall(A1);
disp("É fracamente conexo pois a matriz P2 = (MS+MS´)^1+(MS+MS´)^2+(MS+MS´)^3... contém todos os elementos nao nulos")

disp("e)")
pesos = sheet(2).value;
exec("Material/Algoritmo_Warshall_MIN.sce");
Warshall_MIN(pesos)

disp("i.")
disp("O peso mínimo do caminho minimo de 9 para 3 é de 10. Matriz Q8(9,3)")
disp("ii.")
disp("O peso mínimo do caminho minimo de 9 para 3 é de 10. Matriz Q8(9,3). Este caminho passa nos vértices: 9-8-2-1-3 (M8(9,3)), ou seja, também passa no vértice 2")
disp("f)")
function L=Algoritmo2(A,vi)

    [u,v]=size(A)

    Estado=ones(1,u)

    conjunto=list(vi);
    Estado(vi) = 2;


    Vvisitado=list();

    while length(conjunto)<>0

    Vvisitado($+1)=conjunto(1)

    Estado(conjunto(1))=3;

    for j=1:u
          if A(conjunto(1),j) == 1 & Estado(j)==1 then
              conjunto($+1)=j;
              Estado(j)=2;
          end
    end

    conjunto(1)=null()

    end

    L=string(Vvisitado(1))
    for i=2:length(Vvisitado)
         L=L+" -> "+string(Vvisitado(i))
    end

endfunction

disp(Algoritmo2(MS,5))

disp("4")
disp("a")

Arv = [54 32 66 0 38 0 85 0 0 35 0 0 0 80 0]

disp("b")
exec("Material/Algoritmo_3.sce");
exec("Material/tree_plot.sce");
valores = sheet(3).value;


for i = 1:size(valores,2)
    [pos,encontrado]=localiza(Arv,valores(i))
    if encontrado<>%t then
        Arv(pos) = valores(i);
    end
end


tree_plot(Arv)