Cod sursa(job #147776)
#include <stdio.h>
#define DIM 5000001
int v[DIM];
int n,d,i,j,k,nr,dd;
int main(){
FILE *f = fopen("pascal.in","r");
fscanf(f,"%d %d",&n,&d);
fclose(f);
if ((d!=2)&&(d!=5)&&(d!=3))
exit(0);
for (i=1;i<=n;i++) {
j=i;
k=0;
while (j%d==0) {
k++;
j/=d;
}
v[i]=k+v[i-1];
}
for (i=0;i<=n/2-1;i++) {
//(n!/(i!*(n-i)!))
dd = v[n]-v[i]-v[n-i];
if (dd>0) nr+=2;
}
if (n%2==0) {
i=n/2;
dd = v[n]-v[i]-v[n-i];
if (dd>0) nr+=1;
}
FILE *g = fopen("pascal.out","w");
fprintf(g,"%d",nr);
fclose(g);
return 0;
}