Cod sursa(job #437652)

Utilizator AndreiDDiaconeasa Andrei AndreiD Data 9 aprilie 2010 23:52:47
Problema Patrate2 Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.14 kb
#include <cstdio>
#include <cstring>

#define file_in "patrate2.in"
#define file_out "patrate2.out"

#define Base 1000

int n,rez1[10000],rez2[10100];

void init(int cif, int rez[])
{
    int i,t=0;
    for(i=1;i<=rez[0]||t;++i)
    {
        t+=rez[i]*cif;
        rez[i]=t%Base;
        t/=Base;
    }
    rez[0]=i-1;
}

void mul(int A[], int B[])
{
      int i, j, t, C[10000];
      memset(C, 0, sizeof(C));
      for (i = 1; i <= A[0]; i++)
      {
              for (t=0, j=1; j <= B[0] || t; j++, t/=Base)
                      C[i+j-1]=(t+=C[i+j-1]+A[i]*B[j])%Base;
              if (i + j - 2 > C[0]) C[0] = i + j - 2;
      }
      memcpy(A, C, sizeof(C));
}


int main()
{
    int i;
    freopen(file_in,"r",stdin);
    freopen(file_out,"w",stdout);

    scanf("%d", &n);

    rez1[0]=rez1[1]=1;

    //init(n,rez1);
    for (i=1;i<=n*n;++i)
         init(2,rez1);
	rez2[0]=rez2[1]=1;

	for (i=1;i<=n;++i)
		init(i,rez2);
	mul(rez1,rez2);
    printf("%d",rez1[rez1[0]]);
    for(i=rez1[0]-1;i>0;--i)
        printf("%03d",rez1[i]);


    fclose(stdin);
    fclose(stdout);

    return 0;

}