Pagini recente » Cod sursa (job #664865) | Rating Boiculese Claudiu (yusty95sv) | Cod sursa (job #2132682) | Cod sursa (job #1691183) | Cod sursa (job #1065043)
#include <fstream>
#include <cstring>
using namespace std;
int Fact[100005],Power2[100005];
int N;
void mulLess(int A[],int B)
{
int i, t = 0;
for (i = 1; i <= A[0] || t; i++, t /= 10000)
A[i] = (t += A[i] * B) % 10000;
A[0] = i - 1;
}
void mulLarge(int A[], int B[])
{
int i, j, t, C[100005];
memset(C, 0, sizeof(C));
for (i = 1; i <= A[0]; i++)
{
for (t=0, j=1; j <= B[0] || t; j++, t/=10000)
C[i+j-1]=(t+=C[i+j-1]+A[i]*B[j])%10000;
if (i + j - 2 > C[0]) C[0] = i + j - 2;
}
memcpy(A, C, sizeof(C));
}
void Solve()
{
int i;
freopen ("patrate2.out", "w", stdout);
Fact[0]=Fact[1]=1;
for(i=2;i<=N;i++)
mulLess(Fact,i);
Power2[0]=Power2[1]=1;
for(i=1;i<=N*N;i++)
mulLess(Power2,2);
mulLarge(Fact,Power2);
printf("%d", Fact[Fact[0]]);
for(i=Fact[0]-1;i>0;--i)
printf("%04d", Fact[i]);
printf("\n");
}
int main()
{
freopen ("patrate2.in","r",stdin);
scanf("%d", &N);
Solve();
return 0;
}