float p=1,s=0; int n,k; for(k=1;k<=21;k=k+4) { for(n=1;n<=k;n++) p=p*n; s=s+p; }
for(k=1;k<=21;k=k+4)