usun(X,[X|Reszta],Reszta).
usun(X,[Y|Ogon],[Y|Reszta]) :-
    usun(X,Ogon,Reszta).

wstaw(X,L,Duza) :-
    usun(X,Duza,L).

ciag_arytm([], 0, _, _).
ciag_arytm([S|L1], N, S, K):-
    N1 is N - 1,S1 is S + K,ciag_arytm(L1, N1, S1, K).

permutacja([],[]).
permutacja([X|L],P) :-
    permutacja(L,L1),
    wstaw(X,L1,P).

iloczyn([X],X).
iloczyn([G|O],W):-
    iloczyn(O,W1),W is W1*G.

rodzice(jan,ewa,[adam],[maria]).
rodzice(adam,halina,[norbert],[monika]).
rodzice(piotr,maria,[krzysztof],[]).
ojciec(edward,[],halina).

dziadek(X,Y):-
    ojciec(X,_,Z),rodzice(_,Z,_,[Y]).
dziadek(X,Y):-
    ojciec(X,_,Z),rodzice(_,Z,[Y],_).
dziadek(X,Y):-
    ojciec(X,Z,_),rodzice(Z,_,[Y],_).
dziadek(X,Y):-
    ojciec(X,Z,_),rodzice(Z,_,_,[Y]).

dziadek(X,Y):-
    rodzice(X,_,[Z],_),rodzice(Z,_,[Y],_).
dziadek(X,Y):-
    rodzice(X,_,[Z],_),rodzice(Z,_,_,[Y]).
dziadek(X,Y):-
    rodzice(X,_,_,[Z]),rodzice(_,Z,[Y],_).
dziadek(X,Y):-
    rodzice(X,_,_,[Z]),rodzice(_,Z,_,[Y]).

liczba_wnukow(X,N):-
    findall(1,dziadek(X,_),L),length(L,N).