:- use_module(library(mcintyre)).

:- mc.

:- begin_lpad.

% The folllowing programs represents a
% Barabási Albert preferential attachment model generated
% with python library networkx
% >>> import networkx as nx
% >>> G= nx.barabasi_albert_graph(40,10)

% ?- mc_gibbs_sample(path(0,39),path(0,19),1000,P,[mix(100),block(1)]),
% P = 0.406

connected(A,B):-
    edge(A,B).
connected(A,B):-
    edge(B,A).

path(Start, End) :-
    path(Start, End, [Start]).
path(End, End, _).
path(Start, End, Visited) :-
    connected(Start, Next),
    \+ memberchk(Next, Visited),
    path(Next, End,[Next|Visited]).

edge(0,10):0.1.
edge(0,11):0.1.
edge(0,12):0.1.
edge(0,13):0.1.
edge(0,14):0.1.
edge(0,15):0.1.
edge(0,18):0.1.
edge(0,19):0.1.
edge(0,22):0.1.
edge(0,26):0.1.
edge(0,27):0.1.
edge(0,37):0.1.
edge(0,38):0.1.
edge(1,10):0.1.
edge(1,11):0.1.
edge(1,12):0.1.
edge(1,13):0.1.
edge(1,14):0.1.
edge(1,16):0.1.
edge(1,18):0.1.
edge(1,31):0.1.
edge(1,35):0.1.
edge(2,10):0.1.
edge(2,11):0.1.
edge(2,13):0.1.
edge(2,15):0.1.
edge(2,20):0.1.
edge(2,22):0.1.
edge(2,25):0.1.
edge(2,26):0.1.
edge(2,28):0.1.
edge(2,31):0.1.
edge(2,35):0.1.
edge(3,10):0.1.
edge(3,11):0.1.
edge(3,12):0.1.
edge(3,14):0.1.
edge(3,15):0.1.
edge(3,17):0.1.
edge(3,24):0.1.
edge(3,30):0.1.
edge(3,33):0.1.
edge(3,34):0.1.
edge(3,35):0.1.
edge(3,36):0.1.
edge(4,10):0.1.
edge(4,11):0.1.
edge(4,12):0.1.
edge(4,13):0.1.
edge(4,14):0.1.
edge(4,15):0.1.
edge(4,16):0.1.
edge(4,17):0.1.
edge(4,18):0.1.
edge(4,19):0.1.
edge(4,22):0.1.
edge(4,27):0.1.
edge(4,35):0.1.
edge(4,39):0.1.
edge(5,10):0.1.
edge(5,12):0.1.
edge(5,13):0.1.
edge(5,15):0.1.
edge(5,16):0.1.
edge(5,17):0.1.
edge(5,18):0.1.
edge(5,20):0.1.
edge(5,21):0.1.
edge(5,26):0.1.
edge(5,28):0.1.
edge(5,30):0.1.
edge(5,33):0.1.
edge(6,10):0.1.
edge(6,11):0.1.
edge(6,12):0.1.
edge(6,13):0.1.
edge(6,14):0.1.
edge(6,16):0.1.
edge(6,17):0.1.
edge(6,19):0.1.
edge(6,22):0.1.
edge(6,23):0.1.
edge(6,24):0.1.
edge(6,32):0.1.
edge(6,35):0.1.
edge(6,36):0.1.
edge(6,37):0.1.
edge(6,39):0.1.
edge(7,10):0.1.
edge(7,11):0.1.
edge(7,12):0.1.
edge(7,16):0.1.
edge(7,20):0.1.
edge(7,25):0.1.
edge(7,26):0.1.
edge(8,10):0.1.
edge(8,11):0.1.
edge(8,13):0.1.
edge(8,15):0.1.
edge(8,16):0.1.
edge(8,17):0.1.
edge(8,23):0.1.
edge(9,10):0.1.
edge(9,11):0.1.
edge(9,12):0.1.
edge(9,14):0.1.
edge(9,15):0.1.
edge(9,18):0.1.
edge(9,19):0.1.
edge(9,20):0.1.
edge(9,22):0.1.
edge(9,23):0.1.
edge(9,24):0.1.
edge(9,27):0.1.
edge(9,29):0.1.
edge(9,30):0.1.
edge(9,33):0.1.
edge(9,36):0.1.
edge(9,38):0.1.
edge(9,39):0.1.
edge(10,11):0.1.
edge(10,12):0.1.
edge(10,13):0.1.
edge(10,14):0.1.
edge(10,15):0.1.
edge(10,16):0.1.
edge(10,17):0.1.
edge(10,18):0.1.
edge(10,19):0.1.
edge(10,21):0.1.
edge(10,23):0.1.
edge(10,29):0.1.
edge(10,30):0.1.
edge(10,33):0.1.
edge(10,34):0.1.
edge(10,36):0.1.
edge(10,37):0.1.
edge(10,38):0.1.
edge(11,12):0.1.
edge(11,13):0.1.
edge(11,14):0.1.
edge(11,16):0.1.
edge(11,17):0.1.
edge(11,18):0.1.
edge(11,19):0.1.
edge(11,20):0.1.
edge(11,21):0.1.
edge(11,22):0.1.
edge(11,23):0.1.
edge(11,24):0.1.
edge(11,25):0.1.
edge(11,26):0.1.
edge(11,27):0.1.
edge(11,29):0.1.
edge(11,30):0.1.
edge(11,31):0.1.
edge(11,32):0.1.
edge(11,33):0.1.
edge(12,13):0.1.
edge(12,14):0.1.
edge(12,15):0.1.
edge(12,16):0.1.
edge(12,17):0.1.
edge(12,20):0.1.
edge(12,23):0.1.
edge(12,24):0.1.
edge(12,25):0.1.
edge(12,30):0.1.
edge(12,31):0.1.
edge(12,34):0.1.
edge(12,36):0.1.
edge(12,37):0.1.
edge(12,39):0.1.
edge(13,14):0.1.
edge(13,15):0.1.
edge(13,17):0.1.
edge(13,19):0.1.
edge(13,21):0.1.
edge(13,24):0.1.
edge(13,25):0.1.
edge(13,29):0.1.
edge(13,36):0.1.
edge(13,39):0.1.
edge(14,18):0.1.
edge(14,19):0.1.
edge(14,20):0.1.
edge(14,23):0.1.
edge(14,25):0.1.
edge(14,26):0.1.
edge(14,29):0.1.
edge(14,32):0.1.
edge(14,33):0.1.
edge(14,35):0.1.
edge(14,37):0.1.
edge(14,38):0.1.
edge(15,16):0.1.
edge(15,18):0.1.
edge(15,21):0.1.
edge(15,24):0.1.
edge(15,25):0.1.
edge(15,27):0.1.
edge(15,32):0.1.
edge(15,34):0.1.
edge(15,36):0.1.
edge(15,38):0.1.
edge(16,17):0.1.
edge(16,18):0.1.
edge(16,19):0.1.
edge(16,20):0.1.
edge(16,21):0.1.
edge(16,22):0.1.
edge(16,23):0.1.
edge(16,27):0.1.
edge(16,28):0.1.
edge(16,31):0.1.
edge(16,34):0.1.
edge(17,19):0.1.
edge(17,21):0.1.
edge(17,22):0.1.
edge(17,24):0.1.
edge(17,26):0.1.
edge(17,30):0.1.
edge(17,32):0.1.
edge(17,38):0.1.
edge(18,20):0.1.
edge(18,21):0.1.
edge(18,22):0.1.
edge(18,27):0.1.
edge(18,28):0.1.
edge(18,30):0.1.
edge(18,32):0.1.
edge(18,33):0.1.
edge(18,35):0.1.
edge(18,39):0.1.
edge(19,20):0.1.
edge(19,21):0.1.
edge(19,24):0.1.
edge(19,28):0.1.
edge(19,29):0.1.
edge(19,37):0.1.
edge(19,38):0.1.
edge(20,21):0.1.
edge(20,23):0.1.
edge(20,26):0.1.
edge(20,27):0.1.
edge(20,28):0.1.
edge(20,29):0.1.
edge(20,30):0.1.
edge(20,31):0.1.
edge(20,32):0.1.
edge(20,37):0.1.
edge(21,22):0.1.
edge(21,23):0.1.
edge(21,25):0.1.
edge(21,27):0.1.
edge(21,34):0.1.
edge(21,35):0.1.
edge(22,24):0.1.
edge(22,25):0.1.
edge(22,28):0.1.
edge(22,29):0.1.
edge(22,30):0.1.
edge(22,34):0.1.
edge(23,25):0.1.
edge(23,26):0.1.
edge(23,31):0.1.
edge(23,34):0.1.
edge(24,29):0.1.
edge(24,33):0.1.
edge(24,37):0.1.
edge(24,39):0.1.
edge(25,26):0.1.
edge(25,27):0.1.
edge(25,28):0.1.
edge(25,29):0.1.
edge(25,38):0.1.
edge(26,28):0.1.
edge(26,31):0.1.
edge(26,32):0.1.
edge(27,28):0.1.
edge(27,31):0.1.
edge(27,32):0.1.
edge(27,39):0.1.
edge(28,31):0.1.
edge(28,34):0.1.
edge(28,35):0.1.
edge(28,36):0.1.
edge(28,37):0.1.
edge(29,32):0.1.
edge(29,33):0.1.
edge(29,36):0.1.
edge(30,37):0.1.
edge(32,33):0.1.
edge(32,35):0.1.
edge(32,36):0.1.
edge(33,34):0.1.
edge(33,38):0.1.
edge(33,39):0.1.
edge(37,38):0.1.
edge(37,39):0.1.

:- end_lpad.