Cod sursa(job #1104132)

Utilizator pop_bogdanBogdan Pop pop_bogdan Data 10 februarie 2014 14:50:49
Problema Patrate2 Scor 70
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.08 kb
#include <fstream>
#include <cstring>
using namespace std;

ifstream is("patrate2.in");
ofstream os("patrate2.out");

int S[5001];
int S2[5001];

void Multiply(int A[],int B);
void Multiply2(int A[], int B[]);

int main()
{
    int n;
    is >> n;
    S[0] = 1;
    S[1] = 1;
    for ( int i = 1; i <= n; ++i )
        Multiply(S,i);
    S2[0] = 1;
    S2[1] = 2;
    for ( int i = 1; i <= n*n-1; ++i )
        Multiply(S2,2);
    Multiply2(S,S2);
    for ( int i = S[0]; i >= 1; --i )
        os << S[i];
    return 0;
}

void Multiply(int A[], int B)
{
    int i, t = 0;
    for (i = 1; i <= A[0] || t; i++, t /= 10)
        A[i] = (t += A[i] * B) % 10;
    A[0] = i - 1;
}

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