Combc; ll@n,@m--,@k; Mint z; rep(i,k){ Mint j=i+1; if(n==1){ z=c.C(k,j)*j*(j-1)**m-z; } if(n==2){ z=c.C(k,j)*j*(j-1)*((j-3)*j+3)**m-z; } if(n==3){ Matrixa(2,2); a[0][0]=a[0][1]=((j-6)*j+13)*j-10; a[0][0]+=j-3; a[1][0]=a[1][1]=(j-4)*j+5; a[1][1]+=j-2; a**=m; z=c.C(k,j)*j*(j-1)*((j-2)*(a[0][0]+a[1][0])+(a[0][1]+a[1][1]))-z; } } wt(z);