*****************************************************************************
*       www.FindStat.org - The Combinatorial Statistic Finder               *
*                                                                           *
*       Copyright (C) 2019 The FindStatCrew <info@findstat.org>             *
*                                                                           *
*    This information is distributed in the hope that it will be useful,    *
*    but WITHOUT ANY WARRANTY; without even the implied warranty of         *
*    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.                   *
*****************************************************************************

-----------------------------------------------------------------------------
Statistic identifier: St001140

-----------------------------------------------------------------------------
Collection: Dyck paths

-----------------------------------------------------------------------------
Description: Number of indecomposable modules with projective and injective dimension at least two in the corresponding Nakayama algebra.

-----------------------------------------------------------------------------
References: [1]   Marczinzik, René Upper bounds for the dominant dimension of Nakayama and related algebras. [[zbMATH:06820683]]

-----------------------------------------------------------------------------
Code:
gap('LoadPackage("QPA");')
def NthRadical(M, n):
    if n == 0:
        f = gap.IdentityMapping(M)
    else:
        f = gap.RadicalOfModuleInclusion(M)
        N = gap.Source(f)
        for i in range(n-1):
            h = gap.RadicalOfModuleInclusion(N);
            N = gap.Source(h)
            f = h * f
    return f

def ARQuiverNak(A):
    injA = gap.IndecInjectiveModules(A)
    L = [gap.Source(NthRadical(inj, j))
         for inj in injA
         for j in range(gap.Dimension(inj).sage())]
    return L

def kupisch(D):
    H = D.heights()
    return [1+H[i] for i, s in enumerate(D) if s == 0]+[1]

def statistic(D):
    K = kupisch(D)
    A = gap.NakayamaAlgebra(K, gap.GF(3))
    L = ARQuiverNak(A)
    return sum(1 for x in L
               if gap.ProjDimensionOfModule(x, 30) >= 2 and gap.InjDimensionOfModule(x, 30) >= 2)


-----------------------------------------------------------------------------
Statistic values:

[1,0]                     => 0
[1,0,1,0]                 => 0
[1,1,0,0]                 => 0
[1,0,1,0,1,0]             => 0
[1,0,1,1,0,0]             => 0
[1,1,0,0,1,0]             => 0
[1,1,0,1,0,0]             => 0
[1,1,1,0,0,0]             => 0
[1,0,1,0,1,0,1,0]         => 1
[1,0,1,0,1,1,0,0]         => 0
[1,0,1,1,0,0,1,0]         => 1
[1,0,1,1,0,1,0,0]         => 0
[1,0,1,1,1,0,0,0]         => 0
[1,1,0,0,1,0,1,0]         => 0
[1,1,0,0,1,1,0,0]         => 0
[1,1,0,1,0,0,1,0]         => 0
[1,1,0,1,0,1,0,0]         => 0
[1,1,0,1,1,0,0,0]         => 0
[1,1,1,0,0,0,1,0]         => 0
[1,1,1,0,0,1,0,0]         => 0
[1,1,1,0,1,0,0,0]         => 0
[1,1,1,1,0,0,0,0]         => 0
[1,0,1,0,1,0,1,0,1,0]     => 2
[1,0,1,0,1,0,1,1,0,0]     => 1
[1,0,1,0,1,1,0,0,1,0]     => 1
[1,0,1,0,1,1,0,1,0,0]     => 1
[1,0,1,0,1,1,1,0,0,0]     => 0
[1,0,1,1,0,0,1,0,1,0]     => 1
[1,0,1,1,0,0,1,1,0,0]     => 1
[1,0,1,1,0,1,0,0,1,0]     => 1
[1,0,1,1,0,1,0,1,0,0]     => 2
[1,0,1,1,0,1,1,0,0,0]     => 0
[1,0,1,1,1,0,0,0,1,0]     => 1
[1,0,1,1,1,0,0,1,0,0]     => 2
[1,0,1,1,1,0,1,0,0,0]     => 0
[1,0,1,1,1,1,0,0,0,0]     => 0
[1,1,0,0,1,0,1,0,1,0]     => 1
[1,1,0,0,1,0,1,1,0,0]     => 0
[1,1,0,0,1,1,0,0,1,0]     => 1
[1,1,0,0,1,1,0,1,0,0]     => 0
[1,1,0,0,1,1,1,0,0,0]     => 0
[1,1,0,1,0,0,1,0,1,0]     => 1
[1,1,0,1,0,0,1,1,0,0]     => 0
[1,1,0,1,0,1,0,0,1,0]     => 2
[1,1,0,1,0,1,0,1,0,0]     => 1
[1,1,0,1,0,1,1,0,0,0]     => 0
[1,1,0,1,1,0,0,0,1,0]     => 2
[1,1,0,1,1,0,0,1,0,0]     => 1
[1,1,0,1,1,0,1,0,0,0]     => 0
[1,1,0,1,1,1,0,0,0,0]     => 0
[1,1,1,0,0,0,1,0,1,0]     => 0
[1,1,1,0,0,0,1,1,0,0]     => 0
[1,1,1,0,0,1,0,0,1,0]     => 0
[1,1,1,0,0,1,0,1,0,0]     => 0
[1,1,1,0,0,1,1,0,0,0]     => 0
[1,1,1,0,1,0,0,0,1,0]     => 0
[1,1,1,0,1,0,0,1,0,0]     => 0
[1,1,1,0,1,0,1,0,0,0]     => 0
[1,1,1,0,1,1,0,0,0,0]     => 0
[1,1,1,1,0,0,0,0,1,0]     => 0
[1,1,1,1,0,0,0,1,0,0]     => 0
[1,1,1,1,0,0,1,0,0,0]     => 0
[1,1,1,1,0,1,0,0,0,0]     => 0
[1,1,1,1,1,0,0,0,0,0]     => 0
[1,0,1,0,1,0,1,0,1,0,1,0] => 3
[1,0,1,0,1,0,1,0,1,1,0,0] => 2
[1,0,1,0,1,0,1,1,0,0,1,0] => 2
[1,0,1,0,1,0,1,1,0,1,0,0] => 2
[1,0,1,0,1,0,1,1,1,0,0,0] => 1
[1,0,1,0,1,1,0,0,1,0,1,0] => 1
[1,0,1,0,1,1,0,0,1,1,0,0] => 1
[1,0,1,0,1,1,0,1,0,0,1,0] => 2
[1,0,1,0,1,1,0,1,0,1,0,0] => 3
[1,0,1,0,1,1,0,1,1,0,0,0] => 1
[1,0,1,0,1,1,1,0,0,0,1,0] => 1
[1,0,1,0,1,1,1,0,0,1,0,0] => 2
[1,0,1,0,1,1,1,0,1,0,0,0] => 1
[1,0,1,0,1,1,1,1,0,0,0,0] => 0
[1,0,1,1,0,0,1,0,1,0,1,0] => 2
[1,0,1,1,0,0,1,0,1,1,0,0] => 1
[1,0,1,1,0,0,1,1,0,0,1,0] => 2
[1,0,1,1,0,0,1,1,0,1,0,0] => 1
[1,0,1,1,0,0,1,1,1,0,0,0] => 1
[1,0,1,1,0,1,0,0,1,0,1,0] => 2
[1,0,1,1,0,1,0,0,1,1,0,0] => 1
[1,0,1,1,0,1,0,1,0,0,1,0] => 4
[1,0,1,1,0,1,0,1,0,1,0,0] => 3
[1,0,1,1,0,1,0,1,1,0,0,0] => 2
[1,0,1,1,0,1,1,0,0,0,1,0] => 2
[1,0,1,1,0,1,1,0,0,1,0,0] => 2
[1,0,1,1,0,1,1,0,1,0,0,0] => 2
[1,0,1,1,0,1,1,1,0,0,0,0] => 0
[1,0,1,1,1,0,0,0,1,0,1,0] => 1
[1,0,1,1,1,0,0,0,1,1,0,0] => 1
[1,0,1,1,1,0,0,1,0,0,1,0] => 2
[1,0,1,1,1,0,0,1,0,1,0,0] => 2
[1,0,1,1,1,0,0,1,1,0,0,0] => 2
[1,0,1,1,1,0,1,0,0,0,1,0] => 1
[1,0,1,1,1,0,1,0,0,1,0,0] => 2
[1,0,1,1,1,0,1,0,1,0,0,0] => 3
[1,0,1,1,1,0,1,1,0,0,0,0] => 0
[1,0,1,1,1,1,0,0,0,0,1,0] => 1
[1,0,1,1,1,1,0,0,0,1,0,0] => 2
[1,0,1,1,1,1,0,0,1,0,0,0] => 3
[1,0,1,1,1,1,0,1,0,0,0,0] => 0
[1,0,1,1,1,1,1,0,0,0,0,0] => 0
[1,1,0,0,1,0,1,0,1,0,1,0] => 2
[1,1,0,0,1,0,1,0,1,1,0,0] => 1
[1,1,0,0,1,0,1,1,0,0,1,0] => 1
[1,1,0,0,1,0,1,1,0,1,0,0] => 1
[1,1,0,0,1,0,1,1,1,0,0,0] => 0
[1,1,0,0,1,1,0,0,1,0,1,0] => 1
[1,1,0,0,1,1,0,0,1,1,0,0] => 1
[1,1,0,0,1,1,0,1,0,0,1,0] => 1
[1,1,0,0,1,1,0,1,0,1,0,0] => 2
[1,1,0,0,1,1,0,1,1,0,0,0] => 0
[1,1,0,0,1,1,1,0,0,0,1,0] => 1
[1,1,0,0,1,1,1,0,0,1,0,0] => 2
[1,1,0,0,1,1,1,0,1,0,0,0] => 0
[1,1,0,0,1,1,1,1,0,0,0,0] => 0
[1,1,0,1,0,0,1,0,1,0,1,0] => 2
[1,1,0,1,0,0,1,0,1,1,0,0] => 1
[1,1,0,1,0,0,1,1,0,0,1,0] => 1
[1,1,0,1,0,0,1,1,0,1,0,0] => 1
[1,1,0,1,0,0,1,1,1,0,0,0] => 0
[1,1,0,1,0,1,0,0,1,0,1,0] => 3
[1,1,0,1,0,1,0,0,1,1,0,0] => 2
[1,1,0,1,0,1,0,1,0,0,1,0] => 3
[1,1,0,1,0,1,0,1,0,1,0,0] => 3
[1,1,0,1,0,1,0,1,1,0,0,0] => 1
[1,1,0,1,0,1,1,0,0,0,1,0] => 2
[1,1,0,1,0,1,1,0,0,1,0,0] => 3
[1,1,0,1,0,1,1,0,1,0,0,0] => 1
[1,1,0,1,0,1,1,1,0,0,0,0] => 0
[1,1,0,1,1,0,0,0,1,0,1,0] => 2
[1,1,0,1,1,0,0,0,1,1,0,0] => 2
[1,1,0,1,1,0,0,1,0,0,1,0] => 2
[1,1,0,1,1,0,0,1,0,1,0,0] => 3
[1,1,0,1,1,0,0,1,1,0,0,0] => 1
[1,1,0,1,1,0,1,0,0,0,1,0] => 2
[1,1,0,1,1,0,1,0,0,1,0,0] => 4
[1,1,0,1,1,0,1,0,1,0,0,0] => 2
[1,1,0,1,1,0,1,1,0,0,0,0] => 0
[1,1,0,1,1,1,0,0,0,0,1,0] => 2
[1,1,0,1,1,1,0,0,0,1,0,0] => 4
[1,1,0,1,1,1,0,0,1,0,0,0] => 2
[1,1,0,1,1,1,0,1,0,0,0,0] => 0
[1,1,0,1,1,1,1,0,0,0,0,0] => 0
[1,1,1,0,0,0,1,0,1,0,1,0] => 1
[1,1,1,0,0,0,1,0,1,1,0,0] => 0
[1,1,1,0,0,0,1,1,0,0,1,0] => 1
[1,1,1,0,0,0,1,1,0,1,0,0] => 0
[1,1,1,0,0,0,1,1,1,0,0,0] => 0
[1,1,1,0,0,1,0,0,1,0,1,0] => 1
[1,1,1,0,0,1,0,0,1,1,0,0] => 0
[1,1,1,0,0,1,0,1,0,0,1,0] => 2
[1,1,1,0,0,1,0,1,0,1,0,0] => 1
[1,1,1,0,0,1,0,1,1,0,0,0] => 0
[1,1,1,0,0,1,1,0,0,0,1,0] => 2
[1,1,1,0,0,1,1,0,0,1,0,0] => 1
[1,1,1,0,0,1,1,0,1,0,0,0] => 0
[1,1,1,0,0,1,1,1,0,0,0,0] => 0
[1,1,1,0,1,0,0,0,1,0,1,0] => 1
[1,1,1,0,1,0,0,0,1,1,0,0] => 0
[1,1,1,0,1,0,0,1,0,0,1,0] => 2
[1,1,1,0,1,0,0,1,0,1,0,0] => 1
[1,1,1,0,1,0,0,1,1,0,0,0] => 0
[1,1,1,0,1,0,1,0,0,0,1,0] => 3
[1,1,1,0,1,0,1,0,0,1,0,0] => 2
[1,1,1,0,1,0,1,0,1,0,0,0] => 1
[1,1,1,0,1,0,1,1,0,0,0,0] => 0
[1,1,1,0,1,1,0,0,0,0,1,0] => 3
[1,1,1,0,1,1,0,0,0,1,0,0] => 2
[1,1,1,0,1,1,0,0,1,0,0,0] => 1
[1,1,1,0,1,1,0,1,0,0,0,0] => 0
[1,1,1,0,1,1,1,0,0,0,0,0] => 0
[1,1,1,1,0,0,0,0,1,0,1,0] => 0
[1,1,1,1,0,0,0,0,1,1,0,0] => 0
[1,1,1,1,0,0,0,1,0,0,1,0] => 0
[1,1,1,1,0,0,0,1,0,1,0,0] => 0
[1,1,1,1,0,0,0,1,1,0,0,0] => 0
[1,1,1,1,0,0,1,0,0,0,1,0] => 0
[1,1,1,1,0,0,1,0,0,1,0,0] => 0
[1,1,1,1,0,0,1,0,1,0,0,0] => 0
[1,1,1,1,0,0,1,1,0,0,0,0] => 0
[1,1,1,1,0,1,0,0,0,0,1,0] => 0
[1,1,1,1,0,1,0,0,0,1,0,0] => 0
[1,1,1,1,0,1,0,0,1,0,0,0] => 0
[1,1,1,1,0,1,0,1,0,0,0,0] => 0
[1,1,1,1,0,1,1,0,0,0,0,0] => 0
[1,1,1,1,1,0,0,0,0,0,1,0] => 0
[1,1,1,1,1,0,0,0,0,1,0,0] => 0
[1,1,1,1,1,0,0,0,1,0,0,0] => 0
[1,1,1,1,1,0,0,1,0,0,0,0] => 0
[1,1,1,1,1,0,1,0,0,0,0,0] => 0
[1,1,1,1,1,1,0,0,0,0,0,0] => 0

-----------------------------------------------------------------------------
Created: Apr 09, 2018 at 14:30 by Rene Marczinzik

-----------------------------------------------------------------------------
Last Updated: Aug 24, 2020 at 18:49 by Martin Rubey