Python: get all node descendants in a tree












0












$begingroup$


I have a CSV holding a "flat" table of a tree's edges (NOT binary, but a node cannot have two parents), ~1M edges:



node_id parent_id
1 0
2 1
3 1
4 2
...


The nodes are sorted in a way that a parent_id must always come before any of its children, so a parent_id will always be lower than node_id.



I wish, for each node_id, to get the set of all ancestor nodes (including itself, propagated until root which is node 0 here), and a set of all descendant nodes (including itself, propagated until leaves), and speed is crucial.



Currently what I do at high level:




  1. Read the CSV in pandas, call it nodes_df

  2. Iterate once through nodes_df to get node_ancestors, a {node_id: set(ancestors)} dict adding for each node's ancestors itself and its parent's ancestors (which I know I have seen all by then)

  3. Iterate through nodes_df again in reverse order to get node_descendants, a {node_id: set(ancestors)} dict adding for each node's descendants itself and its child's descendants (which I know I have seen all by then)


Code:



import pandas as pd
from collections import defaultdict

# phase 1
nodes_df = pd.read_csv('input.csv')

# phase 2
node_ancestors = defaultdict(set)
node_ancestors[0] = set([0])
for id, ndata in nodes_df1.iterrows():
node_ancestors[ndata['node_id']].add(ndata['node_id'])
node_ancestors[ndata['node_id']].update(node_ancestors[ndata['parent_id']])

# phase 3
node_descendants = defaultdict(set)
node_descendants[0] = set([0])
for id, ndata in nodes_df1[::-1].iterrows():
node_descendants[ndata['node_id']].add(ndata['node_id'])
node_descendants[ndata['parent_id']].
update(node_descendants[ndata['node_id']])



So, this takes minutes on my laptop, which is ages for my application. How do I improve?



Plausible directions:




  1. Can I use pandas better? Can I get node_ancestors and/or node_descendants by some clever join which is out of my league?

  2. Can I use a python graph library like Networkx or igraph (which in my experience is faster on large graphs)? E.g. in both libraries I have a get_all_shortest_paths methods, which returns something like a {node_id: dist} dictionary, from which I could select the keys, but... I need this for every node, so again a long long loop

  3. Parallelizing - no idea how to do this









share







New contributor




Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.







$endgroup$

















    0












    $begingroup$


    I have a CSV holding a "flat" table of a tree's edges (NOT binary, but a node cannot have two parents), ~1M edges:



    node_id parent_id
    1 0
    2 1
    3 1
    4 2
    ...


    The nodes are sorted in a way that a parent_id must always come before any of its children, so a parent_id will always be lower than node_id.



    I wish, for each node_id, to get the set of all ancestor nodes (including itself, propagated until root which is node 0 here), and a set of all descendant nodes (including itself, propagated until leaves), and speed is crucial.



    Currently what I do at high level:




    1. Read the CSV in pandas, call it nodes_df

    2. Iterate once through nodes_df to get node_ancestors, a {node_id: set(ancestors)} dict adding for each node's ancestors itself and its parent's ancestors (which I know I have seen all by then)

    3. Iterate through nodes_df again in reverse order to get node_descendants, a {node_id: set(ancestors)} dict adding for each node's descendants itself and its child's descendants (which I know I have seen all by then)


    Code:



    import pandas as pd
    from collections import defaultdict

    # phase 1
    nodes_df = pd.read_csv('input.csv')

    # phase 2
    node_ancestors = defaultdict(set)
    node_ancestors[0] = set([0])
    for id, ndata in nodes_df1.iterrows():
    node_ancestors[ndata['node_id']].add(ndata['node_id'])
    node_ancestors[ndata['node_id']].update(node_ancestors[ndata['parent_id']])

    # phase 3
    node_descendants = defaultdict(set)
    node_descendants[0] = set([0])
    for id, ndata in nodes_df1[::-1].iterrows():
    node_descendants[ndata['node_id']].add(ndata['node_id'])
    node_descendants[ndata['parent_id']].
    update(node_descendants[ndata['node_id']])



    So, this takes minutes on my laptop, which is ages for my application. How do I improve?



    Plausible directions:




    1. Can I use pandas better? Can I get node_ancestors and/or node_descendants by some clever join which is out of my league?

    2. Can I use a python graph library like Networkx or igraph (which in my experience is faster on large graphs)? E.g. in both libraries I have a get_all_shortest_paths methods, which returns something like a {node_id: dist} dictionary, from which I could select the keys, but... I need this for every node, so again a long long loop

    3. Parallelizing - no idea how to do this









    share







    New contributor




    Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
    Check out our Code of Conduct.







    $endgroup$















      0












      0








      0





      $begingroup$


      I have a CSV holding a "flat" table of a tree's edges (NOT binary, but a node cannot have two parents), ~1M edges:



      node_id parent_id
      1 0
      2 1
      3 1
      4 2
      ...


      The nodes are sorted in a way that a parent_id must always come before any of its children, so a parent_id will always be lower than node_id.



      I wish, for each node_id, to get the set of all ancestor nodes (including itself, propagated until root which is node 0 here), and a set of all descendant nodes (including itself, propagated until leaves), and speed is crucial.



      Currently what I do at high level:




      1. Read the CSV in pandas, call it nodes_df

      2. Iterate once through nodes_df to get node_ancestors, a {node_id: set(ancestors)} dict adding for each node's ancestors itself and its parent's ancestors (which I know I have seen all by then)

      3. Iterate through nodes_df again in reverse order to get node_descendants, a {node_id: set(ancestors)} dict adding for each node's descendants itself and its child's descendants (which I know I have seen all by then)


      Code:



      import pandas as pd
      from collections import defaultdict

      # phase 1
      nodes_df = pd.read_csv('input.csv')

      # phase 2
      node_ancestors = defaultdict(set)
      node_ancestors[0] = set([0])
      for id, ndata in nodes_df1.iterrows():
      node_ancestors[ndata['node_id']].add(ndata['node_id'])
      node_ancestors[ndata['node_id']].update(node_ancestors[ndata['parent_id']])

      # phase 3
      node_descendants = defaultdict(set)
      node_descendants[0] = set([0])
      for id, ndata in nodes_df1[::-1].iterrows():
      node_descendants[ndata['node_id']].add(ndata['node_id'])
      node_descendants[ndata['parent_id']].
      update(node_descendants[ndata['node_id']])



      So, this takes minutes on my laptop, which is ages for my application. How do I improve?



      Plausible directions:




      1. Can I use pandas better? Can I get node_ancestors and/or node_descendants by some clever join which is out of my league?

      2. Can I use a python graph library like Networkx or igraph (which in my experience is faster on large graphs)? E.g. in both libraries I have a get_all_shortest_paths methods, which returns something like a {node_id: dist} dictionary, from which I could select the keys, but... I need this for every node, so again a long long loop

      3. Parallelizing - no idea how to do this









      share







      New contributor




      Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.







      $endgroup$




      I have a CSV holding a "flat" table of a tree's edges (NOT binary, but a node cannot have two parents), ~1M edges:



      node_id parent_id
      1 0
      2 1
      3 1
      4 2
      ...


      The nodes are sorted in a way that a parent_id must always come before any of its children, so a parent_id will always be lower than node_id.



      I wish, for each node_id, to get the set of all ancestor nodes (including itself, propagated until root which is node 0 here), and a set of all descendant nodes (including itself, propagated until leaves), and speed is crucial.



      Currently what I do at high level:




      1. Read the CSV in pandas, call it nodes_df

      2. Iterate once through nodes_df to get node_ancestors, a {node_id: set(ancestors)} dict adding for each node's ancestors itself and its parent's ancestors (which I know I have seen all by then)

      3. Iterate through nodes_df again in reverse order to get node_descendants, a {node_id: set(ancestors)} dict adding for each node's descendants itself and its child's descendants (which I know I have seen all by then)


      Code:



      import pandas as pd
      from collections import defaultdict

      # phase 1
      nodes_df = pd.read_csv('input.csv')

      # phase 2
      node_ancestors = defaultdict(set)
      node_ancestors[0] = set([0])
      for id, ndata in nodes_df1.iterrows():
      node_ancestors[ndata['node_id']].add(ndata['node_id'])
      node_ancestors[ndata['node_id']].update(node_ancestors[ndata['parent_id']])

      # phase 3
      node_descendants = defaultdict(set)
      node_descendants[0] = set([0])
      for id, ndata in nodes_df1[::-1].iterrows():
      node_descendants[ndata['node_id']].add(ndata['node_id'])
      node_descendants[ndata['parent_id']].
      update(node_descendants[ndata['node_id']])



      So, this takes minutes on my laptop, which is ages for my application. How do I improve?



      Plausible directions:




      1. Can I use pandas better? Can I get node_ancestors and/or node_descendants by some clever join which is out of my league?

      2. Can I use a python graph library like Networkx or igraph (which in my experience is faster on large graphs)? E.g. in both libraries I have a get_all_shortest_paths methods, which returns something like a {node_id: dist} dictionary, from which I could select the keys, but... I need this for every node, so again a long long loop

      3. Parallelizing - no idea how to do this







      python tree graph pandas





      share







      New contributor




      Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.










      share







      New contributor




      Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.








      share



      share






      New contributor




      Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.









      asked 4 mins ago









      Giora SimchoniGiora Simchoni

      1011




      1011




      New contributor




      Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.





      New contributor





      Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.






      Giora Simchoni is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.






















          0






          active

          oldest

          votes











          Your Answer





          StackExchange.ifUsing("editor", function () {
          return StackExchange.using("mathjaxEditing", function () {
          StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
          StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["\$", "\$"]]);
          });
          });
          }, "mathjax-editing");

          StackExchange.ifUsing("editor", function () {
          StackExchange.using("externalEditor", function () {
          StackExchange.using("snippets", function () {
          StackExchange.snippets.init();
          });
          });
          }, "code-snippets");

          StackExchange.ready(function() {
          var channelOptions = {
          tags: "".split(" "),
          id: "196"
          };
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function() {
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled) {
          StackExchange.using("snippets", function() {
          createEditor();
          });
          }
          else {
          createEditor();
          }
          });

          function createEditor() {
          StackExchange.prepareEditor({
          heartbeatType: 'answer',
          autoActivateHeartbeat: false,
          convertImagesToLinks: false,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: null,
          bindNavPrevention: true,
          postfix: "",
          imageUploader: {
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          },
          onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          });


          }
          });






          Giora Simchoni is a new contributor. Be nice, and check out our Code of Conduct.










          draft saved

          draft discarded


















          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fcodereview.stackexchange.com%2fquestions%2f211989%2fpython-get-all-node-descendants-in-a-tree%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown

























          0






          active

          oldest

          votes








          0






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes








          Giora Simchoni is a new contributor. Be nice, and check out our Code of Conduct.










          draft saved

          draft discarded


















          Giora Simchoni is a new contributor. Be nice, and check out our Code of Conduct.













          Giora Simchoni is a new contributor. Be nice, and check out our Code of Conduct.












          Giora Simchoni is a new contributor. Be nice, and check out our Code of Conduct.
















          Thanks for contributing an answer to Code Review Stack Exchange!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          Use MathJax to format equations. MathJax reference.


          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fcodereview.stackexchange.com%2fquestions%2f211989%2fpython-get-all-node-descendants-in-a-tree%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          404 Error Contact Form 7 ajax form submitting

          How to know if a Active Directory user can login interactively

          Refactoring coordinates for Minecraft Pi buildings written in Python