Então eu fiz isso de uma maneira estranha. Percebi que havia dois padrões na forma como a matriz se forma.
O primeiro é como o padrão das linhas superiores tem a diferença entre cada termo, aumentando de 1 -> h, onde h é a altura e l é o comprimento. Então eu construo a linha superior com base nesse padrão
Para uma matriz de dim (3,4) dando a, max RoC = 3
vamos ver a linha superior do formulário
1, (1+1), (2+2), (4+3) = 1, 2, 4, 7
Suponha, em vez disso, que o dim (3,9), dando um max RoC = 3
, verá uma linha superior de
`1, (1+1), (2+2), (4+3), (7+3), (10+3), (13+3), (16+3), (19+3) = 1, 2, 4, 7, 10, 13, 16, 19, 22
O segundo padrão é como as linhas mudam uma da outra. Se considerarmos a matriz:
1 2 4 7 11
3 5 8 12 16
6 9 13 17 20
10 14 18 21 23
15 19 22 24 25
e subtrair cada linha da linha abaixo (ignorando a linha extra), obtemos
2 3 4 5 5
3 4 5 5 4
4 5 5 4 3
5 5 4 3 2
Ao ver essa matriz, podemos notar que essa matriz é a sequência em 2 3 4 5 5 4 3 2
que, para cada linha, 5 termos desse padrão são deslocados por 1 para cada linha. Veja abaixo o visual.
|2 3 4 5 5| 4 3 2
2 |3 4 5 5 4| 3 2
2 3 |4 5 5 4 3| 2
2 3 4 |5 5 4 3 2|
Portanto, para obter a matriz final, pegamos nossa primeira linha que criamos e produzimos essa linha adicionada com os 5 termos necessários desse padrão.
Esse padrão sempre terá as características do início 2-> max value
e do fim, max value -> 2
onde max value = min(h+1, l)
o número de vezes que o valor máximo aparecerá é o appearances of max = h + l -2*c -2
localc = min(h+1, l) - 2
Então, em todo o meu método de criação de novas linhas parece
1 2 3 7 11 + |2 3 4 5 5|4 3 2 = 3 5 8 12 16
3 5 8 12 16 + 2|3 4 5 5 4|3 4 2 = 6 9 13 17 20
6 9 13 17 20 + 2 3|4 5 5 4 3|4 2 = 10 14 18 21 23
10 14 18 21 23 + 2 3 4|5 5 4 3 2| = 15 19 22 24 25
Código relevante abaixo. Não acabou sendo curto, mas ainda gosto do método.
o,r=len,range
def m(l,h):
a,t=[1+sum(([0]+[x for x in r(1,h)]+[h]*(l-h))[:x+1]) for x in r(l)],min(l,h+1);s,c=[x for x in r(2,t)],[a[:]]
for i in r(h-1):
for j in r(o(a)):
a[j]+=(s+[t]*(l+h-2*(t-2)-2)+s[::-1])[0+i:l+i][j]
c+=[a[:]]
for l in c:print(l)
Experimente online!