Cod sursa(job #237952)

Utilizator razyelxrazyelx razyelx Data 30 decembrie 2008 23:16:23
Problema Pascal Scor 20
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.1 kb
#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;
}