We define and compute many-body topological invariants of interacting fermionic symmetryprotected topological phases, protected by an orientation-reversing symmetry, such as time-reversal or reflection symmetry. The topological invariants are given by partition functions obtained by a path integral on unoriented spacetime which, as we show, can be computed for a given ground state wave function by considering a non-local operation, "partial" reflection or transpose. As an application of our scheme, we study the Z8 and Z16 classification of topological superconductors in one and three dimensions.The Thouless-Kohmoto-Nightingale-den Nijs (TKNN) formula [1, 2] is the prototype for topological characterization of phases of matter. It relates the quantized Hall conductance to the (first) Chern number defined for Bloch wave functions. At the level of many-body physics, the quantized Hall conductance can also be formulated in terms of ground state wave functions in the presence of twisted boundary conditions ("the many-body Chern number") [3]. In contrast to local order parameters, the TKNN integer distinguishes different quantum phases of matter by focusing on their global topological properties.More recently, the discovery of topological insulators and superconductors [4, 5] has led to a new research frontier, generally referred to as symmetry protected topological (SPT) phases. These phases are adiabatically connected to topologically trivial states, i.e., atomic insulators which can be represented as simple product states without any entanglement. Nevertheless, they are topologically distinct once a symmetry condition, e.g., timereversal symmetry, is imposed. A complete classification of the noninteracting fermionic SPT phases protected by non-spatial discrete symmetries [6][7][8], as well as crystalline SPT phases protected by spatial symmetries [9][10][11][12], were achieved.However, it was later discovered that the noninteracting topological classification is not the full story, and can be dramatically altered once interaction effects are taken into account [13]. Since then, there have been several works which discuss the breakdown of the noninteracting classification in the presence of interactions [14][15][16][17][18][19][20][21][22][23][24][25].There are various topological invariants for noninteracting fermionic SPT phases using single-particle states (e.g., Bloch wave functions). For example, the Z 2 -valued topological invariants have been introduced for topological insulators both in two and three spatial dimensions [26][27][28][29]. For topological superconductors protected by time-reversal symmetry, the integer-valued topological invariants ("the winding number") have been introduced [6]. However, the discovered breakdown of non-interacting classification clearly indicates that the situation at the interacting level is more intricate, and a general framework to distinguish interacting fermionic SPT phases is lacking. This should be contrasted from the quantized Hall conductance, which can be formulate...