Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- siZe=28;
- a181936_list(DIM)={
- my(E);
- E=matrix(DIM,DIM);
- E[1,1]=1;
- for(n=2,DIM,
- if(!((n-1)%5),
- E[n,1]=0;
- forstep(k=n-1,1,-1,
- E[k,n-k+1]= ( E[k+1,n-k] + E[k,n-k] )
- ),
- E[1,n]=0;
- for(k=2,n,
- E[k,n-k+1]= ( E[k-1,n-k+2] + E[k-1,n-k+1] )
- )
- )
- );
- ans=vector(DIM,k,(k==1)*E[1,1]+(k!=1)*(E[k,1]+E[1,k]));
- };
- table=a181936_list(siZe);
- for(j=1,#table,print((j-1)," ",table[j]));
- quit();
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement