Pagini recente » Cod sursa (job #805650) | Cod sursa (job #533798) | Cod sursa (job #1222174) | Cod sursa (job #1369260) | Cod sursa (job #273382)
Cod sursa(job #273382)
#include<stdio.h>
#include<string.h>
using namespace std;
int n,N,i,sol[1111],doi[1111];
void mul(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 mul2(int A[], int B[])
{
int i, j, t, C[1111];
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));
}
int main(){
FILE *f = fopen("patrate2.in","r");
FILE *g = fopen("patrate2.out","w");
fscanf(f,"%d",&n);
N = n*n;
sol[0] = 1; sol[1] = 1;
doi[0] = 1; doi[1] = 2;
//for(i=2; i<=N; i++)
//mul(sol,2);
while(N){
if(N&1 == 1)
//rez=(rez*n) % MOD;
mul2(sol,doi);
N = N>>1;
//n=(n*n) % MOD;
mul2(doi,doi);
}
for(i=2; i<=n; i++)
mul(sol,i);
for(i=sol[0]; i>=1; i--)
fprintf(g,"%d",sol[i]);
fclose(f);
fclose(g);
return 0;
}