Advertisement
metalni

VI Zadaci za vezbanje prv kolokvium - Balls

Jun 30th, 2021
1,349
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 21.05 KB | None | 0 0
  1. import bisect
  2. import sys
  3.  
  4.  
  5. class Problem:
  6.     def __init__(self, initial, goal=None):
  7.         self.initial = initial
  8.         self.goal = goal
  9.  
  10.     def successor(self, state):
  11.         """За дадена состојба, врати речник од парови {акција : состојба}
  12.        достапни од оваа состојба. Ако има многу следбеници, употребете
  13.        итератор кој би ги генерирал следбениците еден по еден, наместо да
  14.        ги генерирате сите одеднаш.
  15.        :param state: дадена состојба
  16.        :return:  речник од парови {акција : состојба} достапни од оваа
  17.                  состојба
  18.        :rtype: dict
  19.        """
  20.         raise NotImplementedError
  21.  
  22.     def actions(self, state):
  23.         """За дадена состојба state, врати листа од сите акции што може да
  24.        се применат над таа состојба
  25.        :param state: дадена состојба
  26.        :return: листа на акции
  27.        :rtype: list
  28.        """
  29.         raise NotImplementedError
  30.  
  31.     def result(self, state, action):
  32.         """За дадена состојба state и акција action, врати ја состојбата
  33.        што се добива со примена на акцијата над состојбата
  34.        :param state: дадена состојба
  35.        :param action: дадена акција
  36.        :return: резултантна состојба
  37.        """
  38.         raise NotImplementedError
  39.  
  40.     def goal_test(self, state):
  41.         """Врати True ако state е целна состојба. Даденава имплементација
  42.        на методот директно ја споредува state со self.goal, како што е
  43.        специфицирана во конструкторот. Имплементирајте го овој метод ако
  44.        проверката со една целна состојба self.goal не е доволна.
  45.        :param state: дадена состојба
  46.        :return: дали дадената состојба е целна состојба
  47.        :rtype: bool
  48.        """
  49.         return state == self.goal
  50.  
  51.     def path_cost(self, c, state1, action, state2):
  52.         """Врати ја цената на решавачкиот пат кој пристигнува во состојбата
  53.        state2 од состојбата state1 преку акцијата action, претпоставувајќи
  54.        дека цената на патот до состојбата state1 е c. Ако проблемот е таков
  55.        што патот не е важен, оваа функција ќе ја разгледува само состојбата
  56.        state2. Ако патот е важен, ќе ја разгледува цената c и можеби и
  57.        state1 и action. Даденава имплементација му доделува цена 1 на секој
  58.        чекор од патот.
  59.        :param c: цена на патот до состојбата state1
  60.        :param state1: дадена моментална состојба
  61.        :param action: акција која треба да се изврши
  62.        :param state2: состојба во која треба да се стигне
  63.        :return: цена на патот по извршување на акцијата
  64.        :rtype: float
  65.        """
  66.         return c + 1
  67.  
  68.     def value(self):
  69.         """За проблеми на оптимизација, секоја состојба си има вредност.
  70.        Hill-climbing и сличните алгоритми се обидуваат да ја максимизираат
  71.        оваа вредност.
  72.        :return: вредност на состојба
  73.        :rtype: float
  74.        """
  75.         raise NotImplementedError
  76.  
  77.  
  78. """
  79. Дефинирање на класата за структурата на јазел од пребарување.
  80. Класата Node не се наследува
  81. """
  82.  
  83.  
  84. class Node:
  85.     def __init__(self, state, parent=None, action=None, path_cost=0):
  86.         """Креирај јазол од пребарувачкото дрво, добиен од parent со примена
  87.        на акцијата action
  88.        :param state: моментална состојба (current state)
  89.        :param parent: родителска состојба (parent state)
  90.        :param action: акција (action)
  91.        :param path_cost: цена на патот (path cost)
  92.        """
  93.         self.state = state
  94.         self.parent = parent
  95.         self.action = action
  96.         self.path_cost = path_cost
  97.         self.depth = 0  # search depth
  98.         if parent:
  99.             self.depth = parent.depth + 1
  100.  
  101.     def __repr__(self):
  102.         return "<Node %s>" % (self.state,)
  103.  
  104.     def __lt__(self, node):
  105.         return self.state < node.state
  106.  
  107.     def expand(self, problem):
  108.         """Излистај ги јазлите достапни во еден чекор од овој јазол.
  109.        :param problem: даден проблем
  110.        :return: листа на достапни јазли во еден чекор
  111.        :rtype: list(Node)
  112.        """
  113.  
  114.         return [self.child_node(problem, action)
  115.                 for action in problem.actions(self.state)]
  116.  
  117.     def child_node(self, problem, action):
  118.         """Дете јазел
  119.        :param problem: даден проблем
  120.        :param action: дадена акција
  121.        :return: достапен јазел според дадената акција
  122.        :rtype: Node
  123.        """
  124.         next_state = problem.result(self.state, action)
  125.         return Node(next_state, self, action,
  126.                     problem.path_cost(self.path_cost, self.state,
  127.                                       action, next_state))
  128.  
  129.     def solution(self):
  130.         """Врати ја секвенцата од акции за да се стигне од коренот до овој јазол.
  131.        :return: секвенцата од акции
  132.        :rtype: list
  133.        """
  134.         return [node.action for node in self.path()[1:]]
  135.  
  136.     def solve(self):
  137.         """Врати ја секвенцата од состојби за да се стигне од коренот до овој јазол.
  138.        :return: листа од состојби
  139.        :rtype: list
  140.        """
  141.         return [node.state for node in self.path()[0:]]
  142.  
  143.     def path(self):
  144.         """Врати ја листата од јазли што го формираат патот од коренот до овој јазол.
  145.        :return: листа од јазли од патот
  146.        :rtype: list(Node)
  147.        """
  148.         x, result = self, []
  149.         while x:
  150.             result.append(x)
  151.             x = x.parent
  152.         result.reverse()
  153.         return result
  154.  
  155.     """Сакаме редицата од јазли кај breadth_first_search или
  156.    astar_search да не содржи состојби - дупликати, па јазлите што
  157.    содржат иста состојба ги третираме како исти. [Проблем: ова може
  158.    да не биде пожелно во други ситуации.]"""
  159.  
  160.     def __eq__(self, other):
  161.         return isinstance(other, Node) and self.state == other.state
  162.  
  163.     def __hash__(self):
  164.         return hash(self.state)
  165.  
  166.  
  167. """
  168. Дефинирање на помошни структури за чување на листата на генерирани, но непроверени јазли
  169. """
  170.  
  171.  
  172. class Queue:
  173.     """Queue е апстрактна класа / интерфејс. Постојат 3 типа:
  174.        Stack(): Last In First Out Queue (стек).
  175.        FIFOQueue(): First In First Out Queue (редица).
  176.        PriorityQueue(order, f): Queue во сортиран редослед (подразбирливо,од најмалиот кон
  177.                                 најголемиот јазол).
  178.    """
  179.  
  180.     def __init__(self):
  181.         raise NotImplementedError
  182.  
  183.     def append(self, item):
  184.         """Додади го елементот item во редицата
  185.        :param item: даден елемент
  186.        :return: None
  187.        """
  188.         raise NotImplementedError
  189.  
  190.     def extend(self, items):
  191.         """Додади ги елементите items во редицата
  192.        :param items: дадени елементи
  193.        :return: None
  194.        """
  195.         raise NotImplementedError
  196.  
  197.     def pop(self):
  198.         """Врати го првиот елемент од редицата
  199.        :return: прв елемент
  200.        """
  201.         raise NotImplementedError
  202.  
  203.     def __len__(self):
  204.         """Врати го бројот на елементи во редицата
  205.        :return: број на елементи во редицата
  206.        :rtype: int
  207.        """
  208.         raise NotImplementedError
  209.  
  210.     def __contains__(self, item):
  211.         """Проверка дали редицата го содржи елементот item
  212.        :param item: даден елемент
  213.        :return: дали queue го содржи item
  214.        :rtype: bool
  215.        """
  216.         raise NotImplementedError
  217.  
  218.  
  219. class Stack(Queue):
  220.     """Last-In-First-Out Queue."""
  221.  
  222.     def __init__(self):
  223.         self.data = []
  224.  
  225.     def append(self, item):
  226.         self.data.append(item)
  227.  
  228.     def extend(self, items):
  229.         self.data.extend(items)
  230.  
  231.     def pop(self):
  232.         return self.data.pop()
  233.  
  234.     def __len__(self):
  235.         return len(self.data)
  236.  
  237.     def __contains__(self, item):
  238.         return item in self.data
  239.  
  240.  
  241. class FIFOQueue(Queue):
  242.     """First-In-First-Out Queue."""
  243.  
  244.     def __init__(self):
  245.         self.data = []
  246.  
  247.     def append(self, item):
  248.         self.data.append(item)
  249.  
  250.     def extend(self, items):
  251.         self.data.extend(items)
  252.  
  253.     def pop(self):
  254.         return self.data.pop(0)
  255.  
  256.     def __len__(self):
  257.         return len(self.data)
  258.  
  259.     def __contains__(self, item):
  260.         return item in self.data
  261.  
  262.  
  263. class PriorityQueue(Queue):
  264.     """Редица во која прво се враќа минималниот (или максималниот) елемент
  265.    (како што е определено со f и order). Оваа структура се користи кај
  266.    информирано пребарување"""
  267.     """"""
  268.  
  269.     def __init__(self, order=min, f=lambda x: x):
  270.         """
  271.        :param order: функција за подредување, ако order е min, се враќа елементот
  272.                      со минимална f(x); ако order е max, тогаш се враќа елементот
  273.                      со максимална f(x).
  274.        :param f: функција f(x)
  275.        """
  276.         assert order in [min, max]
  277.         self.data = []
  278.         self.order = order
  279.         self.f = f
  280.  
  281.     def append(self, item):
  282.         bisect.insort_right(self.data, (self.f(item), item))
  283.  
  284.     def extend(self, items):
  285.         for item in items:
  286.             bisect.insort_right(self.data, (self.f(item), item))
  287.  
  288.     def pop(self):
  289.         if self.order == min:
  290.             return self.data.pop(0)[1]
  291.         return self.data.pop()[1]
  292.  
  293.     def __len__(self):
  294.         return len(self.data)
  295.  
  296.     def __contains__(self, item):
  297.         return any(item == pair[1] for pair in self.data)
  298.  
  299.     def __getitem__(self, key):
  300.         for _, item in self.data:
  301.             if item == key:
  302.                 return item
  303.  
  304.     def __delitem__(self, key):
  305.         for i, (value, item) in enumerate(self.data):
  306.             if item == key:
  307.                 self.data.pop(i)
  308.  
  309.  
  310. def tree_search(problem, fringe):
  311.     """ Пребарувај низ следбениците на даден проблем за да најдеш цел.
  312.    :param problem: даден проблем
  313.    :type problem: Problem
  314.    :param fringe:  празна редица (queue)
  315.    :type fringe: FIFOQueue or Stack or PriorityQueue
  316.    :return: Node or None
  317.    :rtype: Node
  318.    """
  319.     fringe.append(Node(problem.initial))
  320.     while fringe:
  321.         node = fringe.pop()
  322.         print(node.state)
  323.         if problem.goal_test(node.state):
  324.             return node
  325.         fringe.extend(node.expand(problem))
  326.     return None
  327.  
  328.  
  329. def breadth_first_tree_search(problem):
  330.     """Експандирај го прво најплиткиот јазол во пребарувачкото дрво.
  331.    :param problem: даден проблем
  332.    :type problem: Problem
  333.    :return: Node or None
  334.    :rtype: Node
  335.    """
  336.     return tree_search(problem, FIFOQueue())
  337.  
  338.  
  339. def depth_first_tree_search(problem):
  340.     """Експандирај го прво најдлабокиот јазол во пребарувачкото дрво.
  341.    :param problem: даден проблем
  342.    :type problem: Problem
  343.    :return: Node or None
  344.    :rtype: Node
  345.    """
  346.     return tree_search(problem, Stack())
  347.  
  348.  
  349. """
  350. Неинформирано пребарување во рамки на граф
  351. Основната разлика е во тоа што овде не дозволуваме јамки,
  352. т.е. повторување на состојби
  353. """
  354.  
  355.  
  356. def graph_search(problem, fringe):
  357.     """Пребарувај низ следбениците на даден проблем за да најдеш цел.
  358.     Ако до дадена состојба стигнат два пата, употреби го најдобриот пат.
  359.    :param problem: даден проблем
  360.    :type problem: Problem
  361.    :param fringe:  празна редица (queue)
  362.    :type fringe: FIFOQueue or Stack or PriorityQueue
  363.    :return: Node or None
  364.    :rtype: Node
  365.    """
  366.     closed = set()
  367.     fringe.append(Node(problem.initial))
  368.     while fringe:
  369.         node = fringe.pop()
  370.         if problem.goal_test(node.state):
  371.             return node
  372.         if node.state not in closed:
  373.             closed.add(node.state)
  374.             fringe.extend(node.expand(problem))
  375.     return None
  376.  
  377.  
  378. def breadth_first_graph_search(problem):
  379.     """Експандирај го прво најплиткиот јазол во пребарувачкиот граф.
  380.    :param problem: даден проблем
  381.    :type problem: Problem
  382.    :return: Node or None
  383.    :rtype: Node
  384.    """
  385.     return graph_search(problem, FIFOQueue())
  386.  
  387.  
  388. def depth_first_graph_search(problem):
  389.     """Експандирај го прво најдлабокиот јазол во пребарувачкиот граф.
  390.    :param problem: даден проблем
  391.    :type problem: Problem
  392.    :return: Node or None
  393.    :rtype: Node
  394.    """
  395.     return graph_search(problem, Stack())
  396.  
  397.  
  398. def depth_limited_search(problem, limit=50):
  399.     """Експандирај го прво најдлабокиот јазол во пребарувачкиот граф
  400.    со ограничена длабочина.
  401.    :param problem: даден проблем
  402.    :type problem: Problem
  403.    :param limit: лимит за длабочината
  404.    :type limit: int
  405.    :return: Node or None
  406.    :rtype: Node
  407.    """
  408.  
  409.     def recursive_dls(node, problem, limit):
  410.         """Помошна функција за depth limited"""
  411.         cutoff_occurred = False
  412.         if problem.goal_test(node.state):
  413.             return node
  414.         elif node.depth == limit:
  415.             return 'cutoff'
  416.         else:
  417.             for successor in node.expand(problem):
  418.                 result = recursive_dls(successor, problem, limit)
  419.                 if result == 'cutoff':
  420.                     cutoff_occurred = True
  421.                 elif result is not None:
  422.                     return result
  423.         if cutoff_occurred:
  424.             return 'cutoff'
  425.         return None
  426.  
  427.     return recursive_dls(Node(problem.initial), problem, limit)
  428.  
  429.  
  430. def iterative_deepening_search(problem):
  431.     """Експандирај го прво најдлабокиот јазол во пребарувачкиот граф
  432.    со ограничена длабочина, со итеративно зголемување на длабочината.
  433.    :param problem: даден проблем
  434.    :type problem: Problem
  435.    :return: Node or None
  436.    :rtype: Node
  437.    """
  438.     for depth in range(sys.maxsize):
  439.         result = depth_limited_search(problem, depth)
  440.         if result is not 'cutoff':
  441.             return result
  442.  
  443.  
  444. def uniform_cost_search(problem):
  445.     """Експандирај го прво јазолот со најниска цена во пребарувачкиот граф.
  446.    :param problem: даден проблем
  447.    :type problem: Problem
  448.    :return: Node or None
  449.    :rtype: Node
  450.    """
  451.     return graph_search(problem, PriorityQueue(min, lambda a: a.path_cost))
  452.  
  453.  
  454. class Topki(Problem):
  455.  
  456.     def __init__(self, obstacles, N, initial, goal=None):
  457.         super().__init__(initial, goal)
  458.         self.obstacles = obstacles
  459.         self.N = N
  460.  
  461.     def successor(self, state):
  462.         successors = dict()
  463.  
  464.         for t in state:
  465.             x = t[0]
  466.             y = t[1]
  467.  
  468.             # Up - Left
  469.             if x - 2 >= 0 and y + 2 < self.N:
  470.                 tmp_balls = list(state)
  471.                 if (x - 1, y + 1) in state and (x - 2, y + 2) not in state and (x - 2, y + 2) not in self.obstacles:
  472.                     tmp_balls.remove(t)
  473.                     tmp_balls.remove((x - 1, y + 1))
  474.                     tmp_balls.append((x - 2, y + 2))
  475.                     successors['Gore Levo: (x=' + str(x) + ',y=' + str(y) + ')'] = tuple(tmp_balls)
  476.  
  477.             # Up - Right
  478.             if x + 2 < self.N and y + 2 < self.N:
  479.                 tmp_balls = list(state)
  480.                 if (x + 1, y + 1) in state and (x + 2, y + 2) not in state and (x + 2, y + 2) not in self.obstacles:
  481.                     tmp_balls.remove(t)
  482.                     tmp_balls.remove((x + 1, y + 1))
  483.                     tmp_balls.append((x + 2, y + 2))
  484.                     successors['Gore Desno: (x=' + str(x) + ',y=' + str(y) + ')'] = tuple(tmp_balls)
  485.  
  486.             # Down - Left
  487.             if x - 2 >= 0 and y - 2 >= 0:
  488.                 tmp_balls = list(state)
  489.                 if (x - 1, y - 1) in state and (x - 2, y - 2) not in state and (x - 2, y - 2) not in self.obstacles:
  490.                     tmp_balls.remove(t)
  491.                     tmp_balls.remove((x - 1, y - 1))
  492.                     tmp_balls.append((x - 2, y - 2))
  493.                     successors['Dolu Levo: (x=' + str(x) + ',y=' + str(y) + ')'] = tuple(tmp_balls)
  494.  
  495.             # Down - Right
  496.             if x + 2 < self.N and y - 2 >= 0:
  497.                 tmp_balls = list(state)
  498.                 if (x + 1, y - 1) in state and (x + 2, y - 2) not in state and (x + 2, y - 2) not in self.obstacles:
  499.                     tmp_balls.remove(t)
  500.                     tmp_balls.remove((x + 1, y - 1))
  501.                     tmp_balls.append((x + 2, y - 2))
  502.                     successors['Dolu Desno: (x=' + str(x) + ',y=' + str(y) + ')'] = tuple(tmp_balls)
  503.  
  504.  
  505.             # Right
  506.             if x + 2 < self.N:
  507.                 tmp_balls = list(state)
  508.                 if (x + 1, y) in state and (x + 2, y) not in state and (x + 2, y) not in obstacles:
  509.                     tmp_balls.remove(t)
  510.                     tmp_balls.remove((x + 1, y))
  511.                     tmp_balls.append((x + 2, y))
  512.                     successors['Desno: (x=' + str(x) + ',y=' + str(y) + ')'] = tuple(tmp_balls)
  513.  
  514.             # Left
  515.             if x - 2 >= 0:
  516.                 tmp_balls = list(state)
  517.                 if (x - 1, y) in state and (x - 2, y) not in state and (x - 2, y) not in obstacles:
  518.                     tmp_balls.remove(t)
  519.                     tmp_balls.remove((x - 1, y))
  520.                     tmp_balls.append((x - 2, y))
  521.                     successors['Levo: (x=' + str(x) + ',y=' + str(y) + ')'] = tuple(tmp_balls)
  522.  
  523.         return successors
  524.  
  525.     def actions(self, state):
  526.         return self.successor(state).keys()
  527.  
  528.     def result(self, state, action):
  529.         return self.successor(state)[action]
  530.  
  531.     def goal_test(self, state):
  532.         if len(state) == 1:
  533.             if state[0][0] == self.N // 2 and state[0][1] == self.N - 1:
  534.                 return True
  535.         return False
  536.  
  537.  
  538. if __name__ == "__main__":
  539.     dimenzii = int(input())
  540.     num_balls = int(input())
  541.     balls = list()
  542.     obstacles = list()
  543.  
  544.     for i in range(num_balls):
  545.         balls.append(tuple([int(j) for j in input().split(",")]))
  546.  
  547.     num_obstacles = int(input())
  548.  
  549.     for i in range(num_obstacles):
  550.         obstacles.append(tuple([int(j) for j in input().split(",")]))
  551.  
  552.     topki = Topki(tuple(obstacles), dimenzii, tuple(balls))
  553.     res = breadth_first_graph_search(topki)
  554.     print(res.solution())
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement