Cod sursa(job #123974)

Utilizator FlorianFlorian Marcu Florian Data 17 ianuarie 2008 20:15:46
Problema Patrate2 Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.73 kb
#include<stdio.h>
#include<string.h>
#define B 1000000
#define nr 6
FILE*f=fopen("patrate2.in","r");
FILE*g=fopen("patrate2.out","w");
long n,a[5000],p[5000],z[5000];
void inmult(long c[],long k)
        {
        long i,t,x;
        t=0; i=1;
        while(i<=c[0]||t>0)
                {
                t+=c[i]*k;
                c[i]=t%B;
                t=t/B;
                i++;
                }
        c[0]=i-1;
        }

void inmultire(long a[],long b[])
        {
        long j,i,t=0,c;
        for(i=1;i<=a[0];++i)
                for(j=1;j<=b[0];++j)
                        z[i+j-1]+=a[i]*b[j];
        z[0]=a[0]+b[0]-1;
        for(i=1;i<=z[0];++i)
                {
                c=z[i];
                z[i]=(z[i]+t)%B;
                t=(c+t)/B;
                }
        if(t>0) z[++z[0]]=t;
        }
void afisare()
        {
         long x,pr,j;
         fprintf(g,"%ld",z[z[0]]);
        for(long i=z[0]-1;i>=1;--i)
            {
             for(x=z[i],pr=0;x;++pr,x=x/10);
             for(j=1;j<=nr-pr;++j) fprintf(g,"0");

             fprintf(g,"%ld",z[i]);

            }
        fprintf(g,"\n");
        }
void afis()
        {
         long x,c,i;

        for (i=z[0]; i>=1; i--)
     {
       c=0; x=z[i];
if (z[i])
{
       while (x) {c++; x/=10;}
       for (x=6; x>c; x--) if (i!=z[0]) fprintf(g,"0");
       fprintf(g,"%ld",z[i]);
}
else fprintf(g,"000000");
     }
}
int main()
        {
        //fscanf(f,"%ld",&n);
        n=2;
        long i;
        p[0]=1; p[1]=1;
        for(i=1;i<=n*n;++i)  {inmult(p,2);}
        a[0]=1; a[1]=1;
        for(i=1;i<=n;++i) inmult(a,i);
        inmultire(p,a);
        afis();
        return 0;
        }