Pagini recente » Cod sursa (job #1726712) | Cod sursa (job #2481498) | Cod sursa (job #1279682) | Cod sursa (job #1424495) | Cod sursa (job #237952)
Cod sursa(job #237952)
#include <stdio.h>
FILE*in = fopen("pascal.in","r");
FILE*out= fopen("pascal.out","w");
unsigned long cmmdc(unsigned long a, unsigned long b){
unsigned long r;
while(b!=0){
r=a%b;
a=b;
b=r;
}
return a;
}
int factj=1,c=0,n;
unsigned long rez,numitor;
unsigned long numarator[2000],k,l,i,j;
int main(){
fscanf(in,"%ld%d",&i,&factj);
if((i+1)%2 == 0) n = (i+1)/2;
else n = (i+1)/2+1;
for(j=1;j<n;j++){
// formez fractia
// numitor
if(j!=0) factj *= j;
numitor = factj;
// numarator
for(k=1;k<=j;k++)
numarator[k] = i-j+k;
// operatiile de simplificare
for(l=1;l<k;l++){
rez = cmmdc(numarator[l],numitor);
if(rez == numitor){ c++; break;}
else
if(rez == numarator[l]){
numarator[l] = 1;
numitor /= rez;
}else{
numitor /= rez;
numarator[l] /= rez;
}
}
}
if((i+1)%2 == 0) c = c * 2;
else c = c * 2 -1;
fprintf(out,"%d",c);
return 0;
}