щас что-нибудь попробую сообразить, самому интересно стало

я попытался решить задачу так:
import numpy as np
n=4
psum= n*(n^2+1)//2
def coun_iter(a, k, csum, p_sum,n):# sum[i=1..n](c[n-i]*i)=psum
   #print(csum)
   if (k<n-1):
       for a[k] in np.arange(1,(p_sum-csum) // (n-k)):
           coun_iter(a,k+1,csum+(n-k)*a[k],p_sum,n)
   else:

       a[k]=p_sum-csum
       print(a)
c=np.ones(n, np.uint16)
coun_iter(c,0,0,psum,n)

но оно мне выдает:
[1 1 1 5]
[1 1 2 3]
[1 2 1 2]

а как из условия задачи следует "то есть, n*(n^2+1)/2"?

я вслед за многими математиками насилую проблему нумерации магических квадратов

но это выглядит как маловато решений.

покажите решение

облом.

я просто забыл, что питон в степень возводит по-другому.

и решений было меньше потому что считало неправильно  n*(n^2+1)//2. надо было  n*(n*n+1)//2

Но я проверил на n=3 и оно по-прежнему не находит некоторые решения.

Здравствуйте.  помогите 2733

import numpy as np
n=4
psum= n*(n*n+1)//2
def coun_iter(a, k, csum, p_sum,n):# sum[i=1..n](c[n-i]*i)=psum
   if (k<n-1):
       for a[k] in np.arange(1,1+(p_sum-csum) // (n-k)):
           coun_iter(a,k+1,csum+(n-k)*a[k],p_sum,n)
   else:
        a[k]=p_sum-csum
        if np.cumsum(a)[n-1] < (n * n+1):
            print(a,np.cumsum(a),np.sum(np.cumsum(a)))
c=np.ones(n, np.uint16)
print(psum)
coun_iter(c,0,0,psum,n)

вроде бы теперь работает.

Подскажите пожалуйста, количество фундаментальных решений всегда будет равно количеству свободных неизвестных?

Кто нибудь может пояснить вкратце за рассеяние Ми в фотонных кристаллах ?

Так там сумма двух геометрических прогрессий

Эх, мне бы кто-нибудь за симметрии параметра порядка в купратах пояснил...

Спасибо

🤔🤔